Changeset 5270 for LMDZ6/trunk/libf/dyn3dmem
- Timestamp:
- Oct 24, 2024, 1:55:38 PM (7 weeks ago)
- File:
-
- 1 edited
-
LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.f90 (modified) (31 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.f90
r5268 r5270 11 11 USE getparam, only: ini_getparam, fin_getparam, getpar 12 12 USE Write_Field_loc 13 use netcdf, only: nf90_nowrite, nf90_open, nf90_inq_varid, nf90_close, & 14 nf90_inq_dimid, nf90_inquire_dimension 13 USE netcdf, ONLY: nf90_nowrite, nf90_open, nf90_inq_varid, nf90_close, & 14 nf90_inq_dimid, nf90_inquire_dimension, nf90_inq_dimid, & 15 nf90_inquire_dimension, nf90_enddef, nf90_def_dim, nf90_put_var, nf90_noerr, nf90_close, nf90_inq_varid, & 16 nf90_redef, nf90_write, nf90_unlimited, nf90_float, nf90_clobber, nf90_64bit_offset, nf90_float, & 17 nf90_create, nf90_def_var, nf90_open 15 18 USE parallel_lmdz 16 19 USE pres2lev_mod, only: pres2lev … … 81 84 INCLUDE "dimensions.h" 82 85 INCLUDE "paramet.h" 83 INCLUDE "netcdf.inc"84 86 85 87 INTEGER :: error,ncidpl,rid,rcod … … 177 179 if (ncidpl.eq.-99) then 178 180 rcod=nf90_open('apbp.nc',Nf90_NOWRITe, ncidpl) 179 if (rcod.NE. NF_NOERR) THEN181 if (rcod.NE.nf90_noerr) THEN 180 182 abort_message=' Nudging error -> no file apbp.nc' 181 183 CALL abort_gcm(modname,abort_message,1) … … 185 187 if (ncidpl.EQ.-99) then 186 188 rcod=nf90_open('P.nc',Nf90_NOWRITe,ncidpl) 187 if (rcod.NE. NF_NOERR) THEN189 if (rcod.NE.nf90_noerr) THEN 188 190 abort_message=' Nudging error -> no file P.nc' 189 191 CALL abort_gcm(modname,abort_message,1) … … 194 196 if (ncidpl.eq.-99) then 195 197 rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl) 196 if (rcod.NE. NF_NOERR) THEN198 if (rcod.NE.nf90_noerr) THEN 197 199 abort_message=' Nudging error -> no file u.nc' 198 200 CALL abort_gcm(modname,abort_message,1) … … 205 207 if (ncidpl.eq.-99) then 206 208 rcod=nf90_open('v.nc',nf90_nowrite,ncidpl) 207 if (rcod.NE. NF_NOERR) THEN209 if (rcod.NE.nf90_noerr) THEN 208 210 abort_message=' Nudging error -> no file v.nc' 209 211 CALL abort_gcm(modname,abort_message,1) … … 215 217 if (ncidpl.eq.-99) then 216 218 rcod=nf90_open('T.nc',nf90_nowrite,ncidpl) 217 if (rcod.NE. NF_NOERR) THEN219 if (rcod.NE.nf90_noerr) THEN 218 220 abort_message=' Nudging error -> no file T.nc' 219 221 CALL abort_gcm(modname,abort_message,1) … … 226 228 if (ncidpl.eq.-99) then 227 229 rcod=nf90_open('hur.nc',nf90_nowrite, ncidpl) 228 if (rcod.NE. NF_NOERR) THEN230 if (rcod.NE.nf90_noerr) THEN 229 231 abort_message=' Nudging error -> no file hur.nc' 230 232 CALL abort_gcm(modname,abort_message,1) … … 234 236 235 237 endif 236 error= NF_INQ_DIMID(ncidpl,'LEVEL',rid)237 IF (error.NE. NF_NOERR) error=NF_INQ_DIMID(ncidpl,'PRESSURE',rid)238 IF (error.NE. NF_NOERR) THEN238 error=nf90_inq_dimid(ncidpl,'LEVEL',rid) 239 IF (error.NE.nf90_noerr) error=nf90_inq_dimid(ncidpl,'PRESSURE',rid) 240 IF (error.NE.nf90_noerr) THEN 239 241 abort_message='Nudging: error reading pressure levels' 240 242 CALL abort_gcm(modname,abort_message,1) 241 243 ENDIF 242 error= NF_INQ_DIMLEN(ncidpl,rid,nlevnc)244 error=nf90_inquire_dimension(ncidpl,rid,len=nlevnc) 243 245 write(*,*)trim(modname)//' : number of vertical levels nlevnc', nlevnc 244 246 rcod = nf90_close(ncidpl) … … 1575 1577 USE netcdf, ONLY: nf90_put_var 1576 1578 IMPLICIT NONE 1577 1578 include "netcdf.inc"1579 1579 include "dimensions.h" 1580 1580 include "paramet.h" … … 1605 1605 write(*,*),trim(modname)//' Reading nudging on model levels' 1606 1606 rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl) 1607 IF (rcode.NE. NF_NOERR) THEN1607 IF (rcode.NE.nf90_noerr) THEN 1608 1608 abort_message='Nudging: error -> no file apbp.nc' 1609 1609 CALL abort_gcm(modname,abort_message,1) 1610 1610 ENDIF 1611 1611 rcode = nf90_inq_varid(ncidpl, 'AP', varidap) 1612 IF (rcode.NE. NF_NOERR) THEN1612 IF (rcode.NE.nf90_noerr) THEN 1613 1613 abort_message='Nudging: error -> no AP variable in file apbp.nc' 1614 1614 CALL abort_gcm(modname,abort_message,1) 1615 1615 ENDIF 1616 1616 rcode = nf90_inq_varid(ncidpl, 'BP', varidbp) 1617 IF (rcode.NE. NF_NOERR) THEN1617 IF (rcode.NE.nf90_noerr) THEN 1618 1618 abort_message='Nudging: error -> no BP variable in file apbp.nc' 1619 1619 CALL abort_gcm(modname,abort_message,1) … … 1625 1625 if (guide_plevs.EQ.2) then 1626 1626 rcode = nf90_open('P.nc', nf90_nowrite, ncidp) 1627 IF (rcode.NE. NF_NOERR) THEN1627 IF (rcode.NE.nf90_noerr) THEN 1628 1628 abort_message='Nudging: error -> no file P.nc' 1629 1629 CALL abort_gcm(modname,abort_message,1) 1630 1630 ENDIF 1631 1631 rcode = nf90_inq_varid(ncidp, 'PRES', varidp) 1632 IF (rcode.NE. NF_NOERR) THEN1632 IF (rcode.NE.nf90_noerr) THEN 1633 1633 abort_message='Nudging: error -> no PRES variable in file P.nc' 1634 1634 CALL abort_gcm(modname,abort_message,1) … … 1641 1641 if (guide_u) then 1642 1642 rcode = nf90_open('u.nc', nf90_nowrite, ncidu) 1643 IF (rcode.NE. NF_NOERR) THEN1643 IF (rcode.NE.nf90_noerr) THEN 1644 1644 abort_message='Nudging: error -> no file u.nc' 1645 1645 CALL abort_gcm(modname,abort_message,1) 1646 1646 ENDIF 1647 1647 rcode = nf90_inq_varid(ncidu, 'UWND', varidu) 1648 IF (rcode.NE. NF_NOERR) THEN1648 IF (rcode.NE.nf90_noerr) THEN 1649 1649 abort_message='Nudging: error -> no UWND variable in file u.nc' 1650 1650 CALL abort_gcm(modname,abort_message,1) … … 1673 1673 if (guide_v) then 1674 1674 rcode = nf90_open('v.nc', nf90_nowrite, ncidv) 1675 IF (rcode.NE. NF_NOERR) THEN1675 IF (rcode.NE.nf90_noerr) THEN 1676 1676 abort_message='Nudging: error -> no file v.nc' 1677 1677 CALL abort_gcm(modname,abort_message,1) 1678 1678 ENDIF 1679 1679 rcode = nf90_inq_varid(ncidv, 'VWND', varidv) 1680 IF (rcode.NE. NF_NOERR) THEN1680 IF (rcode.NE.nf90_noerr) THEN 1681 1681 abort_message='Nudging: error -> no VWND variable in file v.nc' 1682 1682 CALL abort_gcm(modname,abort_message,1) … … 1706 1706 if (guide_T) then 1707 1707 rcode = nf90_open('T.nc', nf90_nowrite, ncidt) 1708 IF (rcode.NE. NF_NOERR) THEN1708 IF (rcode.NE.nf90_noerr) THEN 1709 1709 abort_message='Nudging: error -> no file T.nc' 1710 1710 CALL abort_gcm(modname,abort_message,1) 1711 1711 ENDIF 1712 1712 rcode = nf90_inq_varid(ncidt, 'AIR', varidt) 1713 IF (rcode.NE. NF_NOERR) THEN1713 IF (rcode.NE.nf90_noerr) THEN 1714 1714 abort_message='Nudging: error -> no AIR variable in file T.nc' 1715 1715 CALL abort_gcm(modname,abort_message,1) … … 1737 1737 if (guide_Q) then 1738 1738 rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ) 1739 IF (rcode.NE. NF_NOERR) THEN1739 IF (rcode.NE.nf90_noerr) THEN 1740 1740 abort_message='Nudging: error -> no file hur.nc' 1741 1741 CALL abort_gcm(modname,abort_message,1) 1742 1742 ENDIF 1743 1743 rcode = nf90_inq_varid(ncidQ, 'RH', varidQ) 1744 IF (rcode.NE. NF_NOERR) THEN1744 IF (rcode.NE.nf90_noerr) THEN 1745 1745 abort_message='Nudging: error -> no RH variable in file hur.nc' 1746 1746 CALL abort_gcm(modname,abort_message,1) … … 1769 1769 if ((guide_P).OR.(guide_plevs.EQ.1)) then 1770 1770 rcode = nf90_open('ps.nc', nf90_nowrite, ncidps) 1771 IF (rcode.NE. NF_NOERR) THEN1771 IF (rcode.NE.nf90_noerr) THEN 1772 1772 abort_message='Nudging: error -> no file ps.nc' 1773 1773 CALL abort_gcm(modname,abort_message,1) 1774 1774 ENDIF 1775 1775 rcode = nf90_inq_varid(ncidps, 'SP', varidps) 1776 IF (rcode.NE. NF_NOERR) THEN1776 IF (rcode.NE.nf90_noerr) THEN 1777 1777 abort_message='Nudging: error -> no SP variable in file ps.nc' 1778 1778 CALL abort_gcm(modname,abort_message,1) … … 1894 1894 USE netcdf, ONLY: nf90_put_var 1895 1895 IMPLICIT NONE 1896 1897 include "netcdf.inc"1898 1896 include "dimensions.h" 1899 1897 include "paramet.h" … … 1927 1925 write(*,*)trim(modname)//' Reading nudging on model levels' 1928 1926 rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl) 1929 IF (rcode.NE. NF_NOERR) THEN1927 IF (rcode.NE.nf90_noerr) THEN 1930 1928 abort_message='Nudging: error -> no file apbp.nc' 1931 1929 CALL abort_gcm(modname,abort_message,1) 1932 1930 ENDIF 1933 1931 rcode = nf90_inq_varid(ncidpl, 'AP', varidap) 1934 IF (rcode.NE. NF_NOERR) THEN1932 IF (rcode.NE.nf90_noerr) THEN 1935 1933 abort_message='Nudging: error -> no AP variable in file apbp.nc' 1936 1934 CALL abort_gcm(modname,abort_message,1) 1937 1935 ENDIF 1938 1936 rcode = nf90_inq_varid(ncidpl, 'BP', varidbp) 1939 IF (rcode.NE. NF_NOERR) THEN1937 IF (rcode.NE.nf90_noerr) THEN 1940 1938 abort_message='Nudging: error -> no BP variable in file apbp.nc' 1941 1939 CALL abort_gcm(modname,abort_message,1) … … 1946 1944 if (guide_plevs.EQ.2) then 1947 1945 rcode = nf90_open('P.nc', nf90_nowrite, ncidp) 1948 IF (rcode.NE. NF_NOERR) THEN1946 IF (rcode.NE.nf90_noerr) THEN 1949 1947 abort_message='Nudging: error -> no file P.nc' 1950 1948 CALL abort_gcm(modname,abort_message,1) 1951 1949 ENDIF 1952 1950 rcode = nf90_inq_varid(ncidp, 'PRES', varidp) 1953 IF (rcode.NE. NF_NOERR) THEN1951 IF (rcode.NE.nf90_noerr) THEN 1954 1952 abort_message='Nudging: error -> no PRES variable in file P.nc' 1955 1953 CALL abort_gcm(modname,abort_message,1) … … 1961 1959 if (guide_u) then 1962 1960 rcode = nf90_open('u.nc', nf90_nowrite, ncidu) 1963 IF (rcode.NE. NF_NOERR) THEN1961 IF (rcode.NE.nf90_noerr) THEN 1964 1962 abort_message='Nudging: error -> no file u.nc' 1965 1963 CALL abort_gcm(modname,abort_message,1) 1966 1964 ENDIF 1967 1965 rcode = nf90_inq_varid(ncidu, 'UWND', varidu) 1968 IF (rcode.NE. NF_NOERR) THEN1966 IF (rcode.NE.nf90_noerr) THEN 1969 1967 abort_message='Nudging: error -> no UWND variable in file u.nc' 1970 1968 CALL abort_gcm(modname,abort_message,1) … … 1977 1975 if (guide_v) then 1978 1976 rcode = nf90_open('v.nc', nf90_nowrite, ncidv) 1979 IF (rcode.NE. NF_NOERR) THEN1977 IF (rcode.NE.nf90_noerr) THEN 1980 1978 abort_message='Nudging: error -> no file v.nc' 1981 1979 CALL abort_gcm(modname,abort_message,1) 1982 1980 ENDIF 1983 1981 rcode = nf90_inq_varid(ncidv, 'VWND', varidv) 1984 IF (rcode.NE. NF_NOERR) THEN1982 IF (rcode.NE.nf90_noerr) THEN 1985 1983 abort_message='Nudging: error -> no VWND variable in file v.nc' 1986 1984 CALL abort_gcm(modname,abort_message,1) … … 1992 1990 if (guide_T) then 1993 1991 rcode = nf90_open('T.nc', nf90_nowrite, ncidt) 1994 IF (rcode.NE. NF_NOERR) THEN1992 IF (rcode.NE.nf90_noerr) THEN 1995 1993 abort_message='Nudging: error -> no file T.nc' 1996 1994 CALL abort_gcm(modname,abort_message,1) 1997 1995 ENDIF 1998 1996 rcode = nf90_inq_varid(ncidt, 'AIR', varidt) 1999 IF (rcode.NE. NF_NOERR) THEN1997 IF (rcode.NE.nf90_noerr) THEN 2000 1998 abort_message='Nudging: error -> no AIR variable in file T.nc' 2001 1999 CALL abort_gcm(modname,abort_message,1) … … 2007 2005 if (guide_Q) then 2008 2006 rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ) 2009 IF (rcode.NE. NF_NOERR) THEN2007 IF (rcode.NE.nf90_noerr) THEN 2010 2008 abort_message='Nudging: error -> no file hur.nc' 2011 2009 CALL abort_gcm(modname,abort_message,1) 2012 2010 ENDIF 2013 2011 rcode = nf90_inq_varid(ncidQ, 'RH', varidQ) 2014 IF (rcode.NE. NF_NOERR) THEN2012 IF (rcode.NE.nf90_noerr) THEN 2015 2013 abort_message='Nudging: error -> no RH,variable in file hur.nc' 2016 2014 CALL abort_gcm(modname,abort_message,1) … … 2022 2020 if ((guide_P).OR.(guide_plevs.EQ.1)) then 2023 2021 rcode = nf90_open('ps.nc', nf90_nowrite, ncidps) 2024 IF (rcode.NE. NF_NOERR) THEN2022 IF (rcode.NE.nf90_noerr) THEN 2025 2023 abort_message='Nudging: error -> no file ps.nc' 2026 2024 CALL abort_gcm(modname,abort_message,1) 2027 2025 ENDIF 2028 2026 rcode = nf90_inq_varid(ncidps, 'SP', varidps) 2029 IF (rcode.NE. NF_NOERR) THEN2027 IF (rcode.NE.nf90_noerr) THEN 2030 2028 abort_message='Nudging: error -> no SP variable in file ps.nc' 2031 2029 CALL abort_gcm(modname,abort_message,1) … … 2177 2175 INCLUDE "dimensions.h" 2178 2176 INCLUDE "paramet.h" 2179 INCLUDE "netcdf.inc"2180 2177 INCLUDE "comgeom2.h" 2181 2178 … … 2234 2231 ! ---------------------------------------------- 2235 2232 ! Ouverture du fichier 2236 ierr= NF_CREATE("guide_ins.nc",IOR(NF_CLOBBER,NF_64BIT_OFFSET),nid)2233 ierr=nf90_create("guide_ins.nc",IOR(nf90_clobber,nf90_64bit_offset),nid) 2237 2234 ! Definition des dimensions 2238 ierr= NF_DEF_DIM(nid,"LONU",iip1,id_lonu)2239 ierr= NF_DEF_DIM(nid,"LONV",iip1,id_lonv)2240 ierr= NF_DEF_DIM(nid,"LATU",jjp1,id_latu)2241 ierr= NF_DEF_DIM(nid,"LATV",jjm,id_latv)2242 ierr= NF_DEF_DIM(nid,"LEVEL",llm,id_lev)2243 ierr= NF_DEF_DIM(nid,"TIME",NF_UNLIMITED,id_tim)2235 ierr=nf90_def_dim(nid,"LONU",iip1,id_lonu) 2236 ierr=nf90_def_dim(nid,"LONV",iip1,id_lonv) 2237 ierr=nf90_def_dim(nid,"LATU",jjp1,id_latu) 2238 ierr=nf90_def_dim(nid,"LATV",jjm,id_latv) 2239 ierr=nf90_def_dim(nid,"LEVEL",llm,id_lev) 2240 ierr=nf90_def_dim(nid,"TIME",nf90_unlimited,id_tim) 2244 2241 2245 2242 ! Creation des variables dimensions 2246 ierr= NF_DEF_VAR(nid,"LONU",NF_FLOAT,1,id_lonu,vid_lonu)2247 ierr= NF_DEF_VAR(nid,"LONV",NF_FLOAT,1,id_lonv,vid_lonv)2248 ierr= NF_DEF_VAR(nid,"LATU",NF_FLOAT,1,id_latu,vid_latu)2249 ierr= NF_DEF_VAR(nid,"LATV",NF_FLOAT,1,id_latv,vid_latv)2250 ierr= NF_DEF_VAR(nid,"LEVEL",NF_FLOAT,1,id_lev,vid_lev)2251 ierr= NF_DEF_VAR(nid,"cu",NF_FLOAT,2,(/id_lonu,id_latu/),vid_cu)2252 ierr= NF_DEF_VAR(nid,"cv",NF_FLOAT,2,(/id_lonv,id_latv/),vid_cv)2253 ierr= NF_DEF_VAR(nid,"au",NF_FLOAT,2,(/id_lonu,id_latu/),vid_au)2254 ierr= NF_DEF_VAR(nid,"av",NF_FLOAT,2,(/id_lonv,id_latv/),vid_av)2243 ierr=nf90_def_var(nid,"LONU",nf90_float,id_lonu,vid_lonu) 2244 ierr=nf90_def_var(nid,"LONV",nf90_float,id_lonv,vid_lonv) 2245 ierr=nf90_def_var(nid,"LATU",nf90_float,id_latu,vid_latu) 2246 ierr=nf90_def_var(nid,"LATV",nf90_float,id_latv,vid_latv) 2247 ierr=nf90_def_var(nid,"LEVEL",nf90_float,id_lev,vid_lev) 2248 ierr=nf90_def_var(nid,"cu",nf90_float,(/id_lonu,id_latu/),vid_cu) 2249 ierr=nf90_def_var(nid,"cv",nf90_float,(/id_lonv,id_latv/),vid_cv) 2250 ierr=nf90_def_var(nid,"au",nf90_float,(/id_lonu,id_latu/),vid_au) 2251 ierr=nf90_def_var(nid,"av",nf90_float,(/id_lonv,id_latv/),vid_av) 2255 2252 call nf95_def_var(nid, "alpha_T", nf90_float, (/id_lonv, id_latu/), & 2256 2253 varid_alpha_t) … … 2258 2255 varid_alpha_q) 2259 2256 2260 ierr= NF_ENDDEF(nid)2257 ierr=nf90_enddef(nid) 2261 2258 2262 2259 ! Enregistrement des variables dimensions … … 2275 2272 ! Cr�ation des variables sauvegard�es 2276 2273 ! -------------------------------------------------------------------- 2277 ierr = NF_REDEF(nid)2274 ierr = nf90_redef(nid) 2278 2275 ! Pressure (GCM) 2279 2276 dim4=(/id_lonv,id_latu,id_lev,id_tim/) 2280 ierr = NF_DEF_VAR(nid,"SP",NF_FLOAT,4,dim4,varid)2277 ierr = nf90_def_var(nid,"SP",nf90_float,dim4,varid) 2281 2278 ! Surface pressure (guidage) 2282 2279 IF (guide_P) THEN 2283 2280 dim3=(/id_lonv,id_latu,id_tim/) 2284 ierr = NF_DEF_VAR(nid,"ps",NF_FLOAT,3,dim3,varid)2281 ierr = nf90_def_var(nid,"ps",nf90_float,dim3,varid) 2285 2282 ENDIF 2286 2283 ! Zonal wind 2287 2284 IF (guide_u) THEN 2288 2285 dim4=(/id_lonu,id_latu,id_lev,id_tim/) 2289 ierr = NF_DEF_VAR(nid,"u",NF_FLOAT,4,dim4,varid)2290 ierr = NF_DEF_VAR(nid,"ua",NF_FLOAT,4,dim4,varid)2291 ierr = NF_DEF_VAR(nid,"ucov",NF_FLOAT,4,dim4,varid)2286 ierr = nf90_def_var(nid,"u",nf90_float,dim4,varid) 2287 ierr = nf90_def_var(nid,"ua",nf90_float,dim4,varid) 2288 ierr = nf90_def_var(nid,"ucov",nf90_float,dim4,varid) 2292 2289 ENDIF 2293 2290 ! Merid. wind 2294 2291 IF (guide_v) THEN 2295 2292 dim4=(/id_lonv,id_latv,id_lev,id_tim/) 2296 ierr = NF_DEF_VAR(nid,"v",NF_FLOAT,4,dim4,varid)2297 ierr = NF_DEF_VAR(nid,"va",NF_FLOAT,4,dim4,varid)2298 ierr = NF_DEF_VAR(nid,"vcov",NF_FLOAT,4,dim4,varid)2293 ierr = nf90_def_var(nid,"v",nf90_float,dim4,varid) 2294 ierr = nf90_def_var(nid,"va",nf90_float,dim4,varid) 2295 ierr = nf90_def_var(nid,"vcov",nf90_float,dim4,varid) 2299 2296 ENDIF 2300 2297 ! Pot. Temperature 2301 2298 IF (guide_T) THEN 2302 2299 dim4=(/id_lonv,id_latu,id_lev,id_tim/) 2303 ierr = NF_DEF_VAR(nid,"teta",NF_FLOAT,4,dim4,varid)2300 ierr = nf90_def_var(nid,"teta",nf90_float,dim4,varid) 2304 2301 ENDIF 2305 2302 ! Specific Humidity 2306 2303 IF (guide_Q) THEN 2307 2304 dim4=(/id_lonv,id_latu,id_lev,id_tim/) 2308 ierr = NF_DEF_VAR(nid,"q",NF_FLOAT,4,dim4,varid)2305 ierr = nf90_def_var(nid,"q",nf90_float,dim4,varid) 2309 2306 ENDIF 2310 2307 2311 ierr = NF_ENDDEF(nid)2312 ierr = NF_CLOSE(nid)2308 ierr = nf90_enddef(nid) 2309 ierr = nf90_close(nid) 2313 2310 ENDIF ! timestep=0 2314 2311 … … 2317 2314 ! -------------------------------------------------------------------- 2318 2315 2319 ierr= NF_OPEN("guide_ins.nc",NF_WRITE,nid)2316 ierr=nf90_open("guide_ins.nc",nf90_write,nid) 2320 2317 2321 2318 IF (varname=="SP") timestep=timestep+1 2322 2319 2323 ierr = NF_INQ_VARID(nid,varname,varid)2320 ierr = nf90_inq_varid(nid,varname,varid) 2324 2321 SELECT CASE (varname) 2325 2322 CASE ("SP","ps") … … 2360 2357 2361 2358 ierr = nf90_put_var(nid, varid, field_glo, start, count) 2362 ierr = NF_CLOSE(nid)2359 ierr = nf90_close(nid) 2363 2360 2364 2361 DEALLOCATE(field_glo)
Note: See TracChangeset
for help on using the changeset viewer.
