Ignore:
Timestamp:
Jan 22, 2019, 5:02:23 PM (6 years ago)
Author:
mvals
Message:

Mars GCM:

  • Update of rocketduststorm_mod.F90 :

We want to separate both parametrizations related to the formation of the detached dust layers. Therefore, rocketduststorm_mod.F90 now only comprises the rocket dust storm scheme, whereas it contained
also before the calculation of the vertical velocity induced by the presence of the sub-grid scale topography. This latter part is under development and will be integrated as a fully independant
parametrization: the aim is to simulate the entrainment of dust by slope winds, from the boundary layer up to the top of the sub-grid scale topography.

  • Addition of initial surface parameters "summit" and "base" to prepare the previously described slope wind parametrization

MV

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

    r1974 r2079  
    1111     &                     albedice, iceradius, dtemisice, z0,
    1212     &                     zmea, zstd, zsig, zgam, zthe, phisfi,
    13      &                     watercaptag, hmons
     13     &                     watercaptag, hmons, summit, base
    1414      use slope_mod, only: theta_sl, psi_sl
    1515      use phyredem, only: physdem0,physdem1
     
    556556      zsig(1)=0.E+0
    557557      zgam(1)=0.E+0
    558       zthe(1)=0.E+0
    559 c  for rocket dust storm scheme     
     558      zthe(1)=0.E+0     
    560559      hmons(1)=0.E+0
    561 
    562 
     560      summit(1)=0.E+0
     561      base(1)=0.E+0
     562     
    563563c   Specific initializations for "physiq"
    564564c   -------------------------------------
     
    718718      call physdem0("startfi.nc",longitude,latitude,nsoilmx,ngrid,llm,
    719719     &              nq,dtphys,float(day0),time,cell_area,
    720      &              albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe,hmons)
     720     &              albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe,
     721     &              hmons,summit,base)
    721722      call physdem1("startfi.nc",nsoilmx,ngrid,llm,nq,
    722723     &              dtphys,time,
Note: See TracChangeset for help on using the changeset viewer.