/[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.14 - (show annotations) (download)
Tue Dec 14 16:54:08 2004 UTC (19 years, 4 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint57t_post, checkpoint57o_post, checkpoint58e_post, checkpoint57v_post, checkpoint58u_post, checkpoint58w_post, checkpoint57m_post, checkpoint57s_post, checkpoint57k_post, checkpoint57d_post, checkpoint57g_post, checkpoint60, checkpoint61, checkpoint57b_post, checkpoint57c_pre, checkpoint58r_post, checkpoint57i_post, checkpoint57y_post, checkpoint57e_post, checkpoint58n_post, checkpoint58x_post, checkpoint57g_pre, checkpoint58t_post, checkpoint58h_post, checkpoint57y_pre, checkpoint57f_pre, checkpoint58q_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint58j_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint57r_post, checkpoint59, checkpoint58, eckpoint57e_pre, checkpoint57h_done, checkpoint58f_post, checkpoint57x_post, checkpoint57n_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58g_post, checkpoint58o_post, checkpoint57z_post, checkpoint57c_post, checkpoint58y_post, checkpoint58k_post, checkpoint58v_post, checkpoint58s_post, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint58p_post, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint57h_post
Changes since 1.13: +2 -2 lines 
 o every instance of _END_MASTER() has been found and replaced with
   _END_MASTER( myThid ) in order to satisfy certain picky Sun
   preprocessors
1 C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_read_pickup.F,v 1.13 2004/10/14 22:11:42 molod 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 C#ifdef ALLOW_MDSIO
54 IF ( fizhi_mdsio_read_pickup ) THEN
55
56 C Read fields as consecutive records
57 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)
58 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)
59 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)
60 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)
61 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_lw,5,myThid)
62 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_lw,6,myThid)
63 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_lw,7,myThid)
64 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,lwlz,8,myThid)
65 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_sw,9,myThid)
66 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_sw,10,myThid)
67 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_sw,11,myThid)
68 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,swlz,12,myThid)
69 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavelw,13,myThid)
70 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavesw,14,myThid)
71 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavelw,15,myThid)
72 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavesw,16,myThid)
73 CALL MDSREADFIELD(fn,prec,'RL',Nrphys,temp,17,myThid)
74
75 do bj = 1,myByLo(myThid), myByHi(myThid)
76 do bi = 1,myBxLo(myThid), myBxHi(myThid)
77 do j = 1-OLy,sNy+OLy
78 do i = 1-OLx,sNx+OLx
79 raincon(i,j,bi,bj) = temp(i,j,1,bi,bj)
80 rainlsp(i,j,bi,bj) = temp(i,j,2,bi,bj)
81 snowfall(i,j,bi,bj) = temp(i,j,3,bi,bj)
82 enddo
83 enddo
84
85 iras(bi,bj) = nint(temp(1,1,4,bi,bj))
86 nlwcld(bi,bj) = nint(temp(2,1,4,bi,bj))
87 nlwlz(bi,bj) = nint(temp(3,1,4,bi,bj))
88 nswcld(bi,bj) = nint(temp(4,1,4,bi,bj))
89 nswlz(bi,bj) = nint(temp(5,1,4,bi,bj))
90 imstturbsw(bi,bj) = nint(temp(6,1,4,bi,bj))
91 imstturblw(bi,bj) = nint(temp(7,1,4,bi,bj))
92
93 enddo
94 enddo
95
96 ENDIF
97
98 #ifdef ALLOW_MNC
99 IF ( useMNC .and. fizhi_mnc_read_pickup ) THEN
100
101 _BEGIN_MASTER(myThid)
102
103 WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
104
105 CALL MNC_CW_RL_R('D',fn,0,0,'uphy',uphy,myThid)
106 CALL MNC_CW_RL_R('D',fn,0,0,'vphy',vphy,myThid)
107 CALL MNC_CW_RL_R('D',fn,0,0,'thphy',thphy,myThid)
108 CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
109 CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
110
111
112 _END_MASTER(myThid)
113
114 ENDIF
115 #endif
116
117 RETURN
118 END
  ViewVC Help
Powered by ViewVC 1.1.22