Changeset 2435

Timestamp:

Nov 23, 2020, 11:50:01 AM (4 years ago)

Author:

abierjon

Message:

Mars GCM:
Add stormdust in the aerosol opacities computed by aeroptical.F90
+ add a description of aeroptical in util/README

AB

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• util/README (modified) (2 diffs)
• util/aeroptical.F90 (modified) (8 diffs)

trunk/LMDZ.MARS/README

@@ -3235 +3235 @@
 - extend size of long_name character string (ticket #57)
 - remove the use of #ifdef NC_DOUBLE (ticket #67)
+ 
+== 23/11/2020 == AB
+Add stormdust in the aerosol opacities computed by aeroptical.F90
++ add a description of aeroptical in util/README

trunk/LMDZ.MARS/util/README

@@ -22 +22 @@
 This program is used to concatenate consecutive
 output files ("diagfi", and even "stats") for a selection of 
-variable, in order to obtain one single big file. The time simension can
+variable, in order to obtain one single big file. The time dimension can
 be "sols" or "Ls" (in that last case, Ls will not be linear, and
 software like Grads may not be able to use it. To have a linear "Ls"
@@ -187 +187 @@
 extract.e < extract.points.def
 
-  
+--------------------------------------------------------------------
+9) aeroptical.e : compute aerosol opacities at a given wavelength
+--------------------------------------------------------------------
+
+Program to compute opacities [1/km] of all aerosols from dust, water ice
+and stormdust present in the input file, at a wavelength given by the user.
+Computation is made from the aerosol mass mixing ratios and effective radius,
+associated with air density (rho) and files containing the aerosol 
+optical properties (generally present in the GCM datadir/ directory).
+
+NB : water ice effective radius is known to be a bit inaccurate in the GCM,
+so reffwice can be prescribed by the user.
+
+input : diagfi.nc  / concat.nc / stats.nc kind of files
+
+output file is : name_of_input_file_OPA.nc
+aeroptical.e < aeroptical.def

trunk/LMDZ.MARS/util/aeroptical.F90

@@ -8 +8 @@
 ! This program takes as inputs a GCM output file (diagfi,stats,concat) that
 ! contains:
-!      - the Mass Mixing Ratios of dust (dustq) and water ice (h2o_ice)
-!      - their effective radius (reffdust, reffice(*))
+!      - the Mass Mixing Ratios of dust (dustq), water ice (h2o_ice)
+!                                 & stormdust (rdsdustq)
+!      - their effective radius (reffdust, reffice(*), reffstormdust)
 !      - the atmospheric density (rho)
 ! as well as a wavelength to calibrate the opacities, and the directory
 ! containing the ASCII files of the aerosols' optical properties.
 !
-! It outputs a netcdf file containing the opacities [1/km] of the dust and
-! water ice aerosols calibrated at the prescribed wavelength.
+! It outputs a netcdf file containing the opacities [1/km] of the dust,
+! water ice and stormdust aerosols calibrated at the prescribed wavelength.
 !
 !
@@ -56 +57 @@
 real,dimension(:,:,:,:),allocatable :: rho    ! atmospheric density [kg.m-3]
 integer :: naerkind                           ! nb of aerosols to consider
-integer :: iaer                               ! aerosol kind index (1=dust, 2=water ice) ! + NEW AER?
+integer :: iaer                               ! aerosol kind index (1=dust, 2=water ice, 3=stormdust) ! + NEW AER?
 integer :: ilon,ilat,ialt,it                  ! Loop indices for coordinates
 
@@ -174 +175 @@
 
 ! Number of aerosols
-naerkind = 2 ! dust & water ice ! + NEW AER?
+naerkind = 3 ! dust, water ice, stormdust ! + NEW AER?
 
 ! Length allocation of the variables
@@ -288 +289 @@
 
 !==========================================================================
-! --> 1.2.3 + NEW AER? 
-
-!==========================================================================
-! --> 1.2.3 ATMOSPHERIC DENSITY
+! --> 1.2.3 STORMDUST
+
+! STORMDUST MASS MIXING RATIO
+! Check that the GCM file contains that variable
+ierr=nf90_inq_varid(gcmfid,"rdsdustq",tmpvarid)
+error_text="Failed to find variable rdsdustq in "//trim(gcmfile)&
+            //" We'll skip the stormdust aerosol."
+if (ierr.ne.nf90_noerr) then
+  write(*,*)trim(error_text)
+  aerok(3)=.false.
+else
+  ! Load datasets
+  ierr=NF90_GET_VAR(gcmfid,tmpvarid,mmr(3,:,:,:,:))
+  error_text="Failed to load stormdust mass mixing ratio"
+  call status_check(ierr,error_text)
+  write(*,*) "Stormdust mass mixing ratio loaded from "//trim(gcmfile)
+  ! Get missing value
+  ierr=nf90_get_att(gcmfid,tmpvarid,"missing_value",missval(3))
+  if (ierr.ne.nf90_noerr) then
+    missval(3) = 1.e+20
+  endif
+endif
+
+! STORMDUST EFFECTIVE RADIUS
+if (aerok(3)) then
+  ! Check that the GCM file contains that variable
+  ierr=nf90_inq_varid(gcmfid,"reffstormdust",tmpvarid)
+  error_text="Failed to find variable reffstormdust in "//trim(gcmfile)&
+              //" We'll skip the dust aerosol."
+  if (ierr.ne.nf90_noerr) then
+    write(*,*)trim(error_text)
+    aerok(3)=.false.
+  else
+    ! Load datasets
+    ierr=NF90_GET_VAR(gcmfid,tmpvarid,reff(3,:,:,:,:))
+    error_text="Failed to load stormdust effective radius"
+    call status_check(ierr,error_text)
+    write(*,*) "Stormdust effective radius loaded from "//trim(gcmfile)
+  endif
+endif 
+
+
+!==========================================================================
+! --> 1.2.4 + NEW AER? 
+
+
+!==========================================================================
+! --> 1.2.5 ATMOSPHERIC DENSITY
 
 ! Check that the GCM file contains that variable
@@ -319 +364 @@
 rho_aer(1)=2500. ! dust
 rho_aer(2)=920.  ! water ice
+rho_aer(3)=2500. ! stormdust
 ! + NEW AER?
 
@@ -344 +390 @@
         ! Water ice file
         optpropfile = "optprop_icevis_n30.dat"
+      CASE(3) ! STORMDUST
+        ! Dust file
+        optpropfile = "optprop_dustvis_TM_n50.dat"
 ! + NEW AER?
       END SELECT ! iaer
@@ -357 +406 @@
         ! Water ice file
         optpropfile = "optprop_iceir_n30.dat"
+      CASE(3) ! STORMDUST
+        ! Dust file
+        optpropfile = "optprop_dustir_n50.dat"
 ! + NEW AER?
       END SELECT ! iaer
@@ -540 +592 @@
       ! Write the attributes
       ierr=nf90_put_att(outfid,tmpvaridout,"long_name","Water ice extinction opacity at reference wavelength")
+      
+    CASE(3) ! STORMDUST
+      ! Insert the definition of the variable
+      tmpvarname="opastormdust"
+      ierr=NF90_DEF_VAR(outfid,trim(tmpvarname),nf90_float,varshape,tmpvaridout)
+      error_text="ERROR: Couldn't create "//trim(tmpvarname)//" in the outfile"
+      call status_check(ierr,error_text)
+      write(*,*) trim(tmpvarname)//" has been created in the outfile"
+      
+      ! Write the attributes
+      ierr=nf90_put_att(outfid,tmpvaridout,"long_name","Stormdust extinction opacity at reference wavelength")
 
 ! + NEW AER?