Changeset 283 for trunk/LMDZ.MARS/libf/phymars/meso_inc

Timestamp:

Sep 7, 2011, 12:57:29 PM (14 years ago)

Author:

aslmd

Message:

LMDZ.MARS:

AS: J'ai fait ce commit et fait a la main un merge avec les modifications entre temps

24/08/11 == TN

Attempts to tune the water cycle by adding outliers

+ A few structural changes !!

• watercap.h is now obsolete and removed -- all is in surfdat.h
• watercaptag initialized in surfini.F (up to 7 areas defined) instead of initracer.F
  • settings proposed by AS commented
  • experiments by TN decommented. use with caution.
• water ice albedo and thermal inertia in callphys.def and inifis.F
• water ice albedo in surfini.F
• water ice albedo computation in albedocaps.F90
• alb_surfice is now obsolete in physiq.F, albedo_h2o_ice is used instead
• frost_albedo_threshold defined in surfdat.h
• water ice thermal inertia in soil.F

TODO: * calibrate thermal inertia and ice albedo

• have a look at subgrid-scale ice with dryness ?

Location:

trunk/LMDZ.MARS/libf/phymars/meso_inc

Files:

• meso_inc_caps.F (modified) (2 diffs)
• meso_inc_ini.F (modified) (1 diff)
• meso_inc_var.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_caps.F

@@ +1 @@
+
+      !!!! see meso_inc_ini 
+      ! alb_lim = 0.26
+      ! lat_lim = 70.
+      ! inertie_lim = 800.
+
+      !!! This has to go after initracer which change dryness and watercaptag
+
       !!!!!! MARS MESOSCALE MODELING
       !!!!!! TEST TEST TEST TEST  AS+JBM 28/02/11
@@ -14 +22 @@
       !!!!           du GCM lorsqu'ils sont consequents
       !!!! 
-      IF ( caps .and. (igcm_h2o_ice .ne. 0) ) THEN
+      !IF ( caps .and. (igcm_h2o_ice .ne. 0) ) THEN
+      IF ( caps .and. water ) THEN
           PRINT *, 'OVERWRITING watercaptag DEFINITION in INITRACER'
-          PRINT *, 'lat>70 et alb>0.26 => watercaptag=T' 
+          PRINT *, 'lat > lat_lim et alb > alb_lim => watercaptag=T'
+          PRINT *, 'ind for water ice: ', igcm_h2o_ice 
           !! Perennial H20 north cap defined by watercaptag=true (allows surface to be
           !! hollowed by sublimation in vdifc).
           do ig=1,ngridmx
             qsurf(ig,igcm_h2o_ice)=0.  !! on jette les inputs GCM
-            if ( (lati(ig)*180./pi.gt.70.) .and.
-     .           (albedodat(ig).ge.0.26) )  then
-                    watercaptag(ig)=.true.
-                    dryness(ig) = 1.
+            if ( ( lati(ig)*180./pi .gt. lat_lim ) .and.
+     .           ( albedodat(ig) .ge. alb_lim    ) )  then
+                    watercaptag(ig)  = .true.
+                    dryness(ig)      = 1.
             else
-                    watercaptag(ig)=.false.
-                    dryness(ig) = 1.
+                    watercaptag(ig)  = .false.
+                    dryness(ig)      = 1.
             endif  ! (lati, albedodat)
           end do ! (ngridmx)

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_ini.F

@@ -68 +68 @@
 c
 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      !!! see meso_inc_caps
+      !!! this is a test to change outliers' albedo and thermal inertia
+      alb_lim = 0.26
+      lat_lim = 70.
+      inertie_lim = 800. 
+      PRINT *, 'lat_lim ',lat_lim
+      PRINT *, 'alb_lim ',alb_lim
+      PRINT *, 'inertie_lim ',inertie_lim 
+      !!!
+      !!!
+      IF ( caps .and. water ) THEN
+          do ig=1,ngridmx
+            if ( lati(ig)*180./pi .gt. lat_lim ) then
+               if ( albedodat(ig) .ge. alb_lim ) then
+                    albedodat(ig) = alb_surfice
+                    inertiedat(ig,1) = inertie_lim  
+               endif
+               if (inertiedat(ig,1) .ge. inertie_lim ) then
+                    inertiedat(ig,1) = inertie_lim
+               endif
+            endif  ! (lati, albedodat)
+          end do ! (ngridmx)
+      ENDIF ! (caps)

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_var.F

@@ -18 +18 @@
       INTEGER tracerset    !!! this corresponds to config%mars
       CHARACTER (len=20) :: wtnom(nqmx) ! tracer name
+
+      REAL alb_lim 
+      REAL lat_lim 
+      REAL inertie_lim