Changeset 1208 for trunk/LMDZ.MARS/libf/phymars/updatereffrad.F

Timestamp:

Mar 12, 2014, 4:30:18 PM (11 years ago)

Author:

tnavarro

Message:

added tauscaling in startfi + moved start_archive routines from dyn3d to phymars

File:

• trunk/LMDZ.MARS/libf/phymars/updatereffrad.F (modified) (1 diff)

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

@@ -128 +128 @@
         IF (water.AND.activice) THEN 
          IF (microphys) THEN
-         
-          IF (firstcall) THEN
-            DO l=1,nlayer
-              DO ig=1,ngrid
-                call updaterice_micro(pq(ig,l,igcm_h2o_ice),
-     &                                pq(ig,l,igcm_ccn_mass),
-     &                                pq(ig,l,igcm_ccn_number),
-     &                                1.e-3,rice(ig,l),
-     &                                rhocloud(ig,l))
-                nuice(ig,l) = nuice_ref
-              ENDDO
-            ENDDO
-          firstcall = .false.
+        
 c    At firstcall, the true number and true mass of cloud condensation nuclei are not known.
 c    Indeed it is scaled on the prescribed dust opacity via a 'tauscaling' coefficient
-c    computed after radiative transfer.
-          ELSE
-            DO l=1,nlayer
-              DO ig=1,ngrid
-                call updaterice_micro(pq(ig,l,igcm_h2o_ice),
-     &                                pq(ig,l,igcm_ccn_mass),
-     &                                pq(ig,l,igcm_ccn_number),
-     &                                tauscaling(ig),rice(ig,l),
-     &                                rhocloud(ig,l))
-                nuice(ig,l) = nuice_ref
-               ENDDO
-             ENDDO
-          ENDIF ! of if firstcall
+c    computed after radiative transfer. If tauscaling is not in startfi, we make an assumption for its value.
+
+          IF (firstcall) THEN
+            IF (minval(tauscaling).lt.0) tauscaling(:) = 1.e-3 ! default value when non-read in startfi is -1
+            IF (freedust)                tauscaling(:) = 1.    ! if freedust, enforce no rescaling at all
+            firstcall = .false.
+          ENDIF
+ 
+          DO l=1,nlayer
+            DO ig=1,ngrid
+              call updaterice_micro(pq(ig,l,igcm_h2o_ice),
+     &                              pq(ig,l,igcm_ccn_mass),
+     &                              pq(ig,l,igcm_ccn_number),
+     &                              tauscaling(ig),rice(ig,l),
+     &                              rhocloud(ig,l))
+              nuice(ig,l) = nuice_ref
+            ENDDO
+          ENDDO
           
         ELSE ! if not microphys