Changeset 80 for trunk/mars/libf

Timestamp:

Feb 28, 2011, 2:54:20 PM (14 years ago)

Author:

aslmd

Message:

USERS: ajout 2 personnes

LMD_MM_MARS: changements mineurs

• fichiers de reglages dans 'newphys_tracers' pour test THARSIS_new_water_cycle [concluant]
• ajout de callphys.def.mars11.radact
• ajouts des launchers sur GNOME pour sequentiel et parallele
• changement script runmeso suite a revision 79 (correction bug datadir)
• tests d'une option 'caps' mesoscale dans meso_physiq

File:

• trunk/mars/libf/phymars/meso_physiq.F (modified) (2 diffs)

trunk/mars/libf/phymars/meso_physiq.F

@@ -432 +432 @@
       !  PRINT*, tnom(iq), pq(:,:,iq)
       !ENDDO
+
 c
 c ****WRF
@@ -482 +483 @@
             CALL initracer(qsurf,co2ice)
          ENDIF  ! end tracer
+
+      !!!!!! WRF WRF WRF MARS MARS 
+      !!!!!! 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
+      !!!! 
+      !!!! principe: une option 'caps=T' specifique au mesoscale
+      !!!! ... en vue d'un meso_initracer ????
+      !!!! 
+      !!!! depots permanents => albedo TES du PDS
+      !!!! depots saisonniers => alb_surfice (~0.4, cf plus bas)
+      !!!!     [!!!! y compris pour les depots saisonniers sur les depots permanents]
+      !!!!
+      !!!! --> todo: il faut garder les depots saisonniers qui viennent
+      !!!!           du GCM lorsqu'ils sont consequents
+      !!!! 
+      IF ( caps .and. (igcm_h2o_ice .ne. 0) ) THEN
+          PRINT *, 'OVERWRITING watercaptag DEFINITION in INITRACER'
+          PRINT *, 'lat>70 et alb>0.26 => watercaptag=T' 
+          !! 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.
+            else
+                    watercaptag(ig)=.false.
+                    dryness(ig) = 1.
+            endif  ! (lati, albedodat)
+          end do ! (ngridmx)
+      ELSE  ! (caps)
+          print *,'Blork !!!'
+          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