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

Timestamp:

Sep 8, 2011, 1:59:47 AM (13 years ago)

Author:

aslmd

Message:

08/09/11 == AS

LMDZ.MARS + MESOSCALE.
---> Setting up a more realistic water ice source at the poles (notably outliers)

surfini.F ?
Main changes and bug fixes

• reference to comcstfi.h was wrong. big problem because e.g. pi was not known.
• commented about a problem to be fixed, due to surfini being called before initracer.
• MESOSCALE: put the mesoscale north cap definition into a precompiling flag #MESOSCALE

for the moment: if [alb_mean_TES > 0.26 and lat > 70] then outliers
(previously done in meso_inc_caps.F)

inifis.F ?
Just changed a comment with wrong formatting

--> below, only MESOSCALE

soil.F ?
if somewhere IT > IT_outliers, then makes it = IT_outliers

physiq.F ?
meso_inc/meso_inc_caps.F ?
meso_inc/meso_inc_ini.F ?
meso_inc_caps no longer called. keep for reference for the moment.

meso_inc/meso_inc_var.F ?
deleted lines with *_lim variables, now useless

Location:

trunk/LMDZ.MARS/libf/phymars/meso_inc

Files:

• meso_inc_caps.F (modified) (4 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
 
@@ -23 +17 @@
       !!!! 
       !IF ( caps .and. (igcm_h2o_ice .ne. 0) ) THEN
+
       IF ( caps .and. water ) THEN
           PRINT *, 'OVERWRITING watercaptag DEFINITION in INITRACER'
@@ -31 +26 @@
           do ig=1,ngridmx
             qsurf(ig,igcm_h2o_ice)=0.  !! on jette les inputs GCM
-            if ( ( lati(ig)*180./pi .gt. lat_lim ) .and.
-     .           ( albedodat(ig) .ge. alb_lim    ) )  then
+            if ( ( lati(ig)*180./pi .gt. 70. ) .and.
+     .           ( albedodat(ig) .ge. 0.26   ) )  then
                     watercaptag(ig)  = .true.
                     dryness(ig)      = 1.
@@ -44 +39 @@
           print *,'caps=T avec water=F ????'
       ENDIF ! (caps)
-      !!!!!! TEST TEST TEST TEST  AS+JBM 28/02/11
-      !!!!!! TEST TEST TEST TEST  AS+JBM 28/02/11
-      !!!!!! TEST TEST TEST TEST  AS+JBM 28/02/11

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

@@ -69 +69 @@
 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)
+      !!!! 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