Changeset 283 for trunk/LMDZ.MARS/libf/phymars/initracer.F

Timestamp:

Sep 7, 2011, 12:57:29 PM (13 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 ?

File:

• trunk/LMDZ.MARS/libf/phymars/initracer.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/initracer.F

@@ -43 +43 @@
 #include "advtrac.h"
 #include "comgeomfi.h"
-#include "watercap.h"
 #include "chimiedata.h"
 
@@ -50 +49 @@
       real qsurf(ngridmx,nqmx)       ! tracer on surface (e.g.  kg.m-2)
       real co2ice(ngridmx)           ! co2 ice mass on surface (e.g.  kg.m-2)
-      integer iq,ig,count, yeyey
+      integer iq,ig,count
       real r0_lift , reff_lift, nueff_lift
 c     Ratio of small over large dust particles (used when both 
@@ -436 +435 @@
         
 
-
       enddo ! of do iq=1,nqmx
 !      count=count+nbqchem
@@ -587 +585 @@
          alpha_lift(igcm_h2o_vap) =0.
          alpha_devil(igcm_h2o_vap)=0.
-
-c       "Dryness coefficient" controlling the evaporation and
-c        sublimation from the ground water ice (close to 1)
-c        HERE, the goal is to correct for the fact
-c        that the simulated permanent water ice polar caps
-c        is larger than the actual cap and the atmospheric
-c        opacity not always realistic.
-
-         do ig=1,ngridmx
-           if (ngridmx.ne.1) watercaptag(ig)=.false.
-           dryness(ig) = 1.
-         enddo
-
-         IF (caps) THEN
-c Perennial H20 north cap defined by watercaptag=true (allows surface to be
-c hollowed by sublimation in vdifc).
-         yeyey = 0
-         do ig=1,ngridmx
-!          !!! TESTS TESTS outliers
-!          !!! TESTS TESTS outliers
-!          if ( ( lati(ig)*180./pi      .ge.  75 ) .and.
-!     .         ( lati(ig)*180./pi      .le.  77 ) .and.
-!     .         ( ( ( long(ig)*180./pi .ge. 000. ) .and.
-!     .              ( long(ig)*180./pi .le. 120. ) ) 
-!     .             .or.
-!     .             ( ( long(ig)*180./pi .ge. -130. ) .and.
-!     .             ( long(ig)*180./pi .le. -115. ) ) ) ) then
-!             if (yeyey .eq. 0) then  !!! 1/2 en 64x48 sinon trop large en lat
-!              write(*,*) "outliers ", lati(ig)*180./pi, long(ig)*180./pi
-!              if (ngridmx.ne.1) watercaptag(ig)=.true.
-!              dryness(ig) = 1.
-!              albedodat(ig) = 0.45 !! comme alb_surfice
-!              yeyey = 1
-!             else
-!              yeyey = 0
-!             endif
-!          endif
-!          !!! TESTS TESTS outliers
-!          !!! TESTS TESTS outliers
-!
-!          !!! TESTS TESTS addcap
-!          !!! TESTS TESTS addcap     
-!          if ( ( lati(ig)*180./pi      .ge.  82 ) .and.
-!     .         ( lati(ig)*180./pi      .le.  84 ) .and.
-!     .         ( ( long(ig)*180./pi .gt. -030. ) .and.
-!     .              ( long(ig)*180./pi .lt. 090. ) ) ) then
-!              write(*,*) "capadd ", lati(ig)*180./pi, long(ig)*180./pi
-!              if (ngridmx.ne.1) watercaptag(ig)=.true.
-!              albedodat(ig) = 0.45 !! comme alb_surfice
-!              dryness(ig) = 1.
-!          endif
-!          !!! TESTS TESTS addcap
-!          !!! TESTS TESTS addcap 
-    
-           if (lati(ig)*180./pi.gt.84) then
-             if (ngridmx.ne.1) watercaptag(ig)=.true.
-             dryness(ig) = 1.
-c Use the following cap definition for high spatial resolution (latitudinal bin <= 5 deg)
-c             if (lati(ig)*180./pi.lt.85.and.long(ig).ge.0) then
-c               if (ngridmx.ne.1) watercaptag(ig)=.true.
-c               dryness(ig) = 1.
-c             endif
-c             if (lati(ig)*180./pi.ge.85) then
-c               if (ngridmx.ne.1) watercaptag(ig)=.true.
-c               dryness(ig) = 1.
-c             endif
-           endif  ! (lati>80 deg)
-         end do ! (ngridmx)
-        ENDIF ! (caps)
-
          if(water.and.(nqmx.ge.2)) then
            radius(igcm_h2o_ice)=3.e-6