/[MITgcm]/MITgcm/pkg/land/land_read_pickup.F
ViewVC logotype

Contents of /MITgcm/pkg/land/land_read_pickup.F

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

Revision 1.2 - (show annotations) (download)
Thu Mar 11 14:42:00 2004 UTC (20 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57o_post, checkpoint52n_post, checkpoint53d_post, checkpoint54a_pre, checkpoint57m_post, checkpoint55c_post, checkpoint54e_post, checkpoint54a_post, checkpoint53c_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57e_post, checkpoint52l_post, checkpoint55h_post, checkpoint57g_pre, checkpoint54b_post, checkpoint53b_pre, checkpoint55b_post, checkpoint54d_post, checkpoint56c_post, checkpoint52m_post, checkpoint55, checkpoint53a_post, checkpoint57f_pre, checkpoint57a_post, checkpoint54, checkpoint54f_post, checkpoint53b_post, checkpoint55g_post, checkpoint55f_post, checkpoint57r_post, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, checkpoint53, eckpoint57e_pre, checkpoint57h_done, checkpoint53g_post, checkpoint57n_post, checkpoint57p_post, checkpoint57f_post, checkpoint57q_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: +15 -4 lines 
new land formulation:
a) use ground enthalpy as prognostic variable to ensure exact
   energy conservation.
b) account for water temperature and for latent heat of freezing
   in all processes (rain, run-off, ground storage)
c) compute surface and ground temperature implicitly.
1 C $Header: /u/gcmpack/MITgcm/pkg/land/land_read_pickup.F,v 1.1 2003/06/12 17:54:22 jmc Exp $
2 C $Name: $
3
4 #include "LAND_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: LAND_READ_PICKUP
8 C !INTERFACE:
9 SUBROUTINE LAND_READ_PICKUP( myIter, myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | S/R LAND_READ_PICKUP
14 C | o Reads current state of land model from a pickup file
15 C *==========================================================*
16 C \ev
17
18 C !USES:
19 IMPLICIT NONE
20
21 C == Global variables ===
22 #include "LAND_SIZE.h"
23 #include "EEPARAMS.h"
24 #include "PARAMS.h"
25 #include "LAND_PARAMS.h"
26 #include "LAND_VARS.h"
27
28 C == Routine Arguments ==
29 C myIter :: time-step number
30 C myThid :: Number of this instance
31 INTEGER myIter
32 INTEGER myThid
33
34 #ifdef ALLOW_LAND
35
36 C !LOCAL VARIABLES:
37 C fn :: character buffer for creating filename
38 C prec :: precision of pickup files
39 C k :: loop index
40 INTEGER prec, k
41 CHARACTER*(MAX_LEN_FNAM) fn
42 CEOP
43
44 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
45
46 WRITE(fn,'(A,I10.10)') 'pickup_land.',myIter
47 prec = precFloat64
48
49 IF ( land_oldPickup ) THEN
50 C- Read fields as consecutive records
51 CALL MDSREADFIELD(fn,prec,'RL',land_nLev,land_groundT,1,myThid)
52 CALL MDSREADFIELD(fn,prec,'RL',land_nLev,land_groundW,2,myThid)
53
54 ELSE
55
56 C- Read fields as consecutive records
57 CALL MDSREADFIELD(fn,prec,'RL',land_nLev,land_enthalp,1,myThid)
58 CALL MDSREADFIELD(fn,prec,'RL',land_nLev,land_groundW,2,myThid)
59 k=2*land_nLev
60 CALL MDSREADFIELD(fn,prec,'RL', 1, land_skinT, k+1, myThid)
61 CALL MDSREADFIELD(fn,prec,'RL', 1, land_hSnow, k+2, myThid)
62 CALL MDSREADFIELD(fn,prec,'RL', 1, land_snowAge,k+3, myThid)
63
64 ENDIF
65
66 C- jmc: exchange is not really necessary for land model
67 C- and presently exchange S/R cannot work for 3d array where 3rd dim <> Nr
68
69 #endif /* ALLOW_LAND */
70
71 RETURN
72 END
  ViewVC Help
Powered by ViewVC 1.1.22