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

Annotation of /MITgcm/tools/calc_diagnostics_dims

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


Revision 1.1 - (hide annotations) (download)
Fri Feb 13 06:05:49 2004 UTC (16 years, 11 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, 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 edhill 1.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