Changeset 5924 for LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.f90

Timestamp:

Dec 9, 2025, 3:08:05 PM (8 hours ago)

Author:

lebasn

Message:

StratAer?: remove old params and rename others for coherence (flag_nucl and flag_strat_compo). Use integer instead of boolean to choose compo and nucleation methods.

File:

• LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.f90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.f90

@@ -19 +19 @@
   USE aerophys
   USE yomcst_mod_h
-  USE strataer_local_var_mod, ONLY : flag_newclim_file,flag_verbose_strataer
+  USE strataer_local_var_mod, ONLY : flag_verbose_strataer
 
 IMPLICIT NONE
@@ -125 +125 @@
      
     !--init ncdf variables
-    IF(flag_newclim_file) THEN
-      nc_fname = "ocs_so2_h2so4_annual_lmdz.nc"
-      nc_lat   = "LAT"
-      nc_lon   = "LON"
-    ELSE
-      ! old file for retro compatibility
-      nc_fname = "ocs_so2_annual_lmdz.nc"
-      nc_lat   = "lat"
-      nc_lon   = "lon"
-    ENDIF
-
+    nc_fname = "ocs_so2_h2so4_annual_lmdz.nc"
+    nc_lat   = "LAT"
+    nc_lon   = "LON"
+    
     !--reading emission files
     CALL nf95_open(nc_fname, nf90_nowrite, ncid_in)
@@ -180 +173 @@
     IF(flag_verbose_strataer) print *,'code erreur SO2_lifetime_in=', ncerr, varid
     
-    IF(flag_newclim_file) THEN
-       CALL nf95_inq_varid(ncid_in, "O3", varid)
-       ncerr = nf90_get_var(ncid_in, varid, O3_clim_in)
-       IF(flag_verbose_strataer) print *,'code erreur O3=', ncerr, varid
-       
-       CALL nf95_inq_varid(ncid_in, "H2SO4_LIFET", varid)
-       ncerr = nf90_get_var(ncid_in, varid, H2SO4_lifetime_in)
-       IF(flag_verbose_strataer) print *,'code erreur H2SO4_lifetime_in=', ncerr, varid
-    ENDIF
+    CALL nf95_inq_varid(ncid_in, "O3", varid)
+    ncerr = nf90_get_var(ncid_in, varid, O3_clim_in)
+    IF(flag_verbose_strataer) print *,'code erreur O3=', ncerr, varid
+    
+    CALL nf95_inq_varid(ncid_in, "H2SO4_LIFET", varid)
+    ncerr = nf90_get_var(ncid_in, varid, H2SO4_lifetime_in)
+    IF(flag_verbose_strataer) print *,'code erreur H2SO4_lifetime_in=', ncerr, varid
     
     CALL nf95_close(ncid_in)
@@ -215 +206 @@
       
       ! O3 from 2d model is not tracer, in VMR
-      IF(flag_newclim_file) THEN
-         H2SO4_lifetime_mth(:,j,:) = H2SO4_lifetime_in(:,n_lat+1-j,:,mth_cur)
-         ! new input files
-         O3_clim_mth(:,j,:) = 1.e-6*O3_clim_in(:,n_lat+1-j,:,mth_cur)
-      ELSE
-         H2SO4_lifetime_mth(:,j,:) = 1.e-6
-         O3_clim_mth(:,j,:) = 1.e-6
-      ENDIF
+      H2SO4_lifetime_mth(:,j,:) = H2SO4_lifetime_in(:,n_lat+1-j,:,mth_cur)
+      ! new input files
+      O3_clim_mth(:,j,:) = 1.e-6*O3_clim_in(:,n_lat+1-j,:,mth_cur)
     ENDDO
 
@@ -279 +265 @@
            MAX(0.0,MIN(paprs_glo(i,k),paprs_input(kk))-MAX(paprs_glo(i,k+1),paprs_input(kk+1))) &
            *SO2_lifetime_tmp(i,kk)/(paprs_glo(i,k)-paprs_glo(i,k+1))
-      IF(flag_newclim_file) THEN
-         H2SO4_lifetime_glo(i,k)=H2SO4_lifetime_glo(i,k)+ &
-              MAX(0.0,MIN(paprs_glo(i,k),paprs_input(kk)) &
-              -MAX(paprs_glo(i,k+1),paprs_input(kk+1))) &
-              *H2SO4_lifetime_tmp(i,kk)/(paprs_glo(i,k)-paprs_glo(i,k+1))
-      ENDIF
+      H2SO4_lifetime_glo(i,k)=H2SO4_lifetime_glo(i,k)+ &
+           MAX(0.0,MIN(paprs_glo(i,k),paprs_input(kk)) &
+           -MAX(paprs_glo(i,k+1),paprs_input(kk+1))) &
+           *H2SO4_lifetime_tmp(i,kk)/(paprs_glo(i,k)-paprs_glo(i,k+1))
       
       OCS_clim_glo(i,k)=OCS_clim_glo(i,k)+ &
@@ -292 +276 @@
            MAX(0.0,MIN(paprs_glo(i,k),paprs_input(kk))-MAX(paprs_glo(i,k+1),paprs_input(kk+1))) &
            *SO2_clim_tmp(i,kk)/(paprs_glo(i,k)-paprs_glo(i,k+1))
-      IF(flag_newclim_file) THEN
-         O3_clim_glo(i,k)=O3_clim_glo(i,k)+ &
-              MAX(0.0,MIN(paprs_glo(i,k),paprs_input(kk)) &
-              -MAX(paprs_glo(i,k+1),paprs_input(kk+1))) &
-              *O3_clim_tmp(i,kk)/(paprs_glo(i,k)-paprs_glo(i,k+1))
-      ENDIF
+      O3_clim_glo(i,k)=O3_clim_glo(i,k)+ &
+           MAX(0.0,MIN(paprs_glo(i,k),paprs_input(kk)) &
+           -MAX(paprs_glo(i,k+1),paprs_input(kk+1))) &
+           *O3_clim_tmp(i,kk)/(paprs_glo(i,k)-paprs_glo(i,k+1))
       ENDDO
     ENDDO