Changeset 6a678e3 in flexpart.git for src/netcdf_output_mod.f90

Timestamp:

Jan 25, 2016, 2:17:34 PM (8 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

26f6039

Parents:

fddc6ec

Message:

Added option to use double precision for calculating the deposition fields

File:

• src/netcdf_output_mod.f90 (modified) (7 diffs)

src/netcdf_output_mod.f90

@@ -36 +36 @@
   !    switch write_releases
   !  - additional updates for FLEXPART 9.x
+  ! 
+  ! ESO 2016
+  !  - Deposition fields can be calculated in double precision, see variable
+  !    'dep_prec' in par_mod
   !*****************************************************************************
 
@@ -48 +52 @@
                        wetgrid,wetgridsigma,drygrid,drygridsigma,grid,gridsigma,&
                        area,arean,volumen, orooutn
-  use par_mod,   only: dp, maxspec, maxreceptor, nclassunc,&
+  use par_mod,   only: dep_prec, sp, dp, maxspec, maxreceptor, nclassunc,&
                        unitoutrecept,unitoutreceptppt, nxmax,unittmp
   use com_mod,   only: path,length,ldirect,ibdate,ibtime,iedate,ietime, &
@@ -70 +74 @@
                        nested_output, ipout, surf_only, linit_cond, &
                        flexversion,mpi_mode
+
+  use mean_mod
 
   implicit none
@@ -723 +729 @@
   integer, intent(in) :: itime
   real, intent(in)    :: outnum
-  real, intent(out)   :: gridtotalunc,wetgridtotalunc,drygridtotalunc
+  real(dep_prec),intent(out):: wetgridtotalunc,drygridtotalunc
+  real, intent(out)   :: gridtotalunc
   real                :: densityoutrecept(maxreceptor)
   integer             :: ncid,kp,ks,kz,ix,jy,iix,jjy,kzz,kzzm1,ngrid
-  integer             :: nage,i,l, jj
+  integer             :: nage,i,l,jj
   real                :: tot_mu(maxspec,maxpointspec_act)
   real                :: halfheight,dz,dz1,dz2
   real                :: xl,yl,xlrot,ylrot,zagnd,zagndprev
-  real                :: auxgrid(nclassunc),gridtotal,gridsigmatotal
-  real                :: wetgridtotal,wetgridsigmatotal
-  real                :: drygridtotal,drygridsigmatotal
+  real(dep_prec)      :: auxgrid(nclassunc)
+  real(dep_prec)      :: gridtotal,gridsigmatotal
+  real(dep_prec)      :: wetgridtotal,wetgridsigmatotal
+  real(dep_prec)      :: drygridtotal,drygridsigmatotal
+  ! real(sp)            :: gridtotal,gridsigmatotal
+  ! real(sp)            :: wetgridtotal,wetgridsigmatotal
+  ! real(sp)            :: drygridtotal,drygridsigmatotal
 
   real, parameter     :: weightair=28.97
@@ -889 +900 @@
                    wetgridsigma(ix,jy),nclassunc)
               ! Multiply by number of classes to get total concentration
-              wetgrid(ix,jy)=wetgrid(ix,jy)*real(nclassunc)
+              wetgrid(ix,jy)=wetgrid(ix,jy)*real(nclassunc,kind=dep_prec)
               wetgridtotal=wetgridtotal+wetgrid(ix,jy)
               ! Calculate standard deviation of the mean
               wetgridsigma(ix,jy)= &
                    wetgridsigma(ix,jy)* &
-                   sqrt(real(nclassunc))
+                   sqrt(real(nclassunc,kind=dep_prec))
               wetgridsigmatotal=wetgridsigmatotal+ &
                    wetgridsigma(ix,jy)
@@ -1061 +1072 @@
   integer, intent(in) :: itime
   real, intent(in)    :: outnum
-  real, intent(out)   :: gridtotalunc,wetgridtotalunc,drygridtotalunc
+  real(sp), intent(out)   :: gridtotalunc
+  real(dep_prec), intent(out)   :: wetgridtotalunc,drygridtotalunc
 
   print*,'Netcdf output for surface only not yet implemented'
@@ -1117 +1129 @@
   real                :: halfheight,dz,dz1,dz2
   real                :: xl,yl,xlrot,ylrot,zagnd,zagndprev
-  real                :: auxgrid(nclassunc),gridtotal
+  real(dep_prec)      :: auxgrid(nclassunc)
+  real                :: gridtotal
   real, parameter     :: weightair=28.97