Changeset 3147 for trunk/LMDZ.MARS

Timestamp:

Dec 6, 2023, 9:00:08 AM (12 months ago)

Author:

emillour

Message:

Mars PCM:
Minor fixes to ensure outputs of water cycle related variables are only
computed and made if the water cycle is computed.
EM

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/physiq_mod.F (modified) (3 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4374 +4374 @@
 (but used in testphys1d when atm_wat_profile > 0). It is fixed by calling
 initracer in init_testphys1d_mod.  
+
+== 06/12/2023 == EM
+Minor fixes to ensure outputs of water cycle related variables are only
+computed and made if the water cycle is computed.
+ 

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

@@ -2665 +2665 @@
 c        -------------------------------------------------------------------
 
-       if (water) then
          do ig=1,ngrid
            if(mu0(ig).le.0.01) then
             fluxsurf_dir_dn_sw(ig) = 0.
            else
-            fluxsurf_dir_dn_sw(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*
+            if (water) then
+             ! both water and dust contribute
+             fluxsurf_dir_dn_sw(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*
      &                    exp(-(tau(ig,iaer_dust_doubleq)+
      &                          tau(ig,iaer_h2o_ice))/mu0(ig))
-           endif
+            else
+             ! only dust contributes
+             fluxsurf_dir_dn_sw(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*
+     &                    exp(-(tau(ig,iaer_dust_doubleq))/mu0(ig))
+            endif ! of if (water)
+           endif ! of if(mu0(ig).le.0.01)
          enddo
-       endif
 
            ! Density-scaled opacities
@@ -3286 +3291 @@
         call write_output("local_time","Local time",
      &                   'sol',local_time(:))
+        if (water) then
             if (.not.activice) then
                CALL write_output('tauTESap',
@@ -3295 +3301 @@
      &                         '',taucloudtes(:))
              endif
+        endif ! of if (water)
 #else
       !!! this is to ensure correct initialisation of mesoscale model