Changeset 2915 for LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

Timestamp:

Jun 16, 2017, 6:03:16 PM (7 years ago)

Author:

jbmadeleine

Message:

Improved the way thermal inertia works in simulations using imposed SST and sea ice.
Removed inertie_ice which was confusing; Added inertie_sic for sea-ice and inertie_lic for
land-ice. When snow is present, inertie_sno is used instead of inertie_sol, inertie_sic or
inertie_lic. For sea-ice, added a flag called iflag_sic:
If iflag_sic=0, thermal inertia is changed over sea-ice when snow is present
If iflag_sic=1, thermal inertia is kept constant and equal to inertie_sic

File:

• LMDZ5/trunk/libf/phylmd/conf_phys_m.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -173 +173 @@
     REAL,SAVE :: exposant_glace_omp
     REAL,SAVE :: rei_min_omp, rei_max_omp
-    REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_ice_omp
+    INTEGER,SAVE :: iflag_sic_omp
+    REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_sic_omp
+    REAL,SAVE :: inertie_lic_omp
     REAL,SAVE :: qsol0_omp
     REAL,SAVE :: evap0_omp
@@ -1133 +1135 @@
     !-----------------------------------------------------------------------
     !
-    !Config Key  = inertie_ice
+    !Config Key  = iflag_sic
+    !Config Desc =  
+    !Config Def  = 0
+    !Config Help = 
+    !
+    iflag_sic_omp = 0
+    CALL getin('iflag_sic',iflag_sic_omp)
+    !
+    !Config Key  = inertie_sic
     !Config Desc =  
     !Config Def  = 2000.
     !Config Help = 
     !
-    inertie_ice_omp = 2000.
-    CALL getin('inertie_ice',inertie_ice_omp)
+    inertie_sic_omp = 2000.
+    CALL getin('inertie_sic',inertie_sic_omp)
+    !
+    !Config Key  = inertie_lic
+    !Config Desc =  
+    !Config Def  = 2000.
+    !Config Help = 
+    !
+    inertie_lic_omp = 2000.
+    CALL getin('inertie_lic',inertie_lic_omp)
     !
     !Config Key  = inertie_sno
@@ -2132 +2150 @@
     evap0 = evap0_omp
     albsno0 = albsno0_omp
+    iflag_sic = iflag_sic_omp
     inertie_sol = inertie_sol_omp
-    inertie_ice = inertie_ice_omp
+    inertie_sic = inertie_sic_omp
+    inertie_lic = inertie_lic_omp
     inertie_sno = inertie_sno_omp
     rad_froid = rad_froid_omp
@@ -2545 +2565 @@
     write(lunout,*)' evap0 = ', evap0
     write(lunout,*)' albsno0 = ', albsno0
+    write(lunout,*)' iflag_sic = ', iflag_sic
     write(lunout,*)' inertie_sol = ', inertie_sol
-    write(lunout,*)' inertie_ice = ', inertie_ice
+    write(lunout,*)' inertie_sic = ', inertie_sic
+    write(lunout,*)' inertie_lic = ', inertie_lic
     write(lunout,*)' inertie_sno = ', inertie_sno
     write(lunout,*)' f_cdrag_ter = ',f_cdrag_ter