Changeset 1974 for trunk/LMDZ.MARS/libf/phymars/callradite_mod.F

Timestamp:

Jul 18, 2018, 4:48:34 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:
Integration of the detached dust layer parametrizations (rocket dust storm, slope wind lifting, CW, and dust injection scheme, DB).
Still experimental, default behaviour (rdstorm=.false., dustinjection=0) identical to previous revision.
NB: Updated newstart requires an updated "surface.nc" containing the "hmons" field.
EM+MV

File:

• trunk/LMDZ.MARS/libf/phymars/callradite_mod.F (modified) (10 diffs)

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

@@ -7 +7 @@
       SUBROUTINE callradite(icount,ngrid,nlayer,nq,zday,ls,pq,albedo,
      $     emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout,
-     $     dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,fluxtop_sw,
-     &     tauref,tau,aerosol,dsodust,tauscaling,taucloudtes,rdust,rice,
-     &     nuice,co2ice,clearsky,totcloudfrac)
+     $     dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,
+     $     fluxtop_sw,tauref,tau,aerosol,dsodust,tauscaling,
+     $     taucloudtes,rdust,rice,nuice,co2ice,rstormdust,
+     $     totstormfract,clearatm,dsords,
+     $     clearsky,totcloudfrac)
 
       use aeropacity_mod, only: aeropacity
@@ -16 +18 @@
       use dimradmars_mod, only: naerkind, name_iaer,
      &            iaer_dust_conrath,iaer_dust_doubleq,
-     &            iaer_dust_submicron,iaer_h2o_ice
+     &            iaer_dust_submicron,iaer_h2o_ice,
+     &            iaer_stormdust_doubleq 
       use yomlw_h, only: gcp, nlaylte
       use comcstfi_h, only: g,cpp
@@ -154 +157 @@
 c                         the "naerkind" kind of aerosol optical
 c                         properties.
-
 c=======================================================================
 c
@@ -160 +162 @@
 c    -------------
 c
-#include "callkeys.h"
+      include "callkeys.h"
 
 c-----------------------------------------------------------------------
@@ -170 +172 @@
 
       REAL,INTENT(IN) :: pq(ngrid,nlayer,nq)
-      REAL,INTENT(IN) :: tauscaling(ngrid) ! Conversion factor for 
+      REAL,INTENT(INOUT) :: tauscaling(ngrid) ! Conversion factor for 
                                ! qdust and Ndust
       REAL,INTENT(IN) :: albedo(ngrid,2),emis(ngrid)
@@ -194 +196 @@
       REAL,INTENT(OUT) :: nuice(ngrid,nlayer)  ! Estimated effective variance
       REAL,INTENT(IN) :: co2ice(ngrid)           ! co2 ice surface layer (kg.m-2)
+
+c     rocket dust storm
+      LOGICAL,INTENT(IN) :: clearatm ! true for background dust
+      REAL,INTENT(IN) :: totstormfract(ngrid) ! dust storm mesh fraction
+      REAL,INTENT(OUT) :: rstormdust(ngrid,nlayer)  ! Storm dust geometric mean radius (m)
+      REAL dsords(ngrid,nlayer) ! density scaled opacity for rocket dust storm dust
+
 c     sub-grid scale water ice clouds
-      real,intent(out) :: totcloudfrac(ngrid)
-      logical,intent(in) :: clearsky
+      LOGICAL,INTENT(IN) :: clearsky
+      REAL,INTENT(IN) :: totcloudfrac(ngrid)
 
 c
@@ -292 +301 @@
          iaer_dust_submicron=0
          iaer_h2o_ice=0
+         iaer_stormdust_doubleq=0
 
          aer_count=0
@@ -326 +336 @@
            enddo
          endif
+         if (rdstorm.AND.active) then
+           do iaer=1,naerkind
+             if (name_iaer(iaer).eq."stormdust_doubleq") then
+               iaer_stormdust_doubleq = iaer
+               aer_count = aer_count + 1
+             endif
+           enddo
+         end if
 
 c        Check that we identified all tracers:
@@ -372 +390 @@
 c     Updating aerosol size distributions:
       CALL updatereffrad(ngrid,nlayer,
-     &                rdust,rice,nuice,
+     &                rdust,rstormdust,rice,nuice,
      &                reffrad,nueffrad,
      &                pq,tauscaling,tau,pplay)
@@ -385 +403 @@
 c     Computing aerosol optical depth in each layer:
       CALL aeropacity(ngrid,nlayer,nq,zday,pplay,pplev,ls,
-     &     pq,tauscaling,tauref,tau,taucloudtes,aerosol,dsodust,reffrad,
-     &     nueffrad,QREFvis3d,QREFir3d,omegaREFvis3d,omegaREFir3d,
-     &     clearsky,totcloudfrac)
+     &    pq,tauscaling,tauref,tau,taucloudtes,aerosol,dsodust,reffrad,
+     &    QREFvis3d,QREFir3d,omegaREFir3d,
+     &    totstormfract,clearatm,dsords,
+     &    clearsky,totcloudfrac)
 
 c     Starting loop on sub-domain