/[MITgcm]/MITgcm/model/src/ini_forcing.F
ViewVC logotype

Contents of /MITgcm/model/src/ini_forcing.F

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

Revision 1.15 - (show annotations) (download)
Tue Dec 15 00:20:34 1998 UTC (25 years, 5 months ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint19, checkpoint20
Changes since 1.14: +33 -48 lines 
 o Added "natural BCs" as alternative to "virtual salt flux"
 o Re-difined precFloat32 and precFloat64 to be 32 and 64
   so that their values can be meaningfuly set in the data file
 o Modified read_write.F to create an exception if readBinaryPrec
   is not set
 o Replaced CPP control of viscous BCs with run-time control
 o Tidied up input-data precision (ie. ini_depths cnh_dbg...)
 o ini_forcing.F now initialises *all* forcing arrays to zero
 o Definitively tested verification experiments 0,1,2 and 4
   (3 is atmospheric set-up which is in a state of flux)
1 C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_forcing.F,v 1.14 1998/12/09 16:11:52 adcroft Exp $
2
3 #include "CPP_OPTIONS.h"
4
5 CStartOfInterface
6 SUBROUTINE INI_FORCING( myThid )
7 C /==========================================================\
8 C | SUBROUTINE INI_FORCING |
9 C | o Set model initial forcing fields. |
10 C \==========================================================/
11 IMPLICIT NONE
12
13 C === Global variables ===
14 #include "SIZE.h"
15 #include "EEPARAMS.h"
16 #include "PARAMS.h"
17 #include "GRID.h"
18 #include "FFIELDS.h"
19
20 C == Routine arguments ==
21 C myThid - Number of this instance of INI_FORCING
22 INTEGER myThid
23 CEndOfInterface
24
25 C == Local variables ==
26 C iG, jG - Global coordinate index
27 C bi,bj - Loop counters
28 C I,J,K
29 INTEGER iG, jG
30 INTEGER bi, bj
31 INTEGER I, J, K
32
33 _BARRIER
34
35 DO bj = myByLo(myThid), myByHi(myThid)
36 DO bi = myBxLo(myThid), myBxHi(myThid)
37 DO j=1-OLy,sNy+OLy
38 DO i=1-OLx,sNx+OLx
39 fu (i,j,bi,bj) = 0. _d 0
40 fv (i,j,bi,bj) = 0. _d 0
41 Qnet (i,j,bi,bj) = 0. _d 0
42 EmPmR(i,j,bi,bj) = 0. _d 0
43 SST (i,j,bi,bj) = 0. _d 0
44 SSS (i,j,bi,bj) = 0. _d 0
45 ENDDO
46 ENDDO
47 ENDDO
48 ENDDO
49 C
50 _EXCH_XY_R4(fu , myThid )
51 _EXCH_XY_R4(fv , myThid )
52 _EXCH_XY_R4(Qnet , myThid )
53 _EXCH_XY_R4(EmPmR, myThid )
54 _EXCH_XY_R4(SST , myThid )
55 _EXCH_XY_R4(SSS , myThid )
56
57 C CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)
58 C CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)
59
60 DO bj = myByLo(myThid), myByHi(myThid)
61 DO bi = myBxLo(myThid), myBxHi(myThid)
62 DO j=1-OLy,sNy+OLy
63 DO i=1-OLx,sNx+OLx
64 fudot (i,j,bi,bj) = 0. _d 0
65 fvdot (i,j,bi,bj) = 0. _d 0
66 Qnetdot (i,j,bi,bj) = 0. _d 0
67 EmPmRdot(i,j,bi,bj) = 0. _d 0
68 SSTdot (i,j,bi,bj) = 0. _d 0
69 SSSdot (i,j,bi,bj) = 0. _d 0
70 ENDDO
71 ENDDO
72 ENDDO
73 ENDDO
74 C
75 _EXCH_XY_R4(fudot , myThid )
76 _EXCH_XY_R4(fvdot , myThid )
77 _EXCH_XY_R4(Qnetdot , myThid )
78 _EXCH_XY_R4(EmPmRdot, myThid )
79 _EXCH_XY_R4(SSTdot , myThid )
80 _EXCH_XY_R4(SSSdot , myThid )
81
82 RETURN
83 END
  ViewVC Help
Powered by ViewVC 1.1.22