Changes in src/netcdf_output_mod.f90 [4c64400:f28aa0a] in flexpart.git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/netcdf_output_mod.f90
r4c64400 rf28aa0a 19 19 !********************************************************************** 20 20 21 22 21 !***************************************************************************** 23 22 ! * … … 25 24 ! residence time and wet and dry deposition output. * 26 25 ! * 27 ! - writeheader_netcdf generates file including all information previously *26 ! - writeheader_netcdf generates file including all information previously * 28 27 ! 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 * 30 29 ! * 31 30 ! Author: D. Brunner * … … 894 893 gridsigmatotal=0. 895 894 gridtotalunc=0. 896 wetgridtotal=0. _dep_prec897 wetgridsigmatotal=0. _dep_prec898 wetgridtotalunc=0. _dep_prec899 drygridtotal=0. _dep_prec900 drygridsigmatotal=0. _dep_prec901 drygridtotalunc=0. _dep_prec895 wetgridtotal=0. 896 wetgridsigmatotal=0. 897 wetgridtotalunc=0. 898 drygridtotal=0. 899 drygridsigmatotal=0. 900 drygridtotalunc=0. 902 901 903 902 do ks=1,nspec … … 923 922 wetgridsigma(ix,jy),nclassunc) 924 923 ! 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) 926 925 wetgridtotal=wetgridtotal+wetgrid(ix,jy) 927 926 ! Calculate standard deviation of the mean … … 947 946 drygridsigma(ix,jy),nclassunc) 948 947 ! 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) 950 949 drygridtotal=drygridtotal+drygrid(ix,jy) 951 950 ! Calculate standard deviation of the mean 952 951 drygridsigma(ix,jy)= & 953 952 drygridsigma(ix,jy)* & 954 sqrt(real(nclassunc , kind=dep_prec))953 sqrt(real(nclassunc)) 955 954 drygridsigmatotal=drygridsigmatotal+ & 956 955 drygridsigma(ix,jy) … … 1055 1054 if (wetgridtotal.gt.0.) wetgridtotalunc=wetgridsigmatotal/ & 1056 1055 wetgridtotal 1057 if (drygridtotal.gt.0.) drygridtotalunc= real(drygridsigmatotal/ &1058 drygridtotal , kind=dep_prec)1056 if (drygridtotal.gt.0.) drygridtotalunc=drygridsigmatotal/ & 1057 drygridtotal 1059 1058 1060 1059 ! Dump of receptor concentrations … … 1299 1298 wetgridsigma(ix,jy)= & 1300 1299 wetgridsigma(ix,jy)* & 1301 sqrt(real(nclassunc ,kind=dep_prec))1300 sqrt(real(nclassunc)) 1302 1301 endif 1303 1302 … … 1320 1319 drygridsigma(ix,jy)= & 1321 1320 drygridsigma(ix,jy)* & 1322 sqrt(real(nclassunc ,kind=dep_prec))1321 sqrt(real(nclassunc)) 1323 1322 endif 1324 1323
Note: See TracChangeset
for help on using the changeset viewer.