Changeset 1d9681b in flexpart.git
- Timestamp:
- Feb 17, 2017, 7:38:39 PM (7 years ago)
- Branches:
- fp9.3.1-20161214-nc4
- Children:
- 7e5f1dc
- Parents:
- d2d8a96
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
flexpart_code/fpmetbinary_mod.F90
rd2d8a96 r1d9681b 44 44 USE com_mod 45 45 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 47 47 48 48 USE netcdf … … 152 152 ! Zeroes out, in local datastructures, the values dumped/loaded 153 153 ! 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. 154 156 ! The lines here correspond to READ and WRITE in the dump/load routines 155 157 … … 830 832 831 833 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)) 832 844 833 845 … … 841 853 WRITE(iounit) aknew(:) 842 854 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 843 938 844 939 ! Nested, scalar values (for each nest) … … 1156 1251 1157 1252 1158 1159 1253 ! 2d fields 1160 1254 READ(iounit) ps(:,:,:,cm_index) … … 1247 1341 1248 1342 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 1249 1346 1250 1347 … … 1259 1356 READ(iounit) aknew(:) 1260 1357 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 1261 1388 1262 1389
Note: See TracChangeset
for help on using the changeset viewer.