Changeset 1973 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Jul 17, 2018, 4:38:20 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:

• follow-up of previous commit.

MV

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (2 diffs)

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

@@ -408 +408 @@
       real fluxtop_swclf(ngrid,2) 
       real taucloudtesclf(ngrid)
-      real tf, ntf ! tf: fraction of clouds, ntf: fraction without clouds
+      real tf_clf, ntf_clf ! tf: fraction of clouds, ntf: fraction without clouds
       real rave2(ngrid), totrave2(ngrid) ! Mean water ice mean radius (m)
 
@@ -751 +751 @@
                ! cases
                DO ig=1,ngrid
-                  tf=totcloudfrac(ig)
-                  ntf=1.-tf
-                  fluxsurf_lw(ig) = ntf*fluxsurf_lwclf(ig) 
-     &                                          + tf*fluxsurf_lw(ig)
-                  fluxsurf_sw(ig,1) = ntf*fluxsurf_swclf(ig,1) 
-     &                                          + tf*fluxsurf_sw(ig,1)
-                  fluxsurf_sw(ig,2) = ntf*fluxsurf_swclf(ig,2) 
-     &                                          + tf*fluxsurf_sw(ig,2)
-                  fluxtop_lw(ig)  = ntf*fluxtop_lwclf(ig) 
-     &                                          + tf*fluxtop_lw(ig)
-                  fluxtop_sw(ig,1)  = ntf*fluxtop_swclf(ig,1)  
-     &                                          + tf*fluxtop_sw(ig,1)
-                  fluxtop_sw(ig,2)  = ntf*fluxtop_swclf(ig,2)  
-     &                                          + tf*fluxtop_sw(ig,2)
-                  taucloudtes(ig) = ntf*taucloudtesclf(ig) 
-     &                                          + tf*taucloudtes(ig)
-                  zdtlw(ig,1:nlayer) = ntf*zdtlwclf(ig,1:nlayer) 
-     &                                         + tf*zdtlw(ig,1:nlayer)
-                  zdtsw(ig,1:nlayer) = ntf*zdtswclf(ig,1:nlayer) 
-     &                                         + tf*zdtsw(ig,1:nlayer)
+                  tf_clf=totcloudfrac(ig)
+                  ntf_clf=1.-tf_clf
+                  fluxsurf_lw(ig) = ntf_clf*fluxsurf_lwclf(ig) 
+     &                                      + tf_clf*fluxsurf_lw(ig)
+                  fluxsurf_sw(ig,1) = ntf_clf*fluxsurf_swclf(ig,1) 
+     &                                      + tf_clf*fluxsurf_sw(ig,1)
+                  fluxsurf_sw(ig,2) = ntf_clf*fluxsurf_swclf(ig,2) 
+     &                                      + tf_clf*fluxsurf_sw(ig,2)
+                  fluxtop_lw(ig)  = ntf_clf*fluxtop_lwclf(ig) 
+     &                                      + tf_clf*fluxtop_lw(ig)
+                  fluxtop_sw(ig,1)  = ntf_clf*fluxtop_swclf(ig,1)  
+     &                                      + tf_clf*fluxtop_sw(ig,1)
+                  fluxtop_sw(ig,2)  = ntf_clf*fluxtop_swclf(ig,2)  
+     &                                      + tf_clf*fluxtop_sw(ig,2)
+                  taucloudtes(ig) = ntf_clf*taucloudtesclf(ig) 
+     &                                      + tf_clf*taucloudtes(ig)
+                  zdtlw(ig,1:nlayer) = ntf_clf*zdtlwclf(ig,1:nlayer) 
+     &                                      + tf_clf*zdtlw(ig,1:nlayer)
+                  zdtsw(ig,1:nlayer) = ntf_clf*zdtswclf(ig,1:nlayer) 
+     &                                      + tf_clf*zdtsw(ig,1:nlayer)
                ENDDO