Changeset 150

Timestamp:

04/03/13 12:09:10 (11 years ago)

Author:

sdubey

Message:

Added few new routines to read NC files and compute diagnostics to r145.

Few routines of dry physics including radiation module, surface process and convective adjustment in new routine phyparam.f90. dynetat to read start files for dynamics.

check_conserve routine to compute conservation of quatities like mass, energy etc.etat0_heldsz.f90 for held-suarez test case initial conditions.

new Key time_style=lmd or dcmip to use day_step, ndays like in LMDZ

Location:

codes/icosagcm/trunk

Files:

• arch/arch-x86_skd.fcm (added)
• arch/arch-x86_skd.path (added)
• arch/arch-x86_skdmpi.fcm (added)
• arch/arch-x86_skdmpi.path (added)
• bld.cfg (modified) (1 diff)
• clean (added)
• compile (added)
• param_sets/dcmip2012/test4/test4.1/test4.1-0/run.def (modified) (7 diffs)

codes/icosagcm/trunk/bld.cfg

@@ -29 +29 @@
 bld::tool::fflags    %FFLAGS 
 bld::tool::fppkeys   %CPP_KEY %FPP_DEF
-bld::tool::fpp       %FPP
+#bld::tool::fpp       %FPP
 
 # Pre-process code before analysing dependencies

codes/icosagcm/trunk/param_sets/dcmip2012/test4/test4.1/test4.1-0/run.def

@@ -4 +4 @@
 
 # Number of subdivision on a main triangle (nbp) : integer (default=40)
-nbp=80
+nbp=40
 
 # nbp                 20  40  80 160
@@ -10 +10 @@
 
 # sub splitting of main rhombus : integer (default=1)
-nsplit_i=2
-nsplit_j=2
+nsplit_i=1
+nsplit_j=1
 
 # Number of vertical layer (llm) : integer (default=19)
@@ -25 +25 @@
 # optim_it : mesh optimisation : number of iteration : integer (default=0)
 #optim_it=1000
-optim_it=100
+optim_it=0
 
 # ---------------------------------- Time---------------------------------
 
 # scheme type : string ( default='runge_kutta') euler, leapfrog_matsuno, runge_kutta )
-scheme = runge_kutta
+scheme = runge_kutta   leapfrog_matsuno 
 
 # matsuno period : integer ( default=5)
@@ -36 +36 @@
 
 # timestep : real (default=480s)
-dt = 240.
+dt = 480.
 
 # advection called every itau_adv time steps : integer (default 2)
@@ -87 +87 @@
 
 # dissipation time graddiv : real (default=5000)
-tau_graddiv = 2500
+tau_graddiv = 5000
 
 # number of iteration for graddiv : integer (default=1)
@@ -93 +93 @@
 
 # dissipation time nxgradrot (default=5000)
-tau_gradrot = 2500
+tau_gradrot = 5000
 
 # number of iteration for nxgradrot : integer (default=1)
@@ -99 +99 @@
 
 # dissipation time divgrad (theta) (default=5000)
-tau_divgrad=2500
+tau_divgrad= 5000
 
 # number of iteration for divgrad : integer (default=1)