Changeset 83 for trunk/mars

Timestamp:

Mar 1, 2011, 4:53:50 PM (14 years ago)

Author:

aslmd

Message:

mars + LMD_MM_MARS : modifications pour cycle de l'eau, valeurs tunees JBM

---

used settings reached by JBM to obtain his PhD results

alb_surfice = 0.45 --- in physiq.F and meso_physiq.F
ccn_factor = 4.5 --- in watercloud.F
nuice_sed = 0.45 --- in callsedim.F

NB: this is supposed to be further refined in the future

Location:

trunk/mars

Files:

• README (modified) (2 diffs)
• libf/phymars/aeropacity.F (modified) (2 diffs)
• libf/phymars/meso_physiq.F (modified) (1 diff)
• libf/phymars/physiq.F (modified) (1 diff)
• libf/phymars/watercloud.F (modified) (3 diffs)

trunk/mars/README

@@ -588 +588 @@
    to external datafiles (topography, surface properties, etc.)
 
-== 28/02/2011 == JBM+AS
+== 28/02/2011 == JBM + AS
 >> used settings reached by JBM to obtain his PhD results
       alb_surfice = 0.45  ---  in physiq.F and meso_physiq.F
@@ -594 +594 @@
       nuice_sed = 0.45    ---  in callsedim.F
 >> NB: this is supposed to be further refined in the future
+
+== 01/03/2011 == AS + JBM
+>> nasty bug in the water cycle when iradia != 1  [no problem when iradia = 1]
+   --> mesoscale runs w/ water cycle had strange 5-hour fluctuations in RICE, from 80mic to 5mic
+>> PB: calculation of ccn [condensation nuclei] is done in callradite.F
+      * ccn must be saved 
+        --> corrected in physiq.F and meso_physiq.F
+      * ccn must not be modified elsewhere [e.g. in watercloud, when divided by ccn_factor] 
+        --> all calculations on ccn are now moved in callradite

trunk/mars/libf/phymars/aeropacity.F

@@ -106 +106 @@
                                    !   (particules kg-1)
       REAL qtot(ngridmx)           ! Dust column (kg m-2)
+
+c     CCN reduction factor
+      REAL, PARAMETER :: ccn_factor = 4.5  !! comme TESTS_JB // 1. avant
+
 c
 c   local saved variables
@@ -476 +480 @@
       ENDDO ! iaer (loop on aerosol kind)
 
+
+c -----------------------------------------------------------------
+c -----------------------------------------------------------------
+c  Reduce number of nuclei
+!         TEMPORAIRE : r�duction du nombre de nuclei FF 04/200
+!         reduction facteur 3
+!         ccn(ig,l) = ccn(ig,l) / 27.
+!         reduction facteur 2
+!         ccn(ig,l) = ccn(ig,l) / 8.
+c -----------------------------------------------------------------
+       write(*,*) "water_param CCN reduc. fac. ", ccn_factor
+       DO l=1,nlayer
+         DO ig=1,ngrid
+            ccn(ig,l) = ccn(ig,l) / ccn_factor
+         ENDDO
+       ENDDO
+c -----------------------------------------------------------------
+c -----------------------------------------------------------------
+
+
 c -----------------------------------------------------------------
 c Column integrated visible optical depth in each point

trunk/mars/libf/phymars/meso_physiq.F

@@ -330 +330 @@
       REAL ccn(ngridmx,nlayermx)   ! Cloud condensation nuclei
                                    !   (particules kg-1)
+      SAVE ccn  !! in case iradia != 1 
       real rdust(ngridmx,nlayermx) ! dust geometric mean radius (m)
       real qtot1,qtot2 ! total aerosol mass

trunk/mars/libf/phymars/physiq.F

@@ -279 +279 @@
       REAL ccn(ngridmx,nlayermx)   ! Cloud condensation nuclei
                                    !   (particules kg-1)
+      SAVE ccn  !! in case iradia != 1
       real rdust(ngridmx,nlayermx) ! dust geometric mean radius (m)
       real qtot1,qtot2 ! total aerosol mass

trunk/mars/libf/phymars/watercloud.F

@@ -55 +55 @@
       REAL ccn(ngridmx,nlayermx)   ! Cloud condensation nuclei
                                    !   (particules kg-1)
-c     CCN reduction factor
-      REAL, PARAMETER :: ccn_factor = 4.5  !! comme TESTS_JB // 1. avant
       real rdust(ngridmx,nlayermx) ! Dust geometric mean radius (m)
 
@@ -137 +135 @@
         write(*,*) "            i_ice=",i_ice
 
-        write(*,*) "water_param CCN reduction factor:", ccn_factor
- 
         firstcall=.false.
       ENDIF ! of IF (firstcall)
@@ -172 +168 @@
       do l=1,nlay
         do ig=1,ngrid
-          ccn(ig,l) = ccn(ig,l) / ccn_factor
-c         TEMPORAIRE : réduction du nombre de nuclei FF 04/2008 :
-c         reduction facteur 3
-c         ccn(ig,l) = ccn(ig,l) / 27.
-c         reduction facteur 2
-c         ccn(ig,l) = ccn(ig,l) / 8.
 c         Calcul du rayon moyen des particules de glace.
 c         Hypothese : Dans une couche, la glace presente se