Changeset 270 for codes/icosagcm/branches/SATURN_DYNAMICO/TEST

Timestamp:

09/05/14 14:42:27 (10 years ago)

Author:

millour

Message:

Added features for the Saturn case:

• Added possibility to run without startfi or restartfi.nc files
• Added reference temperature "temp_profile.txt" profile to start from
• More XIOS outputs, and put them on "presnivs (pressure) vertical coordinate
• Added "-openmp-threadprivate compat" OpenMP option in Ada arch file

EM

Location:

codes/icosagcm/branches/SATURN_DYNAMICO/TEST

Files:

• lmdz_physics.xml (modified) (3 diffs)
• run_icosa.def (modified) (4 diffs)
• temp_profile.txt (added)

codes/icosagcm/branches/SATURN_DYNAMICO/TEST/lmdz_physics.xml

@@ -13 +13 @@
      <field id="OLR"    long_name="outgoing longwave rad."      unit="W m-2"/>
      
-     <field_group axis_ref="altitude">
+     <field_group axis_ref="presnivs">
        <field id="temp" long_name="temperature"                 unit="K"/>
        <field id="u"    long_name="Zonal wind"                  unit="m.s-1"/>
        <field id="v"    long_name="Meridional wind"             unit="m.s-1"/>
        <field id="p"    long_name="Pressure"                    unit="Pa"/>
+       <field id="input_temp" long_name="Input temperature"  unit="K"/>
+       <field id="input_u" long_name="Input zonal wind"  unit="m/s"/>
+       <field id="input_v" long_name="Input meridional wind"  unit="m/s"/>
+       <field id="dtrad" long_name="Total radiative tendency" unit="K/s"/>
+       <field id="zdtlw" long_name="LW radiative tendency" unit="K/s"/>
+       <field id="zdtsw" long_name="SW radiative tendency" unit="K/s"/>
+       <field id="zdtdyn" long_name="Dyn temperature tendency" unit="K/s"/>
+       <field id="zdtdif" long_name="Turbulent diffusion tendency" unit="K/s"/>
+       <field id="zdtadj" long_name="Convective adjustement tendency" unit="K/s"/>
+       <field id="pdt" long_name="Total temperature tendency" unit="K/s"/>
      </field_group> 
      
@@ -23 +33 @@
     
     
-   <file_definition type="one_file" par_access="collective" output_freq="10ts" sync_freq="10ts" output_level="10" enabled=".TRUE.">
+   <file_definition type="one_file" par_access="collective" output_freq="1ts" sync_freq="1ts" output_level="10" enabled=".TRUE.">
 
      <file id="xios_diagfi" name="xios_diagfi" > 
@@ -34 +44 @@
     
     
+<!--   <axis_definition>
+     <axis id="altitude"/>
+   </axis_definition> -->
    <axis_definition>
-     <axis id="altitude"/>
+     <axis id="presnivs"/>
    </axis_definition>
     

codes/icosagcm/branches/SATURN_DYNAMICO/TEST/run_icosa.def

@@ -9 +9 @@
 
 # Number of subdivision on a main triangle (nbp) : integer (default=40)
-nbp=40
+nbp=10
 
 # nbp                 20  40  80 160
@@ -48 +48 @@
 itau_adv=3
 
-# number of timestep (default 100)
+# number of timestep (default 100) (ignored if run_length is set)
  itaumax = 100
 
@@ -55 +55 @@
 #run_length=0
 #run_length=432000
-# one saturn year
-run_length=929610360
+# one saturn year (24430 days)
+#run_length=929610360
+# 5 saturn day (5*38052=190260)
+#run_length=190260
+# 1 saturn day
+run_length=38052
+# 100 saturn days
+#run_length=3805200
+# 2000 saturn days
+#run_length=76104000
 
 #activate IO (default = true)
@@ -79 +87 @@
 # etat0 : initial state : string (default=jablonowsky06) : 
 # jablonowsky06, academic, ncar
-etat0=isothermal
-etat0_isothermal_temp=175
+#etat0=isothermal
+#etat0_isothermal_temp=175
+
+etat0=temperature_profile
+# for 'etat0=temperature_profile"
+# initial temperature profile provided in file 
+temperature_profile_file=temp_profile.txt
+
+#etat0=start_file
+# for 'etat0=start_file"
+# start file name (default is start.nc)
+# start_file_name=start
+# restart file name (default is restart.nc)
+#restart_file_name=restart
+
 # ------------------------------ Dynamics --------------------------------