#!/usr/bin/perl
use IO::File;
#
# Make default namelist read file for tgcm24:
#
$namelist_file = "tgcm24_default.inp";
$fh_out = new IO::File;
$fh_out->open(">$namelist_file") or die ">>> Can't open new file $namelist_file\n";
print $fh_out <<"EOF";
&tgcm_input
;
; Namelist input file for tgcm24.
;
 LABEL='TIME-GCM tgcm24'
 DATE=2003,80
 CALDAY = 80
;
; -rw----r--  2 ROBLE   28103032  485238856 Jun 28  2004 py3n008.nc
; 03071  71  0  0 to 03080  80  0  0 by 1440 tgcm24 primary
;
; This startup file must be copied from $TCGMDATA (/hao/aim1/tgcm/data) to the 
;   execdir before execution (providing the mss or full path here does not work).
;
 SOURCE='ROBLE.tgcm24rh.py3n008_d80.nc'
 SOURCE_START = 80,0,0
;
 START=  80,0,0
 STOP=   81,0,0
 HIST=    1,0,0
 SAVE=    1,0,0
 STEP= 240
;
 OUTPUT = 'TGCM.tgcm24rh.py3n001.nc'
 MXHIST_PRIM =10 ! default max #histories per primary file = 5
;
 SECOUT = 'TGCM.tgcm24rh.sy3n001.nc'
 SECSTART = 80,1,0
 SECSTOP =  81,0,0
 SECHIST =   0,1,0
 SECSAVE =   0,1,0
 SECFLDS = 'TN','UN','VN','Z'
;SECFLDS = 'TN','UN','VN','W','O3','O1','OX','CO2','CO',
;   'NO','NO2','NOZ','H2O','OH','HO2','H','N4S','N2D',
;   'NE','HOX','Z'
 MXHIST_SECH = 48 ; default max #histories per sech file = 10
;
; OGSWM (APRIL)
 TIDE=6.4E+2,1.1E+2,2.6E+2,0.2E+2,0.8E+2,2.9,4.7,7.2,9.6,7.4
; NGSWM (APRIL)
;TIDE=4.3E+2,1.5E+2,1.9E+2,0.3E+2,0.8E+2,3.4,5.2,-4.6,-0.3,-4.3
; CASE A TIDAL PARAMETERS
TIDE=8.6E+2,3.0E+2,3.8E+2,0.6E+2,1.6E+2,3.4,5.2,-4.6,-0.3,-4.3
; CASE B TIDAL PARAMETERS
;TIDE=7.0E+2,2.0E+2,3.8E+2,0.3E+2,1.2E+2,3.4,5.2,-4.6,-0.3,-4.3
; CASE D TIDAL PARAMETERS
;TIDE=6.4E+2,1.5E+2,6.0E+2,0.3E+2,4.0E+2,3.4,5.2,-4.6,-0.3,-4.3
; CASE E TIDAL PARAMETERS
;TIDE=9.0E+2,1.5E+2,6.0E+2,0.3E+2,4.0E+2,3.4,5.2,-4.6,-0.3,-4.3
;
; CASE J TIDAL PARAMETERS
;TIDE=8.4E+2,1.5E+2,8.0E+2,0.3E+2,4.0E+2,3.4,5.2,-4.6,-0.3,-4.3
;
; BASE CASE TIDAL SIMULATION peqmv.. and seqmv.. run to day 20.0
;TIDE=6.4E+2,1.5E+2,3.0E+2,0.3E+2,4.0E+2,3.4,-1.3,7.2,9.6,7.4
;
; OGSWM (APRIL)
 TIDE2=14.1E+2,15.1
; NGSWM (APRIL)
;TIDE2=12.5E+2,11.4
; CASE J TIDAL PARAMETERS
 TIDE2=25.0E+2,16.4
;
 TIDE3M3=0.,0.
 MAG=-74.5,127.,79.,-70.
 DIFHOR=1
 IUIVI=1
;
; Set INIT_GLBM_AR =1 to initialize argon with global means
; Set INIT_GLBM_NAT=1 to initialize sodium with global means
;
;INIT_GLBM_AR  = 0
;INIT_GLBM_NAT = 0
 AURORA=1
; DISPOSE=0 -> do not dispose to mss.
; DISPOSE=1 -> dispose to mss during model execution.
; DISPOSE=2 -> dispose to mss after model execution.
 DISPOSE=0
 COLFAC=1.5
 F107 =  93.
 F107A = 126.
 POWER= 40.
 CTPOTEN=78.
; alfa30,e30,ed2 are not in tiegcm
 ALFA30= 40.
 E30=0.05
;E30=1.E-20
 ED2=1.E-20
 ALFAD2=10.
 BYIMF = 0.
/
EOF
print "$namelist_file\n";
exit;