Changeset 3147


Ignore:
Timestamp:
Dec 6, 2023, 9:00:08 AM (13 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:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/changelog.txt

    r3146 r3147  
    43744374(but used in testphys1d when atm_wat_profile > 0). It is fixed by calling
    43754375initracer in init_testphys1d_mod. 
     4376
     4377== 06/12/2023 == EM
     4378Minor fixes to ensure outputs of water cycle related variables are only
     4379computed and made if the water cycle is computed.
     4380 
  • trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

    r3144 r3147  
    26652665c        -------------------------------------------------------------------
    26662666
    2667        if (water) then
    26682667         do ig=1,ngrid
    26692668           if(mu0(ig).le.0.01) then
    26702669            fluxsurf_dir_dn_sw(ig) = 0.
    26712670           else
    2672             fluxsurf_dir_dn_sw(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*
     2671            if (water) then
     2672             ! both water and dust contribute
     2673             fluxsurf_dir_dn_sw(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*
    26732674     &                    exp(-(tau(ig,iaer_dust_doubleq)+
    26742675     &                          tau(ig,iaer_h2o_ice))/mu0(ig))
    2675            endif
     2676            else
     2677             ! only dust contributes
     2678             fluxsurf_dir_dn_sw(ig) = flux_1AU/dist_sol/dist_sol*mu0(ig)*
     2679     &                    exp(-(tau(ig,iaer_dust_doubleq))/mu0(ig))
     2680            endif ! of if (water)
     2681           endif ! of if(mu0(ig).le.0.01)
    26762682         enddo
    2677        endif
    26782683
    26792684           ! Density-scaled opacities
     
    32863291        call write_output("local_time","Local time",
    32873292     &                   'sol',local_time(:))
     3293        if (water) then
    32883294            if (.not.activice) then
    32893295               CALL write_output('tauTESap',
     
    32953301     &                         '',taucloudtes(:))
    32963302             endif
     3303        endif ! of if (water)
    32973304#else
    32983305      !!! this is to ensure correct initialisation of mesoscale model
Note: See TracChangeset for help on using the changeset viewer.