Changeset 1985 for trunk/LMDZ.MARS


Ignore:
Timestamp:
Aug 9, 2018, 4:27:52 PM (6 years ago)
Author:
mvals
Message:

Mars GCM:
Correction of the dust storm scheme:

  • in rocketduststorm_mod.F90 the case of not daytime and no storm was not taken into account (it was actually in comment), which lead to NaNs? in the calculation of pdqrds (scheme=0, division by alpha=0.)

Useful:

  • in compute_dtau_mod.F90, the writediagfi of tauref_scenario has been added

MV

Location:
trunk/LMDZ.MARS
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/README

    r1983 r1985  
    26322632Correction:
    26332633- in callsedim_mod, declaration of variable tau(ngrid,nlay) corrected to tau(ngrid,naerkind)
     2634
     2635== 09/08/2018 == MV
     2636Correction of the dust storm scheme:
     2637- in rocketduststorm_mod.F90 the case of not daytime and no storm was not taken into account (it was actually in comment), which lead to NaNs in the calculation of pdqrds (scheme=0, division by alpha=0.)
     2638Useful:
     2639- in compute_dtau_mod.F90, the writediagfi of tauref_scenario has been added

  • trunk/LMDZ.MARS/libf/phymars/compute_dtau_mod.F90

    r1974 r1985  
    5959        call read_dust_scenario(ngrid,nlayer,zday_scenario,pplev,       &
    6060       &                        tauref_scenario)
     61       ! for diagnostics
     62        call WRITEDIAGFI(ngrid,"tauref_scenario","tauref_scenario",     &
     63       &                  "",2,tauref_scenario)
    6164
    6265        ! 2. Compute dtau

  • trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90

    r1974 r1985  
    331331      DO ig=1,ngrid
    332332        !! **********************************************************************
    333         !! 2.1 case 1: Nothing to do when no storm and no slope
    334 !        IF ((mu0(ig) .LE. mu0lim) .AND. .NOT.(storm(ig)) ) then
    335 !          scheme(ig)=1
    336 !          cycle
    337 !        endif
    338         IF ((alpha_hmons(ig) .EQ. 0.) .AND. .NOT.(storm(ig))) then
     333        !! 2.1 case 1: Nothing to do when no storm and no slope, or
     334        !!             no storm and not daytime
     335        IF (((alpha_hmons(ig) .EQ. 0.) .AND. .NOT.(storm(ig))) &
     336             .OR. ((mu0(ig) .LE. mu0lim) .AND. .NOT.(storm(ig))) ) then
    339337          scheme(ig)=1
    340           cycle !!no slope
     338          cycle
    341339        endif
    342340       
Note: See TracChangeset for help on using the changeset viewer.