Changeset d8107c2 in flexpart.git
- Timestamp:
- Jan 13, 2016, 5:02:07 PM (8 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 41d8574
- Parents:
- f75967d
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/netcdf_output_mod.f90
rcbad0f1 rd8107c2 244 244 integer :: nspecDimID, npointDimID, nageclassDimID, ncharDimID, pointspecDimID 245 245 integer :: tID, lonID, latID, levID, poleID, lageID, oroID 246 integer :: volID, areaID 246 247 integer :: rellng1ID, rellng2ID, rellat1ID, rellat2ID, relzz1ID, relzz2ID 247 248 integer :: relcomID, relkindzID, relstartID, relendID, relpartID, relxmassID … … 267 268 ! otherwise give an error which can look confusing). 268 269 ! ********************************************************************* 269 ! open(newunit=test_unit,file=trim(path(2)(1:length(2)))//'test_dir.txt',status='replace',&270 ! &err=100)271 ! close (test_unit, status='delete')272 270 open(unit=unittmp,file=trim(path(2)(1:length(2)))//'test_dir.txt',status='replace',& 273 271 &err=100) … … 341 339 call nf90_err(nf90_def_dim(ncid, 'numpoint', numpoint, npointDimID)) 342 340 341 343 342 ! create variables 344 343 !************************* … … 378 377 call nf90_err(nf90_put_att(ncid, levID, 'standard_name', 'height')) 379 378 call nf90_err(nf90_put_att(ncid, levID, 'long_name', 'height above ground')) 379 380 ! volume 381 call nf90_err(nf90_def_var(ncid, 'volume', nf90_float, (/ lonDimID, latDimID, levDimID /), & 382 volID)) 383 384 ! area 385 call nf90_err(nf90_def_var(ncid, 'area', nf90_float, (/ lonDimID, latDimID /), & 386 areaID)) 387 380 388 381 389 if (write_releases.eqv..true.) then … … 612 620 ! levels 613 621 call nf90_err(nf90_put_var(ncid, levID, outheight(1:numzgrid))) 614 622 623 ! volume 624 call nf90_err(nf90_put_var(ncid, volID, volume(:,:,:))) 625 626 ! area 627 call nf90_err(nf90_put_var(ncid, areaID, area(:,:))) 628 615 629 if (write_releases.eqv..true.) then 616 630 ! release point information … … 889 903 if (mpi_mode.gt.0) then 890 904 do l=1,nclassunc 891 auxgrid(l)=drygridunc (ix,jy,ks,kp,l,nage)905 auxgrid(l)=drygridunc0(ix,jy,ks,kp,l,nage) 892 906 end do 893 907 else
Note: See TracChangeset
for help on using the changeset viewer.