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

Timestamp:

Jul 23, 2024, 5:57:06 PM (2 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/deposition.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/deposition.F) (1 diff)

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

@@ -1 @@
-c Subroutine that estimates the Deposition velocities and the depostion
-C for the different tracers
-      SUBROUTINE deposition(vdep_oce,vdep_sic,vdep_ter,vdep_lic,pctsrf,
-     .                      zrho,zdz,pdtphys,RHcl,masse,t_seri,pplay,
-     .                      paprs,lminmax,qmin,qmax,
-     .                        his_ds,source_tr,tr_seri)
+! Subroutine that estimates the Deposition velocities and the depostion
+! for the different tracers
+SUBROUTINE deposition(vdep_oce, vdep_sic, vdep_ter, vdep_lic, pctsrf, &
+        zrho, zdz, pdtphys, RHcl, masse, t_seri, pplay, &
+        paprs, lminmax, qmin, qmax, &
+        his_ds, source_tr, tr_seri)
 
-      USE dimphy
-      USE infotrac
-      USE indice_sol_mod
+  USE dimphy
+  USE infotrac
+  USE indice_sol_mod
 
-      IMPLICIT NONE
+  IMPLICIT NONE
 
-      INCLUDE "dimensions.h"
-      INCLUDE "chem.h"
-      INCLUDE "YOMCST.h"
-      INCLUDE "paramet.h"
+  INCLUDE "dimensions.h"
+  INCLUDE "chem.h"
+  INCLUDE "YOMCST.h"
+  INCLUDE "paramet.h"
 
-c----------------------------- INPUT -----------------------------------
-      LOGICAL lminmax
-      REAL qmin, qmax
-      REAL vdep_oce(nbtr), vdep_sic(nbtr)
-      REAL vdep_ter(nbtr), vdep_lic(nbtr)      
-      REAL pctsrf(klon,nbsrf)
-      REAL zrho(klon,klev)        !Density of air at mid points of Z (kg/m3)
-      REAL zdz(klon,klev)        
-      REAL pdtphys  ! pas d'integration pour la physique (seconde)
-      REAL RHcl(klon,klev)  ! humidite relativen ciel clair
-      REAL t_seri(klon,klev)  ! temperature
-      REAL pplay(klon,klev)  ! pression pour le mileu de chaque couche (en Pa)
-      REAL paprs(klon, klev+1)    !pressure at interface of layers Z (Pa)
-      REAL masse(nbtr)
-                                          
-c----------------------------- OUTPUT ----------------------------------
-      REAL his_ds(klon,nbtr)                                          
-      REAL source_tr(klon,nbtr)
-      REAL tr_seri(klon, klev,nbtr) !conc of tracers
-c--------------------- INTERNAL VARIABLES ------------------------------      
-      INTEGER i, it
-      REAL vdep        !sed. velocity
+  !----------------------------- INPUT -----------------------------------
+  LOGICAL :: lminmax
+  REAL :: qmin, qmax
+  REAL :: vdep_oce(nbtr), vdep_sic(nbtr)
+  REAL :: vdep_ter(nbtr), vdep_lic(nbtr)
+  REAL :: pctsrf(klon, nbsrf)
+  REAL :: zrho(klon, klev)        !Density of air at mid points of Z (kg/m3)
+  REAL :: zdz(klon, klev)
+  REAL :: pdtphys  ! pas d'integration pour la physique (seconde)
+  REAL :: RHcl(klon, klev)  ! humidite relativen ciel clair
+  REAL :: t_seri(klon, klev)  ! temperature
+  REAL :: pplay(klon, klev)  ! pression pour le mileu de chaque couche (en Pa)
+  REAL :: paprs(klon, klev + 1)    !pressure at interface of layers Z (Pa)
+  REAL :: masse(nbtr)
 
-      DO it=1, nbtr
-      DO i=1, klon
-          vdep=vdep_oce(it)*pctsrf(i,is_oce)+
-     .         vdep_sic(it)*pctsrf(i,is_sic)+
-     .         vdep_ter(it)*pctsrf(i,is_ter)+
-     .         vdep_lic(it)*pctsrf(i,is_lic)
-c--Unit: molec/m2/s for it=1 to nbtr-3, mg/m2/s for it=nbtr-2 to nbtr
-          source_tr(i,it)=source_tr(i,it)
-     .                    -vdep*tr_seri(i,1,it)*zrho(i,1)/1.e2
-          his_ds(i,it)=vdep*tr_seri(i,1,it)*zrho(i,1)/1.e2
-     .                 /RNAVO*masse(it)*1.e3               ! mg/m2/s
-      ENDDO
-      ENDDO
-c
-      END
+  !----------------------------- OUTPUT ----------------------------------
+  REAL :: his_ds(klon, nbtr)
+  REAL :: source_tr(klon, nbtr)
+  REAL :: tr_seri(klon, klev, nbtr) !conc of tracers
+  !--------------------- INTERNAL VARIABLES ------------------------------
+  INTEGER :: i, it
+  REAL :: vdep        !sed. velocity
+
+  DO it = 1, nbtr
+    DO i = 1, klon
+      vdep = vdep_oce(it) * pctsrf(i, is_oce) + &
+              vdep_sic(it) * pctsrf(i, is_sic) + &
+              vdep_ter(it) * pctsrf(i, is_ter) + &
+              vdep_lic(it) * pctsrf(i, is_lic)
+      !--Unit: molec/m2/s for it=1 to nbtr-3, mg/m2/s for it=nbtr-2 to nbtr
+      source_tr(i, it) = source_tr(i, it) &
+              - vdep * tr_seri(i, 1, it) * zrho(i, 1) / 1.e2
+      his_ds(i, it) = vdep * tr_seri(i, 1, it) * zrho(i, 1) / 1.e2 &
+              / RNAVO * masse(it) * 1.e3               ! mg/m2/s
+    ENDDO
+  ENDDO
+  !
+END SUBROUTINE deposition