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

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

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


Revision 1.2 - (show annotations) (download)
Wed Aug 10 19:18:52 2005 UTC (18 years, 9 months ago) by edhill
Branch: MAIN
Changes since 1.1: +70 -60 lines
 o fix SEAICE use of MNC for tave and snapshot output

1 C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_mnc_init.F,v 1.1 2005/04/03 05:16:44 edhill Exp $
2 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
22 C !INPUT PARAMETERS:
23 INTEGER myThid
24 CEOP
25
26 #ifdef ALLOW_SEAICE
27 #ifdef ALLOW_MNC
28 IF (useMNC) THEN
29
30 CALL MNC_CW_ADD_VNAME('si_UWIND','U_xy_Hn__-__t',3,4,myThid)
31 CALL MNC_CW_ADD_VATTR_TEXT('si_UWIND','units','m/s', myThid)
32 CALL MNC_CW_ADD_VATTR_TEXT('si_UWIND','description',
33 & 'Surface (10m) zonal wind velocity', myThid)
34
35 CALL MNC_CW_ADD_VNAME('si_VWIND','V_xy_Hn__-__t',3,4,myThid)
36 CALL MNC_CW_ADD_VATTR_TEXT('si_VWIND','units','m/s', myThid)
37 CALL MNC_CW_ADD_VATTR_TEXT('si_VWIND','description',
38 & 'Surface (10m) meridional wind velocity', myThid)
39
40 CALL MNC_CW_ADD_VNAME('si_EVAP','Cen_xy_Hn__-__t',3,4,myThid)
41 CALL MNC_CW_ADD_VATTR_TEXT('si_EVAP','units','m/s', myThid)
42 CALL MNC_CW_ADD_VATTR_TEXT('si_EVAP','description',
43 & 'Evaporation', myThid)
44
45 CALL MNC_CW_ADD_VNAME('si_PRECIP','Cen_xy_Hn__-__t',3,4,myThid)
46 CALL MNC_CW_ADD_VATTR_TEXT('si_PRECIP','units','m/s', myThid)
47 CALL MNC_CW_ADD_VATTR_TEXT('si_PRECIP','description',
48 & 'Precipitation', myThid)
49
50 CALL MNC_CW_ADD_VNAME('si_RUNOFF','Cen_xy_Hn__-__t',3,4,myThid)
51 CALL MNC_CW_ADD_VATTR_TEXT('si_RUNOFF','units','m/s', myThid)
52 CALL MNC_CW_ADD_VATTR_TEXT('si_RUNOFF','description',
53 & 'River and glacier runoff', myThid)
54
55 CALL MNC_CW_ADD_VNAME('si_ATEMP','Cen_xy_Hn__-__t',3,4,myThid)
56 CALL MNC_CW_ADD_VATTR_TEXT('si_ATEMP','units','deg K', myThid)
57 CALL MNC_CW_ADD_VATTR_TEXT('si_ATEMP','description',
58 & 'Surface (2m) air temperature', myThid)
59
60 CALL MNC_CW_ADD_VNAME('si_AQH','Cen_xy_Hn__-__t',3,4,myThid)
61 CALL MNC_CW_ADD_VATTR_TEXT('si_AQH','units','kg/kg', myThid)
62 CALL MNC_CW_ADD_VATTR_TEXT('si_AQH','description',
63 & 'Surface (2m) specific humidity', myThid)
64
65 CALL MNC_CW_ADD_VNAME('si_SWDOWN','Cen_xy_Hn__-__t',3,4,myThid)
66 CALL MNC_CW_ADD_VATTR_TEXT('si_SWDOWN','units','W/m^2', myThid)
67 CALL MNC_CW_ADD_VATTR_TEXT('si_SWDOWN','description',
68 & 'Downward shortwave radiation', myThid)
69
70 CALL MNC_CW_ADD_VNAME('si_LWDOWN','Cen_xy_Hn__-__t',3,4,myThid)
71 CALL MNC_CW_ADD_VATTR_TEXT('si_LWDOWN','units','W/m^2', myThid)
72 CALL MNC_CW_ADD_VATTR_TEXT('si_LWDOWN','description',
73 & 'Downward longwave radiation', myThid)
74
75 C Time averages
76 CALL MNC_CW_ADD_VNAME('si_UICEtave','U_xy_Hn__-__t',3,4,myThid)
77 CALL MNC_CW_ADD_VATTR_TEXT('si_UICEtave','units','m/s', myThid)
78 CALL MNC_CW_ADD_VATTR_TEXT('si_UICEtave','description',
79 & 'Time-averaged ', myThid)
80
81 CALL MNC_CW_ADD_VNAME('si_VICEtave','V_xy_Hn__-__t',3,4,myThid)
82 CALL MNC_CW_ADD_VATTR_TEXT('si_VICEtave','units','m/s', myThid)
83 CALL MNC_CW_ADD_VATTR_TEXT('si_VICEtave','description',
84 & 'Time-averaged ', myThid)
85
86 CALL MNC_CW_ADD_VNAME('si_FUtave','U_xy_Hn__-__t',3,4,myThid)
87 CALL MNC_CW_ADD_VATTR_TEXT('si_FUtave','units','N/m^2', myThid)
88 CALL MNC_CW_ADD_VATTR_TEXT('si_FUtave','description',
89 & 'Time-averaged zonal surface wind stress', myThid)
90
91 CALL MNC_CW_ADD_VNAME('si_FVtave','V_xy_Hn__-__t',3,4,myThid)
92 CALL MNC_CW_ADD_VATTR_TEXT('si_FVtave','units','N/m^2',myThid)
93 CALL MNC_CW_ADD_VATTR_TEXT('si_FVtave','description',
94 & 'Time-averaged Meridional surface wind stress', myThid)
95
96 CALL MNC_CW_ADD_VNAME(
97 & 'si_EmPmRtave','Cen_xy_Hn__-__t',3,4,myThid)
98 CALL MNC_CW_ADD_VATTR_TEXT(
99 & 'si_EmPmRtave','units','m/s',myThid)
100 CALL MNC_CW_ADD_VATTR_TEXT('si_EmPmRtave','description',
101 & 'Time-averaged net upward freshwater flux', myThid)
102
103 CALL MNC_CW_ADD_VNAME(
104 & 'si_QNETtave','Cen_xy_Hn__-__t',3,4,myThid)
105 CALL MNC_CW_ADD_VATTR_TEXT(
106 & 'si_QNETtave','units','W/m^2', myThid)
107 CALL MNC_CW_ADD_VATTR_TEXT('si_QNETtave','description',
108 & 'Time-averaged net upward surface heat flux', myThid)
109
110 CALL MNC_CW_ADD_VNAME(
111 & 'si_QSWtave','Cen_xy_Hn__-__t',3,4,myThid)
112 CALL MNC_CW_ADD_VATTR_TEXT(
113 & 'si_QSWtave','units','W/m^2',myThid)
114 CALL MNC_CW_ADD_VATTR_TEXT('si_QSWtave','description',
115 & 'Time-averaged net upward shortwave radiation',myThid)
116
117 CALL MNC_CW_ADD_VNAME(
118 & 'si_HEFFtave','Cen_xy_Hn__-__t',3,4,myThid)
119 CALL MNC_CW_ADD_VATTR_TEXT(
120 & 'si_HEFFtave','units','m', myThid)
121 CALL MNC_CW_ADD_VATTR_TEXT('si_HEFFtave','description',
122 & 'Time-averaged effective ice thickness',myThid)
123
124 CALL MNC_CW_ADD_VNAME(
125 & 'si_AREAtave','Cen_xy_Hn__-__t',3,4,myThid)
126 CALL MNC_CW_ADD_VATTR_TEXT(
127 & 'si_AREAtave','units','m^2/m^2',myThid)
128 CALL MNC_CW_ADD_VATTR_TEXT('si_AREAtave','description',
129 & 'Time-averaged fractional ice-covered area',myThid)
130
131 ENDIF
132 #endif /* ALLOW_MNC */
133 #endif /* ALLOW_SEAICE */
134
135 RETURN
136 END
137

  ViewVC Help
Powered by ViewVC 1.1.22