Changeset 1d9681b in flexpart.git


Ignore:
Timestamp:
Feb 17, 2017, 7:38:39 PM (7 years ago)
Author:
Don Morton <Don.Morton@…>
Branches:
fp9.3.1-20161214-nc4
Children:
7e5f1dc
Parents:
d2d8a96
Message:

Continued NC4 for FP development in branch fp9.3.1-20161214-nc4

File:
1 edited

Legend:

Unmodified
Added
Removed
  • flexpart_code/fpmetbinary_mod.F90

    rd2d8a96 r1d9681b  
    4444    USE com_mod
    4545    USE conv_mod
    46     USE par_mod, ONLY : nxmax, nymax, nzmax, nuvzmax, nwzmax, numclass
     46    USE par_mod, ONLY : nxmax, nymax, nzmax, nuvzmax, nwzmax, numclass, maxspec
    4747
    4848    USE netcdf
     
    152152        ! Zeroes out, in local datastructures, the values dumped/loaded
    153153        ! This was written primarily as a testing mechanism.
     154        ! DJM -- 17 February 2017 -- I don't think this routine has been used
     155        ! for anything in recent past.  Might want to consider 86'ing it.
    154156        ! The lines here correspond to READ and WRITE in the dump/load routines
    155157
     
    830832
    831833
     834            dim3dids = (/nxmax_dimid, nymax_dimid, maxspec/)
     835
     836            ncret = nf90_def_var(ncid, 'vdep', NF90_FLOAT, &
     837&                                       dim3dids, ncvarid)
     838            ncret = nf90_def_var_deflate(ncid, ncvarid,   &
     839&                                        shuffle=0,     &
     840&                                        deflate=1,     &
     841&                                        deflate_level=DEF_LEVEL)
     842            ncret = nf90_put_var(ncid, ncvarid, &
     843&                                vdep(0:nxmax-1, 0:nymax-1, 1:maxspec, cm_index))
    832844
    833845
     
    841853            WRITE(iounit) aknew(:)
    842854            WRITE(iounit) bknew(:)
     855
     856            dim1dids = (/numclass/)
     857
     858            ncret = nf90_def_var(ncid, 'z0', NF90_FLOAT, &
     859&                                       dim1dids, ncvarid)
     860            ncret = nf90_def_var_deflate(ncid, ncvarid,   &
     861&                                        shuffle=0,     &
     862&                                        deflate=1,     &
     863&                                        deflate_level=DEF_LEVEL)
     864            ncret = nf90_put_var(ncid, ncvarid, &
     865&                                z0(1:numclass))
     866
     867
     868            dim1dids = (/nwzmax/)
     869
     870            ncret = nf90_def_var(ncid, 'akm', NF90_FLOAT, &
     871&                                       dim1dids, ncvarid)
     872            ncret = nf90_def_var_deflate(ncid, ncvarid,   &
     873&                                        shuffle=0,     &
     874&                                        deflate=1,     &
     875&                                        deflate_level=DEF_LEVEL)
     876            ncret = nf90_put_var(ncid, ncvarid, &
     877&                                akm(1:nwzmax))
     878
     879            ncret = nf90_def_var(ncid, 'bkm', NF90_FLOAT, &
     880&                                       dim1dids, ncvarid)
     881            ncret = nf90_def_var_deflate(ncid, ncvarid,   &
     882&                                        shuffle=0,     &
     883&                                        deflate=1,     &
     884&                                        deflate_level=DEF_LEVEL)
     885            ncret = nf90_put_var(ncid, ncvarid, &
     886&                                bkm(1:nwzmax))
     887
     888
     889            dim1dids = (/nuvzmax/)
     890
     891            ncret = nf90_def_var(ncid, 'akz', NF90_FLOAT, &
     892&                                       dim1dids, ncvarid)
     893            ncret = nf90_def_var_deflate(ncid, ncvarid,   &
     894&                                        shuffle=0,     &
     895&                                        deflate=1,     &
     896&                                        deflate_level=DEF_LEVEL)
     897            ncret = nf90_put_var(ncid, ncvarid, &
     898&                                akz(1:nuvzmax))
     899
     900            ncret = nf90_def_var(ncid, 'bkz', NF90_FLOAT, &
     901&                                       dim1dids, ncvarid)
     902            ncret = nf90_def_var_deflate(ncid, ncvarid,   &
     903&                                        shuffle=0,     &
     904&                                        deflate=1,     &
     905&                                        deflate_level=DEF_LEVEL)
     906            ncret = nf90_put_var(ncid, ncvarid, &
     907&                                bkz(1:nuvzmax))
     908
     909
     910            dim1dids = (/nzmax/)
     911
     912            ncret = nf90_def_var(ncid, 'aknew', NF90_FLOAT, &
     913&                                       dim1dids, ncvarid)
     914            ncret = nf90_def_var_deflate(ncid, ncvarid,   &
     915&                                        shuffle=0,     &
     916&                                        deflate=1,     &
     917&                                        deflate_level=DEF_LEVEL)
     918            ncret = nf90_put_var(ncid, ncvarid, &
     919&                                aknew(1:nzmax))
     920
     921            ncret = nf90_def_var(ncid, 'bknew', NF90_FLOAT, &
     922&                                       dim1dids, ncvarid)
     923            ncret = nf90_def_var_deflate(ncid, ncvarid,   &
     924&                                        shuffle=0,     &
     925&                                        deflate=1,     &
     926&                                        deflate_level=DEF_LEVEL)
     927            ncret = nf90_put_var(ncid, ncvarid, &
     928&                                bknew(1:nzmax))
     929
     930
     931            PRINT *, 'SUM(bknew(1:nzmax)): ', &
     932&                                        SUM(bknew(1:nzmax))
     933
     934
     935
     936
     937
    843938
    844939            ! Nested, scalar values (for each nest)
     
    11561251
    11571252
    1158 
    11591253            ! 2d fields
    11601254            READ(iounit) ps(:,:,:,cm_index)
     
    12471341
    12481342
     1343            ncret = nf90_inq_varid(ncid, 'vdep', ncvarid)
     1344            ncret = nf90_get_var(ncid, ncvarid, vdep(0:nxmax-1,0:nymax-1,1:maxspec, cm_index))
     1345
    12491346
    12501347
     
    12591356            READ(iounit) aknew(:)
    12601357            READ(iounit) bknew(:)
     1358
     1359
     1360            ncret = nf90_inq_varid(ncid, 'z0', ncvarid)
     1361            ncret = nf90_get_var(ncid, ncvarid, z0(1:numclass))
     1362
     1363            ncret = nf90_inq_varid(ncid, 'akm', ncvarid)
     1364            ncret = nf90_get_var(ncid, ncvarid, akm(1:nwzmax))
     1365
     1366            ncret = nf90_inq_varid(ncid, 'bkm', ncvarid)
     1367            ncret = nf90_get_var(ncid, ncvarid, bkm(1:nwzmax))
     1368
     1369            ncret = nf90_inq_varid(ncid, 'akz', ncvarid)
     1370            ncret = nf90_get_var(ncid, ncvarid, akz(1:nuvzmax))
     1371
     1372            ncret = nf90_inq_varid(ncid, 'bkz', ncvarid)
     1373            ncret = nf90_get_var(ncid, ncvarid, bkz(1:nuvzmax))
     1374
     1375            ncret = nf90_inq_varid(ncid, 'aknew', ncvarid)
     1376            ncret = nf90_get_var(ncid, ncvarid, aknew(1:nzmax))
     1377
     1378            ncret = nf90_inq_varid(ncid, 'bknew', ncvarid)
     1379            ncret = nf90_get_var(ncid, ncvarid, bknew(1:nzmax))
     1380
     1381
     1382
     1383            PRINT *, 'SUM(bknew(1:nzmax)): ', &
     1384&                                        SUM(bknew(1:nzmax))
     1385
     1386
     1387
    12611388
    12621389
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG