/[MITgcm]/MITgcm/pkg/seaice/seaice_mnc_init.F
ViewVC logotype

Annotation of /MITgcm/pkg/seaice/seaice_mnc_init.F

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

Revision 1.4 - (hide annotations) (download)
Mon Mar 20 21:36:11 2006 UTC (18 years, 1 month ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint58e_post, checkpoint58u_post, checkpoint58w_post, checkpoint58r_post, checkpoint58n_post, checkpoint58x_post, checkpoint58t_post, checkpoint58h_post, checkpoint58q_post, checkpoint58j_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59i, checkpoint59h, checkpoint59j, checkpoint59, checkpoint58f_post, checkpoint58d_post, checkpoint58c_post, checkpoint58i_post, checkpoint58g_post, checkpoint58o_post, checkpoint58y_post, checkpoint58k_post, checkpoint58v_post, checkpoint58s_post, checkpoint58p_post, checkpoint58m_post
Changes since 1.3: +21 -1 lines 
  seaice: add an EVP solver following Hunke and Dukowicz (1997) and the
  documentation of CICE. Turn on by defining SEAICE_ALLOW_EVP in
  SEAICE_OPTIONS.h and SEAICEuseEVP=.true. in data.seaice. Works only
  with SEAICE_CGRID defined.
  Use at own risk.
1 mlosch 1.4 C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_mnc_init.F,v 1.3 2005/08/11 02:50:12 edhill Exp $
2 edhill 1.1 C $Name: $
3    
4     #include "SEAICE_OPTIONS.h"
5    
6     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7     CBOP
8     C !ROUTINE: SEAICE_MNC_INIT
9    
10     C !INTERFACE:
11     SUBROUTINE SEAICE_MNC_INIT( myThid )
12    
13     C !DESCRIPTION:
14     C Initialize SEAICE variables within MNC
15    
16     C !USES:
17     IMPLICIT NONE
18     #include "SIZE.h"
19     #include "EEPARAMS.h"
20     #include "PARAMS.h"
21 mlosch 1.4 #ifdef SEAICE_ALLOW_EVP
22     #include "SEAICE_PARAMS.h"
23     #endif /* SEAICE_ALLOW_EVP */
24 edhill 1.1
25     C !INPUT PARAMETERS:
26     INTEGER myThid
27     CEOP
28    
29     #ifdef ALLOW_SEAICE
30     #ifdef ALLOW_MNC
31     IF (useMNC) THEN
32    
33 mlosch 1.4 #ifdef SEAICE_ALLOW_EVP
34     IF ( SEAICEuseEVP ) THEN
35     CALL MNC_CW_ADD_VNAME('si_sigma1','Cen_xy_Hn__-__t',3,4,myThid)
36     CALL MNC_CW_ADD_VATTR_TEXT('si_sigma1','units','N/m^2', myThid)
37     CALL MNC_CW_ADD_VATTR_TEXT('si_sigma1','description',
38     & 'divergence of stress tensor', myThid)
39     CALL MNC_CW_ADD_VNAME('si_sigma2','Cen_xy_Hn__-__t',3,4,myThid)
40     CALL MNC_CW_ADD_VATTR_TEXT('si_sigma2','units','N/m^2', myThid)
41     CALL MNC_CW_ADD_VATTR_TEXT('si_sigma2','description',
42     & 'tension of stress tensor', myThid)
43     CALL MNC_CW_ADD_VNAME('si_sigma12','Cor_xy_Hn__-__t',3,4,myThid)
44     CALL MNC_CW_ADD_VATTR_TEXT('si_sigma12','units','N/m^2', myThid)
45     CALL MNC_CW_ADD_VATTR_TEXT('si_sigma12','description',
46     & 'shear of stress tensor', myThid)
47     ENDIF
48     #endif /* SEAICE_ALLOW_EVP */
49    
50 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_UWIND','U_xy_Hn__-__t',3,4,myThid)
51     CALL MNC_CW_ADD_VATTR_TEXT('si_UWIND','units','m/s', myThid)
52     CALL MNC_CW_ADD_VATTR_TEXT('si_UWIND','description',
53 edhill 1.1 & 'Surface (10m) zonal wind velocity', myThid)
54    
55 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_VWIND','V_xy_Hn__-__t',3,4,myThid)
56     CALL MNC_CW_ADD_VATTR_TEXT('si_VWIND','units','m/s', myThid)
57     CALL MNC_CW_ADD_VATTR_TEXT('si_VWIND','description',
58 edhill 1.1 & 'Surface (10m) meridional wind velocity', myThid)
59    
60 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_EVAP','Cen_xy_Hn__-__t',3,4,myThid)
61     CALL MNC_CW_ADD_VATTR_TEXT('si_EVAP','units','m/s', myThid)
62     CALL MNC_CW_ADD_VATTR_TEXT('si_EVAP','description',
63 edhill 1.1 & 'Evaporation', myThid)
64    
65 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_PRECIP','Cen_xy_Hn__-__t',3,4,myThid)
66     CALL MNC_CW_ADD_VATTR_TEXT('si_PRECIP','units','m/s', myThid)
67     CALL MNC_CW_ADD_VATTR_TEXT('si_PRECIP','description',
68 edhill 1.1 & 'Precipitation', myThid)
69    
70 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_RUNOFF','Cen_xy_Hn__-__t',3,4,myThid)
71     CALL MNC_CW_ADD_VATTR_TEXT('si_RUNOFF','units','m/s', myThid)
72     CALL MNC_CW_ADD_VATTR_TEXT('si_RUNOFF','description',
73 edhill 1.1 & 'River and glacier runoff', myThid)
74    
75 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_ATEMP','Cen_xy_Hn__-__t',3,4,myThid)
76     CALL MNC_CW_ADD_VATTR_TEXT('si_ATEMP','units','deg K', myThid)
77     CALL MNC_CW_ADD_VATTR_TEXT('si_ATEMP','description',
78 edhill 1.1 & 'Surface (2m) air temperature', myThid)
79    
80 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_AQH','Cen_xy_Hn__-__t',3,4,myThid)
81     CALL MNC_CW_ADD_VATTR_TEXT('si_AQH','units','kg/kg', myThid)
82     CALL MNC_CW_ADD_VATTR_TEXT('si_AQH','description',
83 edhill 1.1 & 'Surface (2m) specific humidity', myThid)
84    
85 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_SWDOWN','Cen_xy_Hn__-__t',3,4,myThid)
86     CALL MNC_CW_ADD_VATTR_TEXT('si_SWDOWN','units','W/m^2', myThid)
87     CALL MNC_CW_ADD_VATTR_TEXT('si_SWDOWN','description',
88 edhill 1.1 & 'Downward shortwave radiation', myThid)
89    
90 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_LWDOWN','Cen_xy_Hn__-__t',3,4,myThid)
91     CALL MNC_CW_ADD_VATTR_TEXT('si_LWDOWN','units','W/m^2', myThid)
92     CALL MNC_CW_ADD_VATTR_TEXT('si_LWDOWN','description',
93 edhill 1.1 & 'Downward longwave radiation', myThid)
94    
95 edhill 1.3 CALL MNC_CW_ADD_VNAME('si_UICE','U_xy_Hn__-__t',3,4,myThid)
96     CALL MNC_CW_ADD_VATTR_TEXT('si_UICE','units','m/s', myThid)
97     CALL MNC_CW_ADD_VATTR_TEXT('si_UICE','description',
98     & 'Ice U velocity', myThid)
99    
100     CALL MNC_CW_ADD_VNAME('si_VICE','V_xy_Hn__-__t',3,4,myThid)
101     CALL MNC_CW_ADD_VATTR_TEXT('si_VICE','units','m/s', myThid)
102     CALL MNC_CW_ADD_VATTR_TEXT('si_VICE','description',
103     & 'Ice V velocity', myThid)
104    
105     CALL MNC_CW_ADD_VNAME('si_HEFF','Cen_xy_Hn__-__t',3,4,myThid)
106     CALL MNC_CW_ADD_VATTR_TEXT('si_HEFF','units','m', myThid)
107     CALL MNC_CW_ADD_VATTR_TEXT('si_HEFF','description',
108     & 'Effective ice thickness',myThid)
109    
110     CALL MNC_CW_ADD_VNAME('si_AREA','Cen_xy_Hn__-__t',3,4,myThid)
111     CALL MNC_CW_ADD_VATTR_TEXT('si_AREA','units','m^2/m^2',myThid)
112     CALL MNC_CW_ADD_VATTR_TEXT('si_AREA','description',
113     & 'Fractional ice-covered area',myThid)
114    
115 edhill 1.1 C Time averages
116 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_UICEtave','U_xy_Hn__-__t',3,4,myThid)
117     CALL MNC_CW_ADD_VATTR_TEXT('si_UICEtave','units','m/s', myThid)
118     CALL MNC_CW_ADD_VATTR_TEXT('si_UICEtave','description',
119 edhill 1.1 & 'Time-averaged ', myThid)
120    
121 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_VICEtave','V_xy_Hn__-__t',3,4,myThid)
122     CALL MNC_CW_ADD_VATTR_TEXT('si_VICEtave','units','m/s', myThid)
123     CALL MNC_CW_ADD_VATTR_TEXT('si_VICEtave','description',
124 edhill 1.1 & 'Time-averaged ', myThid)
125    
126 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_FUtave','U_xy_Hn__-__t',3,4,myThid)
127     CALL MNC_CW_ADD_VATTR_TEXT('si_FUtave','units','N/m^2', myThid)
128     CALL MNC_CW_ADD_VATTR_TEXT('si_FUtave','description',
129 edhill 1.1 & 'Time-averaged zonal surface wind stress', myThid)
130    
131 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_FVtave','V_xy_Hn__-__t',3,4,myThid)
132     CALL MNC_CW_ADD_VATTR_TEXT('si_FVtave','units','N/m^2',myThid)
133     CALL MNC_CW_ADD_VATTR_TEXT('si_FVtave','description',
134 edhill 1.1 & 'Time-averaged Meridional surface wind stress', myThid)
135    
136 edhill 1.2 CALL MNC_CW_ADD_VNAME(
137     & 'si_EmPmRtave','Cen_xy_Hn__-__t',3,4,myThid)
138     CALL MNC_CW_ADD_VATTR_TEXT(
139     & 'si_EmPmRtave','units','m/s',myThid)
140     CALL MNC_CW_ADD_VATTR_TEXT('si_EmPmRtave','description',
141 edhill 1.1 & 'Time-averaged net upward freshwater flux', myThid)
142    
143 edhill 1.2 CALL MNC_CW_ADD_VNAME(
144     & 'si_QNETtave','Cen_xy_Hn__-__t',3,4,myThid)
145     CALL MNC_CW_ADD_VATTR_TEXT(
146     & 'si_QNETtave','units','W/m^2', myThid)
147     CALL MNC_CW_ADD_VATTR_TEXT('si_QNETtave','description',
148 edhill 1.1 & 'Time-averaged net upward surface heat flux', myThid)
149    
150 edhill 1.2 CALL MNC_CW_ADD_VNAME(
151     & 'si_QSWtave','Cen_xy_Hn__-__t',3,4,myThid)
152     CALL MNC_CW_ADD_VATTR_TEXT(
153     & 'si_QSWtave','units','W/m^2',myThid)
154     CALL MNC_CW_ADD_VATTR_TEXT('si_QSWtave','description',
155     & 'Time-averaged net upward shortwave radiation',myThid)
156    
157     CALL MNC_CW_ADD_VNAME(
158     & 'si_HEFFtave','Cen_xy_Hn__-__t',3,4,myThid)
159     CALL MNC_CW_ADD_VATTR_TEXT(
160     & 'si_HEFFtave','units','m', myThid)
161     CALL MNC_CW_ADD_VATTR_TEXT('si_HEFFtave','description',
162     & 'Time-averaged effective ice thickness',myThid)
163    
164     CALL MNC_CW_ADD_VNAME(
165     & 'si_AREAtave','Cen_xy_Hn__-__t',3,4,myThid)
166     CALL MNC_CW_ADD_VATTR_TEXT(
167     & 'si_AREAtave','units','m^2/m^2',myThid)
168     CALL MNC_CW_ADD_VATTR_TEXT('si_AREAtave','description',
169     & 'Time-averaged fractional ice-covered area',myThid)
170 edhill 1.1
171     ENDIF
172     #endif /* ALLOW_MNC */
173     #endif /* ALLOW_SEAICE */
174    
175     RETURN
176     END
177    
  ViewVC Help
Powered by ViewVC 1.1.22