/[MITgcm]/MITgcm/tools/calc_diagnostics_dims
ViewVC logotype

Contents of /MITgcm/tools/calc_diagnostics_dims

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


Revision 1.1 - (show annotations) (download)
Fri Feb 13 06:05:49 2004 UTC (20 years, 1 month 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, checkpoint52l_pre, hrcube4, checkpoint58e_post, mitgcm_mapl_00, checkpoint52n_post, checkpoint52j_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, checkpoint52l_post, checkpoint55h_post, checkpoint58n_post, checkpoint53b_post, checkpoint58x_post, checkpoint52k_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, checkpoint57v_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, hrcube5, 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, checkpoint52j_pre, 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
 o initial check-in of a script that calculates the exact number of "z"
   levels needed within the diagnostics package
   - script was tested with a data.diagnostics "dummy" file
 o minor tweak to mnc

1 #! /usr/bin/env bash
2 #
3 # $Header: $
4 # $Name: $
5
6 # The purpose of this script is to calculate the exact number of "z"
7 # dimensions needed within a FORTRAN storage array for the MITgcm
8 # diagnostics package.
9
10 usage()
11 {
12 cat <<EOF
13
14 Usage: $0 data_file mitgcm_root var_name [ separator ]
15
16 where:
17 data_file : is the path and file name for the
18 "data.diagnostics" file
19 mitgcm_root : is the path to the "root" directory of
20 the MITgcm source tree
21 var_name : is the name of the variable(s) containing
22 the diagnostic strings
23 separator : an optional argument specifying the
24 separator character within the data
25 file (default='&')
26
27 EOF
28 }
29
30 COMMANDL="$0 $@"
31
32 DATA_FILE="$1"
33 MITGCM_ROOT="$2"
34 VAR_NAME="$3"
35 SEPARATOR="$4"
36 nd_tot=0
37
38 # Check that the arguments were specified and are read-able
39 if test "x$DATA_FILE" = x ; then
40 echo "ERROR: the \"data.diagnostics\" file was not specified"
41 echo " -- please set it using the first argument."
42 usage
43 exit 1
44 fi
45 pack_h="$MITGCM_ROOT"/pkg/diagnostics/diagnostics.h
46 if test ! -r $pack_h ; then
47 echo "ERROR: cannot read file \"$pack_h\" "
48 echo " -- please check that the file exists and that "
49 echo " \$MITGCM_ROOT is correctly set using the "
50 echo " second argument."
51 usage
52 exit 1
53 fi
54 init_vals="$MITGCM_ROOT"/pkg/diagnostics/diagnostics_init_vals.F
55 if test ! -r $init_vals ; then
56 echo "ERROR: cannot read file \"$init_vals\" "
57 echo " -- please check that the file exists and that "
58 echo " \$MITGCM_ROOT is correctly set using the "
59 echo " second argument."
60 usage
61 exit 1
62 fi
63 fizhi_SIZE="$MITGCM_ROOT"/pkg/fizhi/fizhi_SIZE.h
64 if test ! -r $fizhi_SIZE ; then
65 echo "ERROR: cannot read file \"$fizhi_SIZE\" "
66 echo " -- please check that the file exists and that "
67 echo " \$MITGCM_ROOT is correctly set using the "
68 echo " second argument."
69 usage
70 exit 1
71 fi
72 if test "x$VAR_NAME" = x ; then
73 echo "ERROR: \"\$VAR_NAME\" was not specified -- please set it"
74 echo " using the third argument."
75 usage
76 exit 1
77 fi
78 if test "x$SEPARATOR" = x ; then
79 SEPARATOR='&'
80 fi
81 if test ! -r $DATA_FILE ; then
82 echo "ERROR: cannot read file \"$DATA_FILE\""
83 usage
84 exit 1
85 fi
86
87
88 # Get the diagnostic names
89 echo -n "" > ./tmp_diagnostic_names
90 DNAMES=
91 vcode=0
92 cat $DATA_FILE | while read line ; do
93 r0=t
94 r1=t
95 echo $line | grep '^[ ]*#' > /dev/null 2>&1 && r0=f
96 echo $line | grep "$SEPARATOR" > /dev/null 2>&1 && r1=f
97 if test "x$vcode" = x1 ; then
98 echo $line | grep '=' > /dev/null 2>&1 && vcode=0
99 fi
100 echo $line | grep "$VAR_NAME"'[ ]*=' > /dev/null 2>&1 && vcode=1
101 if test $r0 = t -a $r1 = t -a ! "x$vcode" = x0 ; then
102 t1=`echo $line | sed -e "s|$VAR_NAME| |g" | sed -e 's|=| |g'`
103 t2=`echo $t1 | sed -e "s|'| |g" | sed -e 's|,| |g'`
104 echo "$t2" >> ./tmp_diagnostic_names
105 fi
106 done
107 DNAMES=`cat ./tmp_diagnostic_names`
108 rm -rf ./tmp_diagnostic_names
109
110 # Get the size of $NRPHYS
111 t1=`cat $fizhi_SIZE | grep -i Nrphys | grep -i "^[ ]*parameter"`
112 t2=`echo $t1 | sed -e 's|(| |g' | sed -e 's|)| |g' | sed -e 's|=| |g'`
113 NRPHYS=`echo $t2 | awk '{print $3}'`
114
115 # Get the number of "z" dimensions
116 numz=0
117 for dnam in $DNAMES ; do
118 t1=`grep -i "n"$dnam $pack_h | grep -i '^[ ]*EQUIVALENCE'`
119 t2=`echo $t1 | sed -e 's|(| |g' | sed -e 's|)| |g'`
120 t3=`echo $t2 | awk '{print $3}'`
121 if test ! "x$t3" = x ; then
122 t1=`cat $init_vals | grep -i '^[ ]*KDIAG' | grep $t3`
123 t2=`echo $t1 | sed -e 's|(| |g' | sed -e 's|)| |g'`
124 t3=`echo $t2 | sed -e 's|=| |g' | awk '{print $3}'`
125 t1=`echo $t3 | sed -e "s|nrphys|$NRPHYS|g"`
126 numz=$(( $numz + $t1 ))
127 fi
128 done
129 echo "$numz"

  ViewVC Help
Powered by ViewVC 1.1.22