Changeset 2327 for trunk/LMDZ.MARS/util

Timestamp:

May 18, 2020, 6:25:32 PM (5 years ago)

Author:

abierjon

Message:

Mars GCM:
Some improvements in aeroptical.F90 :
1) add the possibility to input a diagfi_P.nc file (with pressure as altitude coordinates) ;
2) replace all "title" attributes into "long_name" attributes in accordance with ticket #48

AB

File:

• trunk/LMDZ.MARS/util/aeroptical.F90 (modified) (18 diffs)

trunk/LMDZ.MARS/util/aeroptical.F90

@@ -74 +74 @@
 character(len=16) :: tmpvarname                  ! temporarily stores a variable name
 real,dimension(:,:,:,:,:),allocatable :: opa_aer ! Aerosols opacities [1/km]
-real :: missval                                  ! Value to put in outfile when the reff is out of bounds
-PARAMETER(missval=1E+20)
+real,dimension(:),allocatable :: missval         ! Value to put in outfile when the reff is out of bounds
+                                                 ! or when there is a mising value in input file
 
 
@@ -179 +179 @@
 allocate(mmr(naerkind,lonlen,latlen,altlen,timelen))
 allocate(reff(naerkind,lonlen,latlen,altlen,timelen))
+
+! Initialize missing value
+allocate(missval(naerkind))
+missval(:)=1.e+20
 
 ! Initialize aerok to .true. for all aerosols
@@ -201 +205 @@
   call status_check(ierr,error_text)
   write(*,*) "Dust mass mixing ratio loaded from "//trim(gcmfile)
+  ! Get missing value
+  ierr=nf90_get_att(gcmfid,tmpvarid,"missing_value",missval(1))
+  if (ierr.ne.nf90_noerr) then
+    missval(1) = 1.e+20
+  endif
 endif
 
@@ -238 +247 @@
   call status_check(ierr,error_text)
   write(*,*) "Water ice mass mixing ratio loaded from "//trim(gcmfile)
+  ! Get missing value
+  ierr=nf90_get_att(gcmfid,tmpvarid,"missing_value",missval(2))
+  if (ierr.ne.nf90_noerr) then
+    missval(2) = 1.e+20
+  endif
 endif
 
@@ -532 +546 @@
     ierr=nf90_put_att(outfid,tmpvaridout,"units","opacity/km")
     ierr=nf90_put_att(outfid,tmpvaridout,"refwavelength",wvl_val)
-    ierr=nf90_put_att(outfid,tmpvaridout,"missing_value",missval)
-    write(*,*)"with missing value = ",missval
+    ierr=nf90_put_att(outfid,tmpvaridout,"missing_value",missval(iaer))
+    write(*,*)"with missing value = ",missval(iaer)
     
     ! End netcdf define mode
@@ -557 +571 @@
 
            ! NB: this should especially handle cases when reff=0
-           opa_aer(iaer,ilon,ilat,ialt,it)=missval
+           opa_aer(iaer,ilon,ilat,ialt,it)=missval(iaer)
+           
+         else if (mmr(iaer,ilon,ilat,ialt,it).eq.missval(iaer)) then
+           ! if there is a missing value in input file
+           opa_aer(iaer,ilon,ilat,ialt,it)=missval(iaer)
+           
          else
            if (radiusdyn(isize).eq.reff(iaer,ilon,ilat,ialt,it)) then
@@ -1049 +1068 @@
 integer, intent(in) :: lonlen ! # of grid points along longitude
 integer, intent(in) :: latlen ! # of grid points along latitude
-integer, intent(in) :: altlen ! # of grid points along latitude
+integer, intent(in) :: altlen ! # of grid points along altitude
 integer, intent(in) :: GCM_layers ! # of GCM atmospheric layers
 integer, intent(in) :: outfid ! NetCDF output file ID
@@ -1204 +1223 @@
 
 !==============================================================================
-! 2.2. Hybrid coordinates ap() , bp(), aps() and bps()
+! 2.1 Hybrid coordinates ap() , bp(), aps() and bps()
 !==============================================================================
 if(hybrid) then 
@@ -1216 +1235 @@
   endif
   ! Write the attributes
-  ierr=nf90_put_att(outfid,tmpvarid,"title","hybrid pressure at midlayers")
+  ierr=nf90_put_att(outfid,tmpvarid,"long_name","hybrid pressure at midlayers")
   ierr=nf90_put_att(outfid,tmpvarid,"units"," ")
   ! End netcdf define mode
@@ -1236 +1255 @@
   endif
   ! Write the attributes
-  ierr=nf90_put_att(outfid,tmpvarid,"title","hybrid sigma at midlayers")
+  ierr=nf90_put_att(outfid,tmpvarid,"long_name","hybrid sigma at midlayers")
   ierr=nf90_put_att(outfid,tmpvarid,"units"," ")
   ! End netcdf define mode
@@ -1258 +1277 @@
     endif
     ! Write the attributes
-    ierr=nf90_put_att(outfid,tmpvarid,"title","hybrid sigma at interlayers")
+    ierr=nf90_put_att(outfid,tmpvarid,"long_name","hybrid sigma at interlayers")
     ierr=nf90_put_att(outfid,tmpvarid,"units"," ")
     ! End netcdf define mode
@@ -1279 +1298 @@
     endif
     ! Write the attributes
-    ierr=nf90_put_att(outfid,tmpvarid,"title","hybrid sigma at interlayers")
+    ierr=nf90_put_att(outfid,tmpvarid,"long_name","hybrid sigma at interlayers")
     ierr=nf90_put_att(outfid,tmpvarid,"units"," ")
     ! End netcdf define mode
@@ -1301 +1320 @@
   endif
   ! Write the attributes
-  ierr=nf90_put_att(outfid,tmpvarid,"title","sigma at midlayers")
+  ierr=nf90_put_att(outfid,tmpvarid,"long_name","sigma at midlayers")
   ierr=nf90_put_att(outfid,tmpvarid,"units"," ")
   ! End netcdf define mode
@@ -1314 +1333 @@
 
 !==============================================================================
-! 2.2. aire() and phisinit()
+! 2.2 aire() and phisinit()
 !==============================================================================
 
@@ -1327 +1346 @@
   endif
   ! Write the attributes
-  ierr=nf90_put_att(outfid,tmpvarid,"title","Mesh area")
+  ierr=nf90_put_att(outfid,tmpvarid,"long_name","Mesh area")
   ierr=nf90_put_att(outfid,tmpvarid,"units","m2")
   ! End netcdf define mode
@@ -1339 +1358 @@
 endif ! of if (area)
 
-IF (phis) THEN
+if (phis) then
 
   ! Switch to netcdf define mode
@@ -1349 +1368 @@
   endif
   ! Write the attributes
-  ierr=nf90_put_att(outfid,tmpvarid,"title","Ground level geopotential")
+  ierr=nf90_put_att(outfid,tmpvarid,"long_name","Ground level geopotential")
   ierr=nf90_put_att(outfid,tmpvarid,"units"," ")
   ! End netcdf define mode
@@ -1360 +1379 @@
   endif
 
-ENDIF ! of IF (phis)
+endif ! of if (phis)