/[MITgcm]/MITgcm/pkg/mnc/mnc_test_001.T
ViewVC logotype

Diff of /MITgcm/pkg/mnc/mnc_test_001.T

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by edhill, Wed Jan 7 19:50:52 2004 UTC revision 1.5 by edhill, Thu Jan 8 07:24:47 2004 UTC
# Line 21  C     Local Variables Line 21  C     Local Variables
21        integer ndim, dims(MAX_DIM)        integer ndim, dims(MAX_DIM)
22        integer myThid, i,j,k,l,ii, n, ind, fid        integer myThid, i,j,k,l,ii, n, ind, fid
23        character*(MAX_LEN_MBUF) msgbuf        character*(MAX_LEN_MBUF) msgbuf
24        character*(MAX_STR) blank, fname, gname, vname, units        character*(MAX_STR) blank, fname, gname1,gname2, vname1,vname2
25        character*(MAX_STR) atname, atval        character*(MAX_STR) atname, atval, units
26        character*(MAX_STR) dnames(MAX_NAMES)        character*(MAX_STR) dnames(MAX_NAMES)
27        REAL*8 fillval        REAL*8 fillval
28        REAL*8 var(MAX_STOR)        REAL*8 var(MAX_STOR), tvar(MAX_STOR)
29    
30        print*, '=== Starting Program:  mnc_test_001 ==='        print*, '=== Starting Program:  mnc_test_001 ==='
31    
# Line 37  C     Local Variables Line 37  C     Local Variables
37          dnames(i)(1:MAX_STR) = blank(1:MAX_STR)          dnames(i)(1:MAX_STR) = blank(1:MAX_STR)
38        ENDDO        ENDDO
39        fname(1:MAX_STR) = blank(1:MAX_STR)        fname(1:MAX_STR) = blank(1:MAX_STR)
40        gname(1:MAX_STR) = blank(1:MAX_STR)        gname1(1:MAX_STR) = blank(1:MAX_STR)
41        vname(1:MAX_STR) = blank(1:MAX_STR)        gname2(1:MAX_STR) = blank(1:MAX_STR)
42          vname1(1:MAX_STR) = blank(1:MAX_STR)
43          vname2(1:MAX_STR) = blank(1:MAX_STR)
44        units(1:MAX_STR) = blank(1:MAX_STR)        units(1:MAX_STR) = blank(1:MAX_STR)
45        atname(1:MAX_STR) = blank(1:MAX_STR)        atname(1:MAX_STR) = blank(1:MAX_STR)
46        atval(1:MAX_STR) = blank(1:MAX_STR)        atval(1:MAX_STR) = blank(1:MAX_STR)
# Line 56  C     Local Variables Line 58  C     Local Variables
58  C     In FORTRAN, the first variable cycles most quickly  C     In FORTRAN, the first variable cycles most quickly
59        DO ii = 1,MAX_STOR        DO ii = 1,MAX_STOR
60          var(ii) = 0.21          var(ii) = 0.21
61            tvar(ii) = 0.23
62        ENDDO        ENDDO
63        ii = 0        ii = 0
64        DO l = 1,dims(3)        DO l = 1,dims(3)
# Line 82  C     In FORTRAN, the first variable cyc Line 85  C     In FORTRAN, the first variable cyc
85        CALL MNC_FILE_CREATE(myThid, fname)        CALL MNC_FILE_CREATE(myThid, fname)
86    
87        print*, ' Calling MNC_GRID_INIT() ...'        print*, ' Calling MNC_GRID_INIT() ...'
88        gname(1:8) = 'XY_Grid '        gname1(1:10) = 'xyzt      '
89        CALL MNC_GRID_INIT(myThid, fname, gname, ndim, dnames)        CALL MNC_GRID_INIT(myThid, fname, gname1, ndim, dnames)
90          gname2(1:10) = 't         '
91          CALL MNC_GRID_INIT(myThid, fname, gname2, 1, dnames(4))
92    
93        print*, ' Calling MNC_VAR_INIT_DBL() ...'        print*, ' Calling MNC_VAR_INIT_DBL() ...'
94        vname(1:11) = 'temperature'        vname1(1:4) = 'time'
95        units(1:3) = 'm/s'        units(1:19) = 'days since 1970-1-1'
96        fillval = -9999.0        CALL MNC_VAR_INIT_DBL(myThid, fname, gname2, vname1, units )
97        CALL MNC_VAR_INIT_DBL(myThid, fname, gname, vname, fillval)        vname2(1:4) = 'temp'
98          units(1:19) = 'm/s                '
99          CALL MNC_VAR_INIT_DBL(myThid, fname, gname1, vname2, units )
100    
101        print*, ' Calling MNC_VAR_ADD_ATTR_STR() ...'        print*, ' Calling MNC_VAR_ADD_ATTR_STR() ...'
102        atname(1:5) = 'units'        atname(1:9) = 'long_name'
103        atval(1:4) = 'cm/s'        atval(1:11) = 'hours since '
104        CALL MNC_VAR_ADD_ATTR_STR(myThid, fname, vname, atname, atval)        CALL MNC_VAR_ADD_ATTR_STR(myThid,fname,vname1, atname, atval)
105          atname(1:11) = 'valid_range'
106          tvar(1) = -23.5
107          tvar(2) = 176.75
108          CALL MNC_VAR_ADD_ATTR_DBL(myThid,fname,vname2, atname, 2, tvar)
109    
110        print*, ' Calling MNC_VAR_WRITE_DBL() ...'        print*, ' Calling MNC_VAR_WRITE_DBL() ...'
111        CALL MNC_VAR_WRITE_DBL(myThid, fname, vname, var)        CALL MNC_VAR_WRITE_DBL(myThid, fname, vname1, var)
112        print*, ' Calling MNC_VAR_WRITE_DBL() ...'        print*, ' Calling MNC_VAR_WRITE_DBL() ...'
113        CALL MNC_VAR_WRITE_DBL(myThid, fname, vname, var)        CALL MNC_VAR_WRITE_DBL(myThid, fname, vname2, var)
114    
115    C     CALL MNC_DUMP()
116    
117        print*, ' Calling MNC_FILE_CLOSE() ...'        print*, ' Calling MNC_FILE_CLOSE() ...'
118        CALL MNC_FILE_CLOSE(myThid, fname)        CALL MNC_FILE_CLOSE(myThid, fname)
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22