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

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

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

Revision 1.15 - (hide 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 adcroft 1.15 C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_forcing.F,v 1.14 1998/12/09 16:11:52 adcroft Exp $
2 cnh 1.1
3 cnh 1.12 #include "CPP_OPTIONS.h"
4 cnh 1.1
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 adcroft 1.14 IMPLICIT NONE
12 cnh 1.1
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 cnh 1.6 _BARRIER
34    
35 adcroft 1.10 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 adcroft 1.15 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 adcroft 1.10 ENDDO
46     ENDDO
47     ENDDO
48     ENDDO
49 adcroft 1.15 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 adcroft 1.10
57 adcroft 1.15 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 cnh 1.6
60 adcroft 1.15 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 cnh 1.7 ENDDO
71     ENDDO
72     ENDDO
73 adcroft 1.15 ENDDO
74 cnh 1.1 C
75 adcroft 1.15 _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 cnh 1.1
82     RETURN
83     END
  ViewVC Help
Powered by ViewVC 1.1.22