/[MITgcm]/MITgcm/pkg/aim/phy_lscond.F
ViewVC logotype

Contents of /MITgcm/pkg/aim/phy_lscond.F

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


Revision 1.4 - (show annotations) (download)
Fri Sep 27 20:05:11 2002 UTC (21 years, 7 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint48f_post, checkpoint51k_post, checkpoint53f_post, checkpoint47j_post, checkpoint54a_pre, checkpoint55c_post, checkpoint53b_pre, checkpoint48d_pre, checkpoint51l_post, checkpoint51j_post, branch-exfmods-tag, checkpoint47e_post, checkpoint57m_post, checkpoint52l_pre, checkpoint48i_post, checkpoint52e_pre, hrcube4, hrcube5, checkpoint57g_pre, checkpoint52j_post, checkpoint47f_post, checkpoint48d_post, checkpoint51o_pre, checkpoint57f_post, checkpoint46j_post, checkpoint47c_post, checkpoint50e_post, checkpoint52e_post, checkpoint50c_post, checkpoint46i_post, checkpoint51n_pre, checkpoint47d_post, checkpoint57j_post, checkpoint47a_post, checkpoint57b_post, checkpoint52d_pre, checkpoint53c_post, checkpoint53d_post, checkpoint57f_pre, checkpoint48a_post, checkpoint55d_pre, checkpoint51f_pre, checkpoint57g_post, checkpoint48e_post, checkpoint57c_pre, checkpoint48h_post, checkpoint55j_post, checkpoint56b_post, checkpoint50c_pre, checkpoint57h_pre, branchpoint-genmake2, checkpoint46k_post, checkpoint52j_pre, checkpoint54a_post, branch-netcdf, checkpoint50d_pre, checkpoint55h_post, checkpoint51r_post, checkpoint47i_post, checkpoint52b_pre, checkpoint52n_post, checkpoint54b_post, checkpoint46l_pre, checkpoint46j_pre, checkpoint51i_post, checkpoint57e_post, checkpoint54d_post, checkpoint47h_post, checkpoint48c_post, checkpoint46l_post, checkpoint56c_post, checkpoint54e_post, checkpoint55b_post, checkpoint51e_post, checkpoint51b_post, checkpoint51l_pre, checkpoint52m_post, checkpoint51c_post, checkpoint55, checkpoint53a_post, checkpoint55a_post, checkpoint57a_post, checkpoint48, checkpoint49, checkpoint47b_post, checkpoint56, checkpoint57o_post, checkpoint55g_post, checkpoint57h_done, checkpoint51o_post, checkpoint48g_post, checkpoint57k_post, checkpoint57d_post, checkpoint55f_post, checkpoint57i_post, checkpoint51q_post, checkpoint52l_post, checkpoint52k_post, checkpoint57h_post, checkpoint57a_pre, checkpoint54, checkpoint57, checkpoint53b_post, checkpoint51, checkpoint50, checkpoint53, checkpoint52, checkpoint50d_post, checkpoint52d_post, checkpoint51b_pre, checkpoint52a_post, checkpoint47g_post, checkpoint52b_post, checkpoint53g_post, checkpoint52f_post, checkpoint57n_post, checkpoint52c_post, checkpoint46m_post, checkpoint57p_post, checkpoint51h_pre, checkpoint50g_post, checkpoint50b_pre, checkpoint51g_post, ecco_c52_e35, checkpoint54f_post, checkpoint51f_post, checkpoint48b_post, checkpoint50b_post, eckpoint57e_pre, checkpoint57c_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint52a_pre, checkpoint47d_pre, checkpoint51d_post, checkpoint48c_pre, checkpoint51m_post, checkpoint51t_post, checkpoint53d_pre, checkpoint47, checkpoint55e_post, checkpoint54c_post, checkpoint50h_post, checkpoint52i_post, checkpoint51a_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51p_post, checkpoint51n_post, checkpoint55i_post, checkpoint51i_pre, checkpoint57l_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint56a_post, checkpoint51s_post, checkpoint55d_post
Branch point for: netcdf-sm0, branch-genmake2, branch-nonh, tg2-branch, checkpoint51n_branch, branch-exfmods-curt
Changes since 1.3: +22 -14 lines
Clean up AIM package (and keep the results unchanged):
a) include CPP_OPTION and use IMPLICT NONE in all routines ;
  declare all the variables _RL ;
b) use _d 0 for all numerical constants in Physics package,
  so that the code works with g77 (and give the right answer)
c) use ifdef ALLOW_AIM everywhere so that the package can be
 compiled without increasing the memory size.
d) clean-up the AIM interface (remove commented lines, unused
  variables ...)

1 C $Header: /u/gcmpack/MITgcm/pkg/aim/phy_lscond.F,v 1.3 2001/05/29 19:28:53 cnh Exp $
2 C $Name: $
3
4 #include "AIM_OPTIONS.h"
5
6 SUBROUTINE LSCOND (PSA,QA,QSAT,
7 * PRECLS,DTLSC,DQLSC,myThid)
8 C--
9 C-- SUBROUTINE LSCOND (PSA,QA,QSAT,
10 C-- * PRECLS,DTLSC,DQLSC)
11 C--
12 C-- Purpose: Compute large-scale precipitation and
13 C-- associated tendencies of temperature and moisture
14 C-- Input: PSA = norm. surface pressure [p/p0] (2-dim)
15 C-- QA = specific humidity [g/kg] (3-dim)
16 C-- QSAT = saturation spec. hum. [g/kg] (3-dim)
17 C-- Output: PRECLS = large-scale precipitation [g/(m^2 s)] (2-dim)
18 C-- DTLSC = temperature tendency from l.s. cond (3-dim)
19 C-- DQLSC = hum. tendency [g/(kg s)] from l.s. cond (3-dim)
20 C--
21
22 IMPLICIT NONE
23
24 C Resolution parameters
25
26 C-- size for MITgcm & Physics package :
27 #include "AIM_SIZE.h"
28
29 #include "EEPARAMS.h"
30
31 #include "AIM_GRID.h"
32
33 C Physical constants + functions of sigma and latitude
34
35 #include "com_physcon.h"
36
37 C Large-scale condensation constants
38
39 #include "com_lsccon.h"
40
41 C-- Routine arguments:
42 INTEGER myThid
43 _RL PSA(NGP), QA(NGP,NLEV), QSAT(NGP,NLEV)
44 _RL PRECLS(NGP), DTLSC(NGP,NLEV), DQLSC(NGP,NLEV)
45
46 #ifdef ALLOW_AIM
47
48 C-- Local variables:
49 INTEGER J, K
50
51 C- jmc: declare all local variables:
52 _RL RTLSC, TFACT, PRG, PFACT
53 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
54
55 C-- 1. Initialization
56
57 RTLSC = 1./(TRLSC*3600)
58 TFACT = ALHC/CP
59 PRG = P0/GG
60
61 DO J=1,NGP
62 DTLSC(J,1) = 0.
63 DQLSC(J,1) = 0.
64 PRECLS(J) = 0.
65 ENDDO
66
67 C-- 2. Tendencies of temperature and moisture
68 C
69 DO K=2,NLEV
70 DO J=1,NGP
71 DQLSC(J,K) = MIN(0. _d 0,(RHLSC*QSAT(J,K)-QA(J,K)))*RTLSC
72 DTLSC(J,K) = -TFACT*DQLSC(J,K)
73 ENDDO
74 ENDDO
75
76 C-- 3. Large-scale precipitation
77
78 DO J=1,NGP
79 DO K=2,NLEVxy(J,myThid)
80 PFACT = DSIG(K)*PRG
81 PRECLS(J) = PRECLS(J)-PFACT*DQLSC(J,K)
82 ENDDO
83 ENDDO
84
85 DO J=1,NGP
86 PRECLS(J) = PRECLS(J)*PSA(J)
87 ENDDO
88
89 C--
90 #endif /* ALLOW_AIM */
91
92 RETURN
93 END

  ViewVC Help
Powered by ViewVC 1.1.22