Changeset 5104 for LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/seasalt.f90

Timestamp:

Jul 23, 2024, 5:57:06 PM (4 months ago)

Author:

abarral

Message:

Replace 1DUTILS.h by module lmdz_1dutils.f90
Replace 1DConv.h by module lmdz_old_1dconv.f90 (it's only used by old_* files)
Convert *.F in DUST to *.f90

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/seasalt.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/seasalt.F) (1 diff)

LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/seasalt.f90

@@ -1 @@
-c       This SUBROUTINE estimateis Sea Salt emission fluxes over
-c       Oceanic surfaces.
-c 
-      SUBROUTINE seasalt(v_10m, u_10m, pct_ocean, lmt_sea_salt)
-    
-      USE dimphy
-      IMPLICIT NONE 
-c
-      INCLUDE "dimensions.h"
-      INCLUDE "chem.h"
-      INCLUDE "chem_spla.h"
-      INCLUDE "YOMCST.h"
-      INCLUDE "YOECUMF.h"
-c
-      INTEGER i, bin                 !local variables
-      REAL pct_ocean(klon)           !hfraction of Ocean in each grid
-      REAL v_10m(klon), u_10m(klon)  !V&H components of wind @10 m 
-      REAL w_speed_10m(klon)         !wind speed at 10m from surface
-      REAL lmt_sea_salt(klon,ss_bins)!sea salt emission flux - mg/m2/s
-      REAL sea_salt_flux(ss_bins)    !sea salt emission flux per unit wind speed
+! This SUBROUTINE estimateis Sea Salt emission fluxes over
+! Oceanic surfaces.
+!
+SUBROUTINE seasalt(v_10m, u_10m, pct_ocean, lmt_sea_salt)
 
-      REAL wind, ocean
-c
-c------Sea salt emission fluxes for each size bin calculated 
-c------based on on parameterisation of Gong et al. (1997).
-c------Fluxes of sea salt for each size bin are given in mg/m^2/sec
-c------at wind speed of 1 m/s at 10m height (at 80% RH).
-c------Fluxes at various wind speeds (@10 m from sea 
-c------surfaces are estimated using relationship: F=flux*U_10^3.14
-c
-cnhl for size bin of 0.03-0.5 and 0.5-20
-      DATA sea_salt_flux/4.5E-09,8.7E-7/
+  USE dimphy
+  IMPLICIT NONE
+  !
+  INCLUDE "dimensions.h"
+  INCLUDE "chem.h"
+  INCLUDE "chem_spla.h"
+  INCLUDE "YOMCST.h"
+  INCLUDE "YOECUMF.h"
+  !
+  INTEGER :: i, bin                 !local variables
+  REAL :: pct_ocean(klon)           !hfraction of Ocean in each grid
+  REAL :: v_10m(klon), u_10m(klon)  !V&H components of wind @10 m
+  REAL :: w_speed_10m(klon)         !wind speed at 10m from surface
+  REAL :: lmt_sea_salt(klon, ss_bins)!sea salt emission flux - mg/m2/s
+  REAL :: sea_salt_flux(ss_bins)    !sea salt emission flux per unit wind speed
 
-      DO i=1, klon
-      w_speed_10m(i)= (v_10m(i)**2.0+u_10m(i)**2.0)**0.5
-      ENDDO
-c
-      DO bin=1,ss_bins
-      wind=0.0
-      ocean=0.0
-      DO i=1, klon
-      lmt_sea_salt(i,bin)=sea_salt_flux(bin)*(w_speed_10m(i)**3.41)
-     . *pct_ocean(i)*1.e-4*1.e-3                       !g/cm2/s
-      wind=wind+w_speed_10m(i)
-      ocean=ocean+pct_ocean(i)
-      ENDDO
-!      print *,'Sea Salt flux = ',sea_salt_flux(bin)
-      ENDDO
-!      print *,'SUM OF WIND = ',wind
-!      print *,'SUM OF OCEAN SURFACE = ',ocean
-      RETURN    
-      END
+  REAL :: wind, ocean
+  !
+  !------Sea salt emission fluxes for each size bin calculated
+  !------based on on parameterisation of Gong et al. (1997).
+  !------Fluxes of sea salt for each size bin are given in mg/m^2/sec
+  !------at wind speed of 1 m/s at 10m height (at 80% RH).
+  !------Fluxes at various wind speeds (@10 m from sea
+  !------surfaces are estimated using relationship: F=flux*U_10^3.14
+  !
+  !nhl for size bin of 0.03-0.5 and 0.5-20
+  DATA sea_salt_flux/4.5E-09, 8.7E-7/
+
+  DO i = 1, klon
+    w_speed_10m(i) = (v_10m(i)**2.0 + u_10m(i)**2.0)**0.5
+  ENDDO
+  !
+  DO bin = 1, ss_bins
+    wind = 0.0
+    ocean = 0.0
+    DO i = 1, klon
+      lmt_sea_salt(i, bin) = sea_salt_flux(bin) * (w_speed_10m(i)**3.41) &
+              * pct_ocean(i) * 1.e-4 * 1.e-3                       !g/cm2/s
+      wind = wind + w_speed_10m(i)
+      ocean = ocean + pct_ocean(i)
+    ENDDO
+    ! print *,'Sea Salt flux = ',sea_salt_flux(bin)
+  ENDDO
+  ! print *,'SUM OF WIND = ',wind
+  ! print *,'SUM OF OCEAN SURFACE = ',ocean
+  RETURN
+END SUBROUTINE seasalt