/[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.7 - (hide annotations) (download)
Fri Nov 9 22:15:18 2012 UTC (11 years, 6 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint65, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, HEAD
Changes since 1.6: +2 -1 lines 
Merge SEAICE_SIZE.h inclusion from MITgcm_contrib/torge/itd/code/
into main branch
1 heimbach 1.7 C $Header: /u/gcmpack/MITgcm_contrib/torge/itd/code/seaice_mnc_init.F,v 1.1 2012/10/24 21:48:53 torge 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 heimbach 1.7 #include "SEAICE_SIZE.h"
23 mlosch 1.4 #include "SEAICE_PARAMS.h"
24     #endif /* SEAICE_ALLOW_EVP */
25 edhill 1.1
26     C !INPUT PARAMETERS:
27     INTEGER myThid
28     CEOP
29    
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 mlosch 1.6 & 'si_EmPmRtave','units','kg/m^2/s',myThid)
140 edhill 1.2 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    
174     RETURN
175     END
176    
  ViewVC Help
Powered by ViewVC 1.1.22