Changeset 3524 for palm/trunk/SOURCE/chemistry_model_mod.f90

Timestamp:

Nov 14, 2018 1:36:44 PM (5 years ago)

Author:

raasch

Message:

unused variables removed, missing working precision added, missing preprocessor directives added, bugfix concerning allocation of t_surf_wall_v in nopointer case, declaration statements rearranged to avoid compile time errors, mpi_abort arguments replaced to avoid compile errors

File:

• palm/trunk/SOURCE/chemistry_model_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! working precision added to make code Fortran 2008 conform
+! 
+! 3458 2018-10-30 14:51:23Z kanani
 ! from chemistry branch r3443, banzhafs, basit:
 ! replace surf_lsm_h%qv1(m) by q(k,j,i) for mixing ratio in chem_depo
@@ -4780 +4783 @@
 
          ! vapour pressure deficit influence
-         F_vpd = min(1.,((1.-F_min(lu))*(vpd_min(lu)-vpd)/(vpd_min(lu)-vpd_max(lu)) + F_min(lu) ))
+         F_vpd = MIN( 1.0_wp, ((1.0_wp-F_min(lu))*(vpd_min(lu)-vpd)/(vpd_min(lu)-vpd_max(lu)) + F_min(lu) ) )
          F_vpd = max(F_vpd,F_min(lu))
 
@@ -4861 +4864 @@
 
       ! Compute visible and near-infrared radiation
-      rv=max(0.1,rdu+rdv)
-      rn=max(0.01,rdm+rdn)
+      rv = MAX( 0.1_wp, rdu+rdv )
+      rn = MAX( 0.01_wp, rdm+rdn )
 
       ! Compute ratio between input global radiation and total radiation computed here
-      ratio=min(0.9,glrad/(rv+rn))
+      ratio = MIN( 0.9_wp, glrad/(rv+rn) )
 
       !     Calculate total visible radiation
@@ -4871 +4874 @@
 
       !     Calculate fraction of PAR in the direct beam
-      fv=min(0.99, (0.9-ratio)/0.7)            ! help variable
-      fv=max(0.01,rdu/rv*(1.0-fv**0.6667))     ! fraction of PAR in the direct beam
+      fv = MIN( 0.99_wp, (0.9_wp-ratio)/0.7_wp )            ! help variable
+      fv = MAX( 0.01_wp, rdu/rv*(1.0_wp-fv**0.6667_wp) )     ! fraction of PAR in the direct beam
 
       ! Compute direct and diffuse parts of PAR