/[MITgcm]/MITgcm/pkg/fizhi/fizhi_read_pickup.F
ViewVC logotype

Contents of /MITgcm/pkg/fizhi/fizhi_read_pickup.F

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

Revision 1.15 - (show annotations) (download)
Thu Aug 27 18:00:01 2009 UTC (14 years, 8 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, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint64, checkpoint65, checkpoint62, checkpoint63, 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, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint61v, checkpoint61w, checkpoint61z, checkpoint61x, checkpoint61y, HEAD
Changes since 1.14: +26 -27 lines 
use type specific (RL or RS) S/R from rw pkg instead of calling old S/R
 MDSREADFIELD,MDSWRITEFIELD from mdsio pkg.
1 C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_read_pickup.F,v 1.14 2004/12/14 16:54:08 edhill Exp $
2 C $Name: $
3
4 #include "FIZHI_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: FIZHI_READ_PICKUP
8 C !INTERFACE:
9 SUBROUTINE FIZHI_READ_PICKUP( myIter, myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | S/R FIZHI_READ_PICKUP
14 C | o Reads current fizhi state from a pickup file
15 C *==========================================================*
16 C \ev
17
18 C !USES:
19 CEOP
20 IMPLICIT NONE
21
22 C == Global variables ===
23 #include "SIZE.h"
24 #include "fizhi_SIZE.h"
25 #include "fizhi_land_SIZE.h"
26 #include "fizhi_coms.h"
27 #include "fizhi_land_coms.h"
28 #include "fizhi_io_comms.h"
29 #include "EEPARAMS.h"
30 #include "PARAMS.h"
31
32 C == Routine Arguments ==
33 C myIter :: time-step number
34 C myThid :: Number of this instance
35 INTEGER myIter
36 INTEGER myThid
37
38 C !LOCAL VARIABLES:
39 C fn :: character buffer for creating filename
40 C prec :: precision of pickup files
41 C k :: loop index
42 INTEGER i, prec
43 CHARACTER*(MAX_LEN_FNAM) fn
44 _RL temp(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy)
45 integer j,bi,bj
46
47 DO i = 1,MAX_LEN_FNAM
48 fn(i:i) = ' '
49 ENDDO
50 WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
51 prec = precFloat64
52
53 IF ( fizhi_mdsio_read_pickup ) THEN
54
55 C Read fields as consecutive records
56 CALL READ_REC_3D_RL( fn,prec,Nrphys,uphy, 1,myIter,myThid )
57 CALL READ_REC_3D_RL( fn,prec,Nrphys,vphy, 2,myIter,myThid )
58 CALL READ_REC_3D_RL( fn,prec,Nrphys,thphy, 3,myIter,myThid )
59 CALL READ_REC_3D_RL( fn,prec,Nrphys,sphy, 4,myIter,myThid )
60 CALL READ_REC_3D_RL( fn,prec,Nrphys,cldtot_lw, 5,myIter,myThid )
61 CALL READ_REC_3D_RL( fn,prec,Nrphys,cldras_lw, 6,myIter,myThid )
62 CALL READ_REC_3D_RL( fn,prec,Nrphys,cldlsp_lw, 7,myIter,myThid )
63 CALL READ_REC_3D_RL( fn,prec,Nrphys,lwlz, 8,myIter,myThid )
64 CALL READ_REC_3D_RL( fn,prec,Nrphys,cldtot_sw, 9,myIter,myThid )
65 CALL READ_REC_3D_RL( fn,prec,Nrphys,cldras_sw,10,myIter,myThid )
66 CALL READ_REC_3D_RL( fn,prec,Nrphys,cldlsp_sw,11,myIter,myThid )
67 CALL READ_REC_3D_RL( fn,prec,Nrphys,swlz, 12,myIter,myThid )
68 CALL READ_REC_3D_RL( fn,prec,Nrphys,qliqavelw,13,myIter,myThid )
69 CALL READ_REC_3D_RL( fn,prec,Nrphys,qliqavesw,14,myIter,myThid )
70 CALL READ_REC_3D_RL( fn,prec,Nrphys,fccavelw, 15,myIter,myThid )
71 CALL READ_REC_3D_RL( fn,prec,Nrphys,fccavesw, 16,myIter,myThid )
72 CALL READ_REC_3D_RL( fn,prec,Nrphys,temp, 17,myIter,myThid )
73
74 do bj = 1,myByLo(myThid), myByHi(myThid)
75 do bi = 1,myBxLo(myThid), myBxHi(myThid)
76 do j = 1-OLy,sNy+OLy
77 do i = 1-OLx,sNx+OLx
78 raincon(i,j,bi,bj) = temp(i,j,1,bi,bj)
79 rainlsp(i,j,bi,bj) = temp(i,j,2,bi,bj)
80 snowfall(i,j,bi,bj) = temp(i,j,3,bi,bj)
81 enddo
82 enddo
83
84 iras(bi,bj) = NINT(temp(1,1,4,bi,bj))
85 nlwcld(bi,bj) = NINT(temp(2,1,4,bi,bj))
86 nlwlz(bi,bj) = NINT(temp(3,1,4,bi,bj))
87 nswcld(bi,bj) = NINT(temp(4,1,4,bi,bj))
88 nswlz(bi,bj) = NINT(temp(5,1,4,bi,bj))
89 imstturbsw(bi,bj) = NINT(temp(6,1,4,bi,bj))
90 imstturblw(bi,bj) = NINT(temp(7,1,4,bi,bj))
91
92 enddo
93 enddo
94
95 ENDIF
96
97 #ifdef ALLOW_MNC
98 IF ( useMNC .and. fizhi_mnc_read_pickup ) THEN
99
100 _BEGIN_MASTER(myThid)
101
102 WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
103
104 CALL MNC_CW_RL_R('D',fn,0,0,'uphy',uphy,myThid)
105 CALL MNC_CW_RL_R('D',fn,0,0,'vphy',vphy,myThid)
106 CALL MNC_CW_RL_R('D',fn,0,0,'thphy',thphy,myThid)
107 CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
108 CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
109
110
111 _END_MASTER(myThid)
112
113 ENDIF
114 #endif
115
116 RETURN
117 END
  ViewVC Help
Powered by ViewVC 1.1.22