/[MITgcm]/MITgcm/pkg/thsice/thsice_map_exf.F
ViewVC logotype

Annotation of /MITgcm/pkg/thsice/thsice_map_exf.F

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

Revision 1.7 - (hide annotations) (download)
Thu Apr 4 16:44:34 2013 UTC (11 years, 2 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64g, 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: +28 -3 lines 
In case Energy-Reference-Level is used (temp_EvPrRn=0), account for
 energy content of Precip + RunOff. Assumes: 1) Rain has same temp as Air
 2) Snow has no heat capacity (consistent with seaice & thsice & exf pkgs)
 3) No distinction between sea-water Cp and fresh-water Cp
 4) Run-Off comes at the temp of surface water (with same Cp)
Note: since variations of Lf+Lv with Tsurf are neglected in Evap over ice
 and  snow (since snow has no Cp); this means that water vapor is considered
 to be released to Atmos @ 0^oC (=ERL).
1 jmc 1.7 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_map_exf.F,v 1.6 2008/11/07 20:43:16 jmc Exp $
2 mlosch 1.1 C $Name: $
3    
4     #include "THSICE_OPTIONS.h"
5     #ifdef ALLOW_EXF
6     #include "EXF_OPTIONS.h"
7     #endif
8    
9     CBOP
10     C !ROUTINE: THSICE_MAP_EXF
11     C !INTERFACE:
12     SUBROUTINE THSICE_MAP_EXF(
13 jmc 1.7 I iceMsk, locSST,
14     O totPrc, snowPrc, qPrcRnO, flxSW,
15 mlosch 1.1 I iMin,iMax,jMin,jMax, bi,bj, myThid )
16     C !DESCRIPTION: \bv
17     C *==========================================================*
18     C | S/R THSICE_MAP_EXF
19 jmc 1.5 C | Interface S/R : map Precip, Snow and shortwave fluxes
20 mlosch 1.1 C | from pkg EXF to thsice variables
21     C *==========================================================*
22     C \ev
23    
24     C !USES:
25     IMPLICIT NONE
26    
27     C == Global data ==
28     #include "SIZE.h"
29     #include "EEPARAMS.h"
30     #include "PARAMS.h"
31     #include "FFIELDS.h"
32     #ifdef ALLOW_EXF
33 jmc 1.2 # include "EXF_CONSTANTS.h"
34     # include "EXF_PARAM.h"
35     # include "EXF_FIELDS.h"
36 mlosch 1.1 #endif
37    
38     C !INPUT/OUTPUT PARAMETERS:
39     C === Routine arguments ===
40     C iceMsk :: sea-ice fraction: no ice=0, grid all ice 1 []
41 jmc 1.7 C locSST :: local Sea-Surface Temperature [deg.C]
42 mlosch 1.1 C totPrc :: Total Precipitation (including run-off) [kg/m2/s]
43     C snowPrc :: Snow Precipitation [kg/m2/s]
44 jmc 1.7 C qPrcRnO :: Energy content of Precip+RunOff (+=down) [W/m2]
45 jmc 1.5 C flxSW :: Downward short-wave surface flux (+=down) [W/m2]
46 mlosch 1.1 C iMin,iMax :: range of indices of computation domain
47     C jMin,jMax :: range of indices of computation domain
48     C bi,bj :: current tile indices
49     C myThid :: Thread no. that called this routine.
50     _RL iceMsk (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
51 jmc 1.7 _RL locSST (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
52 mlosch 1.1 _RL totPrc (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
53     _RL snowPrc(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
54 jmc 1.7 _RL qPrcRnO(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
55 mlosch 1.1 _RL flxSW (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
56     INTEGER iMin,iMax
57     INTEGER jMin,jMax
58     INTEGER bi,bj
59     INTEGER myThid
60     CEOP
61    
62     #ifdef ALLOW_EXF
63    
64     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
65     C === Local variables ===
66     C i,j :: current grid point indices
67     INTEGER i,j
68    
69     DO j = jMin, jMax
70     DO i = iMin, iMax
71 jmc 1.3 #ifdef ALLOW_ATM_TEMP
72 mlosch 1.1 totPrc(i,j) = precip(i,j,bi,bj)*rhoConstFresh
73 jmc 1.3 #endif
74 jmc 1.4 #ifdef ALLOW_RUNOFF
75 mlosch 1.1 totPrc(i,j) = totPrc(i,j) + runoff(i,j,bi,bj)*rhoConstFresh
76 jmc 1.4 #else
77     STOP 'ABNORMAL END: S/R THSICE_MAP_EXF: ALLOW_RUNOFF undef'
78 mlosch 1.1 #endif
79 jmc 1.3 #ifdef ALLOW_DOWNWARD_RADIATION
80 jmc 1.5 flxSW (i,j) = swdown(i,j,bi,bj)
81 jmc 1.3 #else
82     STOP 'ABNORMAL END: S/R THSICE_MAP_EXF: DOWNWARD_RADIATION undef'
83     #endif
84 mlosch 1.1 ENDDO
85     ENDDO
86 jmc 1.5
87 jmc 1.3 #ifdef ALLOW_ATM_TEMP
88 mlosch 1.1 IF ( snowPrecipFile .NE. ' ' ) THEN
89     DO j = jMin, jMax
90     DO i = iMin, iMax
91     snowPrc(i,j) = snowPrecip(i,j,bi,bj)*rhoConstFresh
92     ENDDO
93     ENDDO
94     ELSE
95 jmc 1.6 C If specific snow precipitiation is not available, use
96 mlosch 1.1 C precipitation when ever the air temperature is below 0 degC
97     DO j = jMin, jMax
98     DO i = iMin, iMax
99     IF ( iceMsk(i,j,bi,bj).GT.0. _d 0
100     & .AND. atemp(i,j,bi,bj).LE.cen2kel ) THEN
101 jmc 1.6 cML & .AND. atemp(i,j,bi,bj).LE.Tf0kel ) THEN
102 mlosch 1.1 snowPrc(i,j) = precip(i,j,bi,bj)*rhoConstFresh
103     ENDIF
104     ENDDO
105     ENDDO
106     ENDIF
107 jmc 1.7 IF ( temp_EvPrRn .NE. UNSET_RL ) THEN
108     C-- Account for energy content of Precip + RunOff :
109     C assume 1) rain has same temp as Air (higher altitude, e.g., 850.mb would
110     C be better); 2) Snow has no heat capacity (+ is counted separately)
111     C 3) no distinction between sea-water Cp and fresh-water Cp
112     C 4) Run-Off comes at the temp of surface water (with same Cp)
113     DO j = jMin, jMax
114     DO i = iMin, iMax
115     qPrcRnO(i,j) = HeatCapacity_Cp
116     & *( atemp(i,j,bi,bj) - cen2kel - temp_EvPrRn )
117     & *( precip(i,j,bi,bj)*rhoConstFresh - snowPrc(i,j) )
118     &
119     #ifdef ALLOW_RUNOFF
120     qPrcRnO(i,j) = qPrcRnO(i,j)
121     & + HeatCapacity_Cp
122     & *( locSST(i,j,bi,bj) - temp_EvPrRn )
123     & *runoff(i,j,bi,bj)*rhoConstFresh
124     #endif
125     ENDDO
126     ENDDO
127     ENDIF
128 jmc 1.3 #else /* ALLOW_ATM_TEMP */
129 jmc 1.6 STOP 'ABNORMAL END: S/R THSICE_MAP_EXF: ATM_TEMP undef'
130 jmc 1.3 #endif /* ALLOW_ATM_TEMP */
131 mlosch 1.1
132     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
133    
134     #endif /* ALLOW_EXF */
135    
136     RETURN
137     END
  ViewVC Help
Powered by ViewVC 1.1.22