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

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

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

Revision 1.11 - (show annotations) (download) (as text)
Fri Mar 19 03:28:37 2004 UTC (20 years, 2 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint58l_post, checkpoint64z, checkpoint57t_post, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint57o_post, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint58e_post, checkpoint57v_post, checkpoint52n_post, checkpoint53d_post, checkpoint58u_post, checkpoint58w_post, checkpoint54a_pre, checkpoint57m_post, checkpoint55c_post, checkpoint54e_post, checkpoint57s_post, checkpoint54a_post, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint53c_post, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint64, checkpoint65, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint57b_post, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint57e_post, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint55h_post, checkpoint58n_post, checkpoint58x_post, checkpoint57g_pre, checkpoint54b_post, checkpoint53b_pre, checkpoint55b_post, checkpoint58t_post, checkpoint58h_post, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint54d_post, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint56c_post, checkpoint52m_post, checkpoint57y_pre, checkpoint55, checkpoint53a_post, checkpoint57f_pre, checkpoint57a_post, checkpoint54, checkpoint58q_post, checkpoint54f_post, checkpoint53b_post, checkpoint59q, checkpoint59p, checkpoint55g_post, checkpoint59r, checkpoint58j_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, checkpoint53, eckpoint57e_pre, checkpoint57h_done, checkpoint58f_post, checkpoint53g_post, checkpoint57x_post, checkpoint57n_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58g_post, checkpoint58o_post, checkpoint57z_post, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint57c_post, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint58y_post, checkpoint55e_post, checkpoint58k_post, checkpoint58v_post, checkpoint53f_post, checkpoint55a_post, checkpoint53d_pre, checkpoint54c_post, checkpoint58s_post, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint58p_post, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint57j_post, checkpoint61z, checkpoint61x, checkpoint61y, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint57h_post, checkpoint56a_post, checkpoint55d_post, HEAD
Changes since 1.10: +16 -16 lines
File MIME type: application/x-troff 
 o edit all MNC subroutines so that myThid is the _last_ argument
1 C $Header: /u/u3/gcmpack/MITgcm/pkg/mnc/mnc_test_001.T,v 1.10 2004/02/04 05:45:09 edhill Exp $
2 C $Name: $
3
4 C Simple test program for MNC
5
6 PROGRAM mnc_test_001
7
8 implicit none
9 #include "MNC_OPTIONS.h"
10 #include "EEPARAMS.h"
11
12 C Functions
13 integer ILNBLNK
14
15 C Local Variables
16 integer MAX_STR, MAX_NAMES, MAX_DIM, MAX_STOR
17 PARAMETER ( MAX_STR = 100 )
18 PARAMETER ( MAX_NAMES = 10 )
19 PARAMETER ( MAX_DIM = 10 )
20 PARAMETER ( MAX_STOR = 1000 )
21 integer ndim, dims(MAX_DIM)
22 integer myThid, i,j,k,l,ii, n, ind, fid
23 character*(MAX_LEN_MBUF) msgbuf
24 character*(MAX_STR) blank, fname, gname1,gname2, vname1,vname2
25 character*(MAX_STR) atname, atval, units
26 character*(MAX_STR) dnames(MAX_NAMES)
27 REAL*8 fillval
28 REAL*8 var(MAX_STOR), tvar(MAX_STOR)
29
30 print*, '=== Starting Program: mnc_test_001 ==='
31
32 print*, ' Initalizing variables ...'
33 DO i = 1,MAX_STR
34 blank(i:i) = ' '
35 ENDDO
36 DO i = 1,MAX_NAMES
37 dnames(i)(1:MAX_STR) = blank(1:MAX_STR)
38 ENDDO
39 fname(1:MAX_STR) = blank(1:MAX_STR)
40 gname1(1:MAX_STR) = blank(1:MAX_STR)
41 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)
45 atname(1:MAX_STR) = blank(1:MAX_STR)
46 atval(1:MAX_STR) = blank(1:MAX_STR)
47 fname(1:12) = 'test_001.nc '
48 ndim = 4
49 dims(1) = 7
50 dims(2) = 4
51 dims(3) = 2
52 dims(4) = -1
53 dnames(1)(1:1) = 'X'
54 dnames(2)(1:1) = 'Y'
55 dnames(3)(1:1) = 'Z'
56 dnames(4)(1:1) = 'T'
57 myThid = 1
58 C In FORTRAN, the first variable cycles most quickly
59 DO ii = 1,MAX_STOR
60 var(ii) = 0.21
61 tvar(ii) = 0.23
62 ENDDO
63 ii = 0
64 DO l = 1,dims(3)
65 DO k = 1,dims(3)
66 DO j = 1,dims(2)
67 DO i = 1,dims(1)
68 ii = ii + 1
69 var(ii) = ((i*100 + j)*100 + k)*100 + l
70 ENDDO
71 ENDDO
72 ENDDO
73 ENDDO
74
75 print*, ' Calling MNC_INIT() ...'
76 CALL MNC_INIT(myThid)
77
78 print*, ' Calling MNC_FILE_CREATE() ...'
79 CALL MNC_FILE_CREATE(fname, myThid)
80
81 print*, ' Calling MNC_DIM_INIT() ...'
82 DO i = 1,ndim
83 CALL MNC_DIM_INIT(fname, dnames(i),dims(i), myThid)
84 ii = ILNBLNK(dnames(i))
85 print*, ' ... created dimension: ', dnames(i)(1:ii)
86 ENDDO
87
88 print*, ' Calling MNC_GRID_INIT() ...'
89 gname1(1:10) = 'xyzt '
90 CALL MNC_GRID_INIT(fname, gname1, ndim, dnames, myThid)
91 gname2(1:10) = 't '
92 CALL MNC_GRID_INIT(fname, gname2, 1, dnames(4), myThid)
93
94 print*, ' Calling MNC_VAR_INIT_DBL() ...'
95 vname1(1:4) = 'time'
96 units(1:19) = 'days since 1970-1-1'
97 CALL MNC_VAR_INIT_DBL(fname, gname2, vname1, myThid)
98 vname2(1:4) = 'temp'
99 units(1:19) = 'deg C '
100 CALL MNC_VAR_INIT_DBL(fname, gname1, vname2, myThid)
101
102 print*, ' Calling MNC_VAR_ADD_ATTR_STR() ...'
103 atname(1:9) = 'grid_type'
104 CALL MNC_VAR_ADD_ATTR_INT(fname,vname1, atname, 1, 1, myThid)
105 atname(1:11) = 'valid_range'
106 tvar(1) = -23.5
107 tvar(2) = 176.75
108 CALL MNC_VAR_ADD_ATTR_DBL(fname,vname2, atname, 2, tvar, myThid)
109
110 print*, ' Calling MNC_VAR_WRITE_DBL() ...'
111 CALL MNC_VAR_WRITE_DBL(fname, vname1, var, myThid)
112 print*, ' Calling MNC_VAR_WRITE_DBL() ...'
113 CALL MNC_VAR_WRITE_DBL(fname, vname2, var, myThid)
114
115 CALL MNC_DUMP(myThid)
116
117 print*, ' Calling MNC_FILE_CLOSE() ...'
118 CALL MNC_FILE_CLOSE(fname, myThid)
119
120 CALL MNC_DUMP(myThid)
121
122 print*, '=========================================='
123
124 print*, ' Calling MNC_FILE_READALL() ...'
125 CALL MNC_FILE_READALL(fname, myThid)
126
127 CALL MNC_DUMP(myThid)
128
129 print*, '=== Ending Program: mnc_test_001 ==='
130
131 END
132
133 CEH3 ;;; Local Variables: ***
134 CEH3 ;;; mode:fortran ***
135 CEH3 ;;; End: ***
  ViewVC Help
Powered by ViewVC 1.1.22