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

Contents of /MITgcm/pkg/thsice/thsice_get_bulkf.F

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

Revision 1.2 - (show annotations) (download)
Wed Apr 7 23:40:34 2004 UTC (20 years, 1 month ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint57o_post, checkpoint52n_post, checkpoint53d_post, checkpoint54a_pre, checkpoint57m_post, checkpoint55c_post, checkpoint54e_post, checkpoint57s_post, checkpoint54a_post, checkpoint53c_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint57e_post, checkpoint55h_post, checkpoint53b_post, checkpoint57g_pre, checkpoint54b_post, checkpoint53b_pre, checkpoint55b_post, checkpoint54d_post, checkpoint56c_post, checkpoint52m_post, checkpoint57y_pre, checkpoint55, checkpoint53a_post, checkpoint57f_pre, checkpoint57a_post, checkpoint54, checkpoint54f_post, checkpoint57v_post, checkpoint55g_post, checkpoint55f_post, checkpoint57r_post, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, checkpoint53, eckpoint57e_pre, checkpoint57h_done, checkpoint53g_post, checkpoint57x_post, checkpoint57n_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint57q_post, checkpoint57z_post, checkpoint57c_post, checkpoint55e_post, checkpoint53f_post, checkpoint55a_post, checkpoint53d_pre, checkpoint54c_post, checkpoint57j_post, checkpoint57h_pre, checkpoint57l_post, checkpoint57h_post, checkpoint56a_post, checkpoint55d_post
Changes since 1.1: +22 -11 lines 
major changes in pkg/thsice: allows atmospheric model (AIM) to use thsice.
- split thsice_therm.F in 2 S/R: thsice_solve4temp.F & thsice_calc_thickn.F
- move most of the ocean & bulk_force interface in thsice_main.F
- add a "slab ocean" component to be used with atmospheric model
1 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_get_bulkf.F,v 1.1 2003/11/23 01:20:13 jmc Exp $
2 C $Name: $
3
4 #include "THSICE_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: THSICE_GET_BULKF
8 C !INTERFACE:
9 SUBROUTINE THSICE_GET_BULKF(
10 I iceornot, Tsf,
11 O flxExceptSw, df0dT, evap, dEvdT,
12 I i,j,bi,bj,myThid )
13 C !DESCRIPTION: \bv
14 C *==========================================================*
15 C | S/R THSICE_GET_BULKF
16 C *==========================================================*
17 C | Interface S/R : get Surface Fluxes from pkg BULK_FORCE
18 C *==========================================================*
19 C \ev
20
21 C !USES:
22 IMPLICIT NONE
23
24 C == Global data ==
25 #ifdef ALLOW_BULK_FORCE
26 #include "SIZE.h"
27 #include "EEPARAMS.h"
28 #include "BULKF.h"
29 #endif
30
31 C !INPUT/OUTPUT PARAMETERS:
32 C === Routine arguments ===
33 C iceornot :: 0=open water, 1=ice cover
34 C Tsf :: surface (ice or snow) temperature (oC)
35 C flxExceptSw :: net (downward) surface heat flux, except short-wave [W/m2]
36 C df0dT :: deriv of flx with respect to Tsf [W/m/K]
37 C evap :: surface evaporation (>0 if evaporate) [kg/m2/s]
38 C dEvdT :: deriv of evap. with respect to Tsf [kg/m2/s/K]
39 C i,j, bi,bj :: current grid point indices
40 C myThid :: Thread no. that called this routine.
41 INTEGER i,j, bi,bj
42 INTEGER myThid
43 INTEGER iceornot
44 _RL Tsf
45 _RL flxExceptSw
46 _RL df0dT
47 _RL evap
48 _RL dEvdT
49 CEOP
50
51 #ifdef ALLOW_THSICE
52 #ifdef ALLOW_BULK_FORCE
53
54 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
55 C === Local variables ===
56
57 _RL flwup ! upward LW at surface (W m-2)
58 _RL flwNet_dwn ! net (downward) LW at surface (W m-2)
59 _RL fsh ! surface downward sensible heat (W m-2)
60 _RL flh ! surface downward latent heat (W m-2)
61 _RL ust, vst, ssq
62
63 ust = 0.
64 vst = 0.
65 ssq = 0.
66
67 CALL BULKF_FORMULA_LANL(
68 I uwind(i,j,bi,bj), vwind(i,j,bi,bj), wspeed(i,j,bi,bj),
69 I Tair(i,j,bi,bj), Qair(i,j,bi,bj), cloud(i,j,bi,bj), Tsf,
70 O flwup, flh, fsh, df0dT, ust, vst, evap, ssq, dEvdT,
71 I iceornot, readwindstress)
72
73 flwNet_dwn = flw(i,j,bi,bj) - flwup
74 flxExceptSw = flwNet_dwn + fsh + flh
75
76 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
77
78 #endif /* ALLOW_BULK_FORCE */
79 #endif /* ALLOW_THSICE */
80
81 RETURN
82 END
  ViewVC Help
Powered by ViewVC 1.1.22