Changeset 2704 for LMDZ5/trunk/libf/phylmd/StratAer/traccoag_mod.F90

Timestamp:

Nov 20, 2016, 2:15:32 PM (8 years ago)

Author:

oboucher

Message:

This revision concerns the StratAer? module and should not impact the rest of LMDz
Bug correction in interp_sulf_input.F90
Update of miecalc_aer.F90 and traccoag_mod.F90
Phytrac tracers are now dealt with in XIOS through the Fortran interface with minimal input in the xml
Making tracer groups in DefLists? for StratAer?

File:

• LMDZ5/trunk/libf/phylmd/StratAer/traccoag_mod.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/StratAer/traccoag_mod.F90

@@ -73 +73 @@
     REAL,PARAMETER    :: xlat_sai=0.0           ! latitude of SAI in degree
     REAL,PARAMETER    :: xlon_sai=120.35        ! longitude of SAI in degree
-
-!--be careful - this needs to be changed with resolution - here for 96x95 - seems to work for 48x36 as well
-    REAL,PARAMETER    :: dlat=0.9474   ! d latitude in degree
-    REAL,PARAMETER    :: dlon=1.875    ! d longitude in degree
 
 !--other local variables
@@ -95 +91 @@
     REAL                                   :: zdz                 ! thickness of atm. model layer in m
     REAL,DIMENSION(klon,klev)              :: dens_aer            ! density of aerosol particles [kg/m3 aerosol] with default H2SO4 mass fraction
+    REAL                                   :: dlat, dlon          ! d latitude and d longitude of grid in degree
 
     IF (is_mpi_root) THEN
       PRINT *,'in traccoag: date from phys_cal_mod =',year_cur,'-',mth_cur,'-',day_cur,'-',hour
     ENDIF
+
+    dlat=180./2./FLOAT(nbp_lat)   ! d latitude in degree
+    dlon=360./2./FLOAT(nbp_lon)   ! d longitude in degree
 
     DO it=1, nbtr_bin
@@ -168 +168 @@
         DO i=1,klon
           !Pinatubo eruption at 15.14N, 120.35E
-          IF  ( xlat(i).GT.xlat_vol-dlat .AND. xlat(i).LT.xlat_vol+dlat .AND. &
-                xlon(i).GT.xlon_vol-dlon .AND. xlon(i).LT.xlon_vol+dlon ) THEN
+          IF  ( xlat(i).GE.xlat_vol-dlat .AND. xlat(i).LT.xlat_vol+dlat .AND. &
+                xlon(i).GE.xlon_vol-dlon .AND. xlon(i).LT.xlon_vol+dlon ) THEN
 !         compute altLMDz
             altLMDz(:)=0.0
@@ -209 +209 @@
 !       IF  ((mth_cur==4 .AND. &
 !       SAI continuous emission o
-        IF  ( xlat(i).GT.xlat_sai-dlat .AND. xlat(i).LT.xlat_sai+dlat .AND. &
-          &   xlon(i).GT.xlon_sai-dlon .AND. xlon(i).LT.xlon_sai+dlon ) THEN
+        IF  ( xlat(i).GE.xlat_sai-dlat .AND. xlat(i).LT.xlat_sai+dlat .AND. &
+          &   xlon(i).GE.xlon_sai-dlon .AND. xlon(i).LT.xlon_sai+dlon ) THEN
 !         compute altLMDz
           altLMDz(:)=0.0