/[MITgcm]/MITgcm/verification/exp2/code/external_fields_scale.F
ViewVC logotype

Contents of /MITgcm/verification/exp2/code/external_fields_scale.F

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

Revision 1.2 - (show annotations) (download)
Wed Nov 29 22:31:47 2000 UTC (23 years, 4 months ago) by adcroft
Branch: MAIN
CVS Tags: HEAD
Changes since 1.1: +0 -0 lines
FILE REMOVED 
Fixed confusion about units of forcing arrays in FFIELDS.h
namely Fu,Fv,Qnet,Qsw,EmPmR:
  - Removed verification/natl_box/code/external_fields_scale.F
        (did not differ from that in model/src)
  - Changed units of fu,fv,Qnet,Qsw,EmPmR back to proper units
     (see FFIELDS.h for description)
  - Scale fu,fv,Qnet,Qsw,EmPmR when used in external_forcing_surf.F,
    kpp_calc.F and kpp_transport_t.F
  - Removed model/src/external_fields_scale.F and calls to it
  - verification/natl_box uses flux data with "atmospheric" sign so
    a special version of external_fields_load.F is used to
    change the data as it's read in. This way, the arrays
    have the right units and signs at all times tha a user could
    possibly use them.
1
2 #include "CPP_OPTIONS.h"
3
4 CStartOfInterface
5 SUBROUTINE EXTERNAL_FIELDS_SCALE( myThid )
6 C /==========================================================\
7 C | SUBROUTINE EXTERNAL_FIELDS_SCALE |
8 C | o Rescale fields from external source . |
9 C |==========================================================|
10 C | Notes |
11 C | Conversion of flux fields are described in FFIELDS.h |
12 C \==========================================================/
13 IMPLICIT NONE
14
15 C === Global variables ===
16 #include "SIZE.h"
17 #include "EEPARAMS.h"
18 #include "PARAMS.h"
19 #include "FFIELDS.h"
20 #include "GRID.h"
21
22 C === Routine arguments ===
23 C myThid - Thread no. that called this routine.
24 INTEGER myThid
25 CEndOfInterface
26
27 C === Functions ===
28 LOGICAL DIFFERENT_MULTIPLE
29 EXTERNAL DIFFERENT_MULTIPLE
30
31 C === Local arrays ===
32
33 C === Local variables ===
34 INTEGER bi,bj,i,j
35
36 _BARRIER
37
38 DO bj = myByLo(myThid), myByHi(myThid)
39 DO bi = myBxLo(myThid), myBxHi(myThid)
40 DO j=1-oLy,sNy+oLy
41 DO i=1-oLx,sNx+oLx
42 c Zonal wind stress fu:
43 fu(i,j,bi,bj) = fu(i,j,bi,bj)
44 & *horiVertRatio*recip_rhoNil*recip_dRf(1)
45 c Meridional wind stress fv:
46 fv(i,j,bi,bj) = fv(i,j,bi,bj)
47 & *horiVertRatio*recip_rhoNil*recip_dRf(1)
48 #ifdef USE_NATURAL_BCS
49 c Freshwater flux EmPmR:
50 EmPmR(i,j,bi,bj) =
51 & -EmPmR(i,j,bi,bj)*recip_dRf(1)*salt(i,j,1,bi,bj)
52 #else
53 c Freshwater flux EmPmR:
54 EmPmR(i,j,bi,bj) =
55 & -EmPmR(i,j,bi,bj)*recip_dRf(1)*35.
56 #endif
57 c Net heat flux Qnet:
58 Qnet(i,j,bi,bj) =
59 & - Qnet(i,j,bi,bj)*recip_Cp*recip_rhoNil*recip_dRf(1)
60 #ifdef SHORTWAVE_HEATING
61 c Short wave radiative flux Qsw:
62 Qsw(i,j,bi,bj) =
63 & - Qsw(i,j,bi,bj)*recip_Cp*recip_rhoNil*recip_dRf(1)
64 #endif
65 ENDDO
66 ENDDO
67 ENDDO
68 ENDDO
69
70 RETURN
71 END
  ViewVC Help
Powered by ViewVC 1.1.22