Changes in src/netcdf_output_mod.f90 [4c64400:f28aa0a] in flexpart.git


Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/netcdf_output_mod.f90

    r4c64400 rf28aa0a  
    1919!**********************************************************************
    2020
    21 
    2221  !*****************************************************************************
    2322  !                                                                            *
     
    2524  !  residence time and wet and dry deposition output.                         *
    2625  !                                                                            *
    27   !  - writeheader_netcdf generates file including all information previously  *
     26  !  - writeheader_netcdf generates file including all information previously    *
    2827  !    stored in separate header files                                         *
    29   !  - concoutput_netcdf write concentration output and wet and dry deposition *
     28  !  - concoutput_netcdf write concentration output and wet and dry deposition   *
    3029  !                                                                            *
    3130  !     Author: D. Brunner                                                     *
     
    894893  gridsigmatotal=0.
    895894  gridtotalunc=0.
    896   wetgridtotal=0._dep_prec
    897   wetgridsigmatotal=0._dep_prec
    898   wetgridtotalunc=0._dep_prec
    899   drygridtotal=0._dep_prec
    900   drygridsigmatotal=0._dep_prec
    901   drygridtotalunc=0._dep_prec
     895  wetgridtotal=0.
     896  wetgridsigmatotal=0.
     897  wetgridtotalunc=0.
     898  drygridtotal=0.
     899  drygridsigmatotal=0.
     900  drygridtotalunc=0.
    902901
    903902  do ks=1,nspec
     
    923922                   wetgridsigma(ix,jy),nclassunc)
    924923              ! Multiply by number of classes to get total concentration
    925               wetgrid(ix,jy)=wetgrid(ix,jy)*real(nclassunc,kind=sp)
     924              wetgrid(ix,jy)=wetgrid(ix,jy)*real(nclassunc,kind=dep_prec)
    926925              wetgridtotal=wetgridtotal+wetgrid(ix,jy)
    927926              ! Calculate standard deviation of the mean
     
    947946                   drygridsigma(ix,jy),nclassunc)
    948947              ! Multiply by number of classes to get total concentration
    949               drygrid(ix,jy)=drygrid(ix,jy)*real(nclassunc,kind=sp)
     948              drygrid(ix,jy)=drygrid(ix,jy)*real(nclassunc)
    950949              drygridtotal=drygridtotal+drygrid(ix,jy)
    951950              ! Calculate standard deviation of the mean
    952951              drygridsigma(ix,jy)= &
    953952                   drygridsigma(ix,jy)* &
    954                    sqrt(real(nclassunc, kind=dep_prec))
     953                   sqrt(real(nclassunc))
    955954              drygridsigmatotal=drygridsigmatotal+ &
    956955                   drygridsigma(ix,jy)
     
    10551054  if (wetgridtotal.gt.0.) wetgridtotalunc=wetgridsigmatotal/ &
    10561055       wetgridtotal
    1057   if (drygridtotal.gt.0.) drygridtotalunc=real(drygridsigmatotal/ &
    1058        drygridtotal, kind=dep_prec)
     1056  if (drygridtotal.gt.0.) drygridtotalunc=drygridsigmatotal/ &
     1057       drygridtotal
    10591058
    10601059  ! Dump of receptor concentrations
     
    12991298              wetgridsigma(ix,jy)= &
    13001299                   wetgridsigma(ix,jy)* &
    1301                    sqrt(real(nclassunc,kind=dep_prec))
     1300                   sqrt(real(nclassunc))
    13021301            endif
    13031302
     
    13201319              drygridsigma(ix,jy)= &
    13211320                   drygridsigma(ix,jy)* &
    1322                    sqrt(real(nclassunc,kind=dep_prec))
     1321                   sqrt(real(nclassunc))
    13231322            endif
    13241323
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG