source: codes/icosagcm/trunk/run_adv.def @ 210

Last change on this file since 210 was 25, checked in by dubos, 12 years ago

Minor changes :
caldyn_sw.f90, advect_tracer.f90
icosa_mod.f90 : added parameters for NCAR test cases needing global scope
guided_mod.f90 : CALL to guided_ncar now takes tt=it*dt instead of it as input

Significant changes :
timeloop_gcm.f90 : re-activated CALL to advection scheme
disvert_ncar.f90,
etat0_ncar.f90
guided_ncar_mod.f90 : simplification, introduced several getin(...), update due to recent changes in advection test cases (deformational flow, Hadley cell)
run_adv.def : new keys, reorganized for legibility

Tests :
icosa_gcm.exe tested with ncar_adv_shape=const and ncar_adv_wind=solid,deform,hadley.
q1=1 maintained to machine accuracy. Surface pressure slightly oscillates as expected.

FIXME : Tests by Sarvesh with revision 24 show incorrect advection of cosine bell by solid-body rotation. Not fixed.

File size: 2.2 KB
Line 
1#------------------------------- Mesh ---------------------------------
2
3# Number of subdivision on a main triangle (nbp) : integer (default=40)
4nbp=20
5
6# optim_it : mesh optimisation : number of iteration : integer (default=0)
7optim_it=0
8
9# Number of vertical layers (llm) : integer (default=19)
10llm=19
11
12# disvert : vertical discretisation : string (default='std') : std, ncar
13disvert=ncar
14
15#---------------------------------- Misc --------------------------------
16
17# number of tracer (nqtot) : integer (default 1)
18nqtot=1
19
20# timestep : real (default=90s)
21dt = 180.
22
23# number of timestep (default 100)
24itaumax = 38400
25
26# output field period : integer (default none)
27write_period=3600
28
29# etat0 : initial state : string (default=jablonowsky06) :
30# jablonowsky06, academic, ncar
31etat0=ncar
32
33# caldyn : computation type for gcm equation : string (default=gcm) : gcm, adv
34caldyn=adv
35
36#------------------------------ Dynamics --------------------------------
37
38# scheme type : string ( default='adam_bashforth')
39# euler, leapfrog, leapfrog_matsuno, adam_bashforth)
40scheme = euler
41
42# matsuno period : integer ( default=5)
43matsuno_period = 5
44
45# dissipation time graddiv : real (default=5000)
46tetagdiv = 1800
47tau_graddiv = 1800
48
49# dissipation time nxgradrot (default=5000)
50tetagrot=1800
51tau_gradrot=1800
52tau_divgrad=1800
53
54# guided_type : string (default=none) : none, ncar
55guided_type=ncar
56
57#-------------------- parameters for NCAR test cases ------------------------
58
59# NCAR advection test, initial tracer : string ( default='cos_bell')
60# const, slotted_cyl, cos_bell, dbl_cos_bell_q1, dbl_cos_bell_q2, complement, hadley
61ncar_adv_shape=complement
62
63# NCAR advection test, wind field : string (default='deform') : solid, deform, hadley
64ncar_adv_wind=deform
65
66# ncar_dz : model layer thickness in meters: real (default=400)
67# used if disvert=ncar
68ncar_dz=400
69
70# ncar_T0 : reference temperature for NCAR test cases : real (default=300)
71# also used by disvert if disvert=ncar
72ncar_T0=300
73
74# ncar_p0 : reference pressure for NCAR test cases : real (default=1e5)
75# also used by disvert if disvert=ncar
76ncar_p0=1e5
77
78# ncar_disvert_c : exponent for B(eta) : integer (default=1)
79# used by disvert if disvert=ncar
80ncar_disvert_c=1
Note: See TracBrowser for help on using the repository browser.