Changeset 2830 for trunk/LMDZ.MARS/util/aeropt_mod.F90

Timestamp:

Nov 23, 2022, 4:09:32 PM (2 years ago)

Author:

abierjon

Message:

Mars GCM:
The program util/aeroptical.F90 can now compute column-integrated optical depths
of the aerosols, in addition to the opacity profiles. The user has to specify it
('yes' or 'no') in aeroptical.def
Detailed changes :

• update of the aeroptical.def file to ask the user if the column optical depth should be computed (non-retrocompatible change)
• add in the init2 subroutine a computation of the layers' height delta_z, with adaptable method depending on the availability of some variables in the input file and on wether or not the input file has been zrecasted before. Preliminary validation shows that these different methods yield a +/-3% error on the output column optical depth tau_[aer]. The delta_z variable is also written in the output file.
• add a log file for warnings in the interpolation subroutine in aeropt_mod.F90 + add some comments in the code
• add the ouput "zzlev" (= interlayer altitudes) in libf/phymars/physiq_mod.F so that it can be directly used by aeroptical.F90 to compute delta_z

AB

File:

• trunk/LMDZ.MARS/util/aeropt_mod.F90 (modified) (6 diffs)

trunk/LMDZ.MARS/util/aeropt_mod.F90

@@ -12 +12 @@
  real,dimension(:,:),save,allocatable :: Qext    ! Extinction efficiency coefficient [/]
  real,dimension(:,:),save,allocatable :: omeg    ! Single scattering albedo [/]
+ 
+ integer,save :: aeroptlogfileID = 100           ! Log file's unit ID
 
  CONTAINS 
@@ -46 +48 @@
   integer :: read_ok                         ! to know if the line is well read
   integer :: iwvl,isize                      ! Wavelength and Particle size loop indices
+  
+  !==============================================================================
+  ! aeropt_mod module variables used here:
+  !==============================================================================
+  ! nwvl,nsize,wvl,radiusdyn,Qext,omeg
 
   !==========================================================================
@@ -228 +235 @@
   real,dimension(2) :: reffQext         ! Qext after reff interpolation
   real,dimension(2) :: reffomeg         ! omeg after reff interpolation
+  
+  !==============================================================================
+  ! aeropt_mod module variables used here:
+  !==============================================================================
+  ! nwvl,nsize,wvl,radiusdyn,Qext,omeg,aeroptlogfileID
 
   !==========================================================================
@@ -262 +274 @@
   ENDDO
 
-  if ((radiusdyn(nsize).lt.reff_val).or.(radiusdyn(1).gt.reff_val)) then
-    write(*,*) "WARNING: the effective radius (",reff_val,"m) is out of the bounds"
-    write(*,*) "of the optical properties' tables"
-    write(*,*) "(reff_inf=",radiusdyn(1),"m ; reff_sup=",radiusdyn(nsize),"m)"
-    write(*,*) "A missing value will be written for opacity."
-
-    ! NB: this should especially handle cases when reff=0
+  if (reff_val.eq.missval) then
+    ! if the effective radius is a missing value,
+    ! put a missing value for the optical properties
+    Qext_val = missval
+    omeg_val = missval
+!    write(*,*) "reff missing value detected"
+    
+  else if ((radiusdyn(nsize).lt.reff_val).or.(radiusdyn(1).gt.reff_val)) then
+    write(aeroptlogfileID,*) "WARNING: the effective radius (",reff_val,"m) is out of the bounds"
+    write(aeroptlogfileID,*) "of the optical properties' tables"
+    write(aeroptlogfileID,*) "(reff_inf=",radiusdyn(1),"m ; reff_sup=",radiusdyn(nsize),"m)"
+    write(aeroptlogfileID,*) "A missing value will be written for optical properties."
+
     Qext_val = missval
     omeg_val = missval
@@ -333 +351 @@
 
   implicit none
+  
+  !==============================================================================
+  ! aeropt_mod module variables used here:
+  !==============================================================================
+  ! wvl,radiusdyn,Qext,omeg
 
   if (allocated(wvl)) deallocate(wvl)
@@ -342 +365 @@
  
 !*******************************************************************************
+ SUBROUTINE create_logfile(filename)
+  !==============================================================================
+  ! Purpose:
+  ! Deallocate module variables
+  !==============================================================================
+
+  implicit none
+  
+  !==============================================================================
+  ! Arguments:
+  !==============================================================================
+  character(len=128),intent(in) :: filename ! name of the aeropt log file
+  !==============================================================================
+  ! aeropt_mod module variables used here:
+  !==============================================================================
+  ! aeroptlogfileID
+  
+  ! Create output file
+  write(*,*) "Creating aeropt log file: ",trim(filename)
+  write(*,*) "It will overwrite any existing file with the same name."
+  OPEN(UNIT=aeroptlogfileID,FILE=trim(filename),FORM='formatted',STATUS='replace')
+  ! NB: STATUS='replace' will overwrite any existing file with the same name
+  
+ END SUBROUTINE create_logfile 
+ 
+!*******************************************************************************
 
 END MODULE aeropt_mod