Changeset 283 for trunk/LMDZ.MARS/libf/phymars/physiq.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/physiq.F (modified) (7 diffs)

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

@@ -128 +128 @@
 
 #include "chimiedata.h"
-#include "watercap.h"
 #include "param.h"
 #include "param_v3.h"
@@ -191 +190 @@
       REAL fluxgrd(ngridmx)          ! surface conduction flux (W.m-2)
       REAL qsurf(ngridmx,nqmx)       ! tracer on surface (e.g. kg.m-2)
-      REAL q2(ngridmx,nlayermx+1)    ! Turbulent Kinetic Energy 
+      REAL q2(ngridmx,nlayermx+1)    ! Turbulent Kinetic Energy
+      
+      REAL watercapflag(ngridmx)     ! water cap flag
 
 c     Variables used by the water ice microphysical scheme:
@@ -197 +198 @@
       REAL nuice(ngridmx,nlayermx)   ! Estimated effective variance
                                      !   of the size distribution
-c     Albedo of deposited surface ice
-      !!REAL, PARAMETER :: alb_surfice = 0.4 ! 0.45
-      REAL, PARAMETER :: alb_surfice = 0.45 !!TESTS_JB 
 
 c     Variables used by the slope model
@@ -422 +420 @@
 
         IF (tracer.AND.water.AND.(ngridmx.NE.1)) THEN
-          write(*,*)"physiq: water_param Surface ice alb:",alb_surfice
+          write(*,*)"physiq: water_param Surface water ice albedo:", 
+     .                  albedo_h2o_ice
         ENDIF
                    
@@ -1095 +1094 @@
 #endif
 c       -------------------------------------------------------------
-c       Change of surface albedo (set to 0.4) in case of ground frost
+c       Change of surface albedo in case of ground frost
 c       everywhere except on the north permanent cap and in regions
 c       covered by dry ice. 
@@ -1102 +1101 @@
          do ig=1,ngrid
            if ((co2ice(ig).eq.0).and.
-     &        (qsurf(ig,igcm_h2o_ice).gt.0.005)) then
-              albedo(ig,1) = alb_surfice
-              albedo(ig,2) = alb_surfice
+     &        (qsurf(ig,igcm_h2o_ice).gt.frost_albedo_threshold)) then
+              albedo(ig,1) = albedo_h2o_ice
+              albedo(ig,2) = albedo_h2o_ice
+c              write(*,*) "frost thickness", qsurf(ig,igcm_h2o_ice)
+c              write(*,*) "physiq.F frost :"
+c     &        ,lati(ig)*180./pi, long(ig)*180./pi
            endif
          enddo  ! of do ig=1,ngrid
@@ -1512 +1514 @@
      &                       'surface h2o_ice',
      &                       'kg.m-2',2,qsurf(1,igcm_h2o_ice))
+
+            if (caps) then
+             do ig=1,ngridmx
+                if (watercaptag(ig)) watercapflag(ig) = 1
+             enddo
+             CALL WRITEDIAGFI(ngridmx,'watercaptag',
+     &                         'Ice water caps',
+     &                         '',2,watercapflag)
+            endif
+            CALL WRITEDIAGFI(ngridmx,'albedo',
+     &                         'albedo',
+     &                         '',2,albedo(1:ngridmx,1))
            endif !(water)