Changeset 1202 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd

Timestamp:

Mar 7, 2014, 11:51:22 AM (11 years ago)

Author:

aslmd

Message:

MESOSCALE. old physics for polar runs: better handling of Titus cap and handling of Tyler surface properties

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars

Files:

• inifis.F (modified) (1 diff)
• physiq.F (modified) (4 diffs)
• surfini.F (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/inifis.F

@@ -504 +504 @@
          endif
 
+         READ(99,*)
+         READ(99,*)
+
+         !! prescribed crocus line for caps
+         READ(99,fmt='(a)') ch1
+         READ(99,*,iostat=ierr) tituscap
+         if (ierr.ne.0) then
+           write(*,*) 'no titus cap requested. set it to false'
+           tituscap = .false.
+         endif 
+         WRITE(*,8000) ch1,tituscap
+
          PRINT*,'--------------------------------------------'
          PRINT*

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/physiq.F

@@ -330 +330 @@
       integer igmin, lmin
       logical tdiag
-      logical tituscap
 
       REAL zplev(ngrid,nlayermx+1),zplay(ngrid,nlayermx)
@@ -345 +344 @@
      
 c=======================================================================
+
+
 
 
@@ -1064 +1065 @@
 c   -------------------------------------------
 
-      tituscap = .true.
+      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+      !!! get the actual co2 seasonal cap from Titus observations
+      !!! (inherited from GCM only at first step, but not a big deal)
       IF (tituscap) THEN
-         PRINT*,'TITUS CAP !!!'
-         !!! get the actual co2 seasonal cap from Titus observations
          CALL geticecover( ngrid, 180.*zls/pi,
      .                  180.*long/pi, 180.*lati/pi, co2ice )
-         co2ice = co2ice * 10000.
+         co2ice(:) = co2ice(:)*10000.
+         emis(:) = 0.95 ! so that points outside the cap are indeed at 0.95
+                        ! avoid unwanted patchiness from GCM initial state
       ENDIF
+      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
       IF(callcond) THEN
@@ -1081 +1085 @@
      $              fluxsurf_sw) 
 
+         !! Titus cap security. 
+         !! -- neglect pressure changes caused by sublimation/condensation
          IF (tituscap) THEN 
            DO ig=1,ngrid

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/surfini.F

@@ -37 +37 @@
                           !! si elles sont consequentes
          enddo
-         if ( ( lati(ig)*180./pi .gt. 70. ) .and.
-     .        ( albedodat(ig) .ge. 0.26   ) )  then
-                 write(*,*)"outlier ",ig,albedodat(ig),inertiedat(ig)
+
+!! AS: my previous stuff
+!         if ( ( lati(ig)*180./pi .gt. 70. ) .and.
+!     .        ( albedodat(ig) .ge. 0.26   ) )  then
+!                 write(*,*)"outlier ",ig,albedodat(ig),inertiedat(ig)
+!                 watercaptag(ig)  = .true.
+!                 dryness(ig)      = 1.
+!                 albedodat(ig)    = 0.45 !albedo_h2o_ice  !! pour output 
+!                 inertiedat(ig)   = 800.
+!                 write(*,*)"new values ",ig,albedodat(ig),inertiedat(ig)
+!         else
+!                 watercaptag(ig)  = .false.
+!                 dryness(ig)      = 1.
+!         endif 
+!         if ( inertiedat(ig) .ge. 800. ) then
+!            write(*,*)"change inertie from ",inertiedat(ig)," to 800."
+!            inertiedat(ig)   = 800.
+!         endif
+
+!! using Tyler and Barnes maps.
+         dryness(ig)      = 1.
+         if ( inertiedat(ig) .ge. 600. ) then
+                 write(*,*)"ice ",ig,albedodat(ig),inertiedat(ig)
                  watercaptag(ig)  = .true.
-                 dryness(ig)      = 1.
-                 albedodat(ig)    = 0.45 !albedo_h2o_ice  !! pour output 
-                 inertiedat(ig)   = 800.
-                 write(*,*)"new values ",ig,albedodat(ig),inertiedat(ig)
          else
                  watercaptag(ig)  = .false.
-                 dryness(ig)      = 1.
-         endif 
-         if ( inertiedat(ig) .ge. 800. ) then
-            write(*,*)"change inertie from ",inertiedat(ig)," to 800."
-            inertiedat(ig)   = 800.
          endif