File indexing completed on 2025-07-08 05:11:10 UTC
view on githubraw file Latest commit 00c7090d on 2025-07-07 16:10:22 UTC
cf336ab6c5 Ryan*0001 #include "LAYERS_OPTIONS.h"
0002
0003
0004
0005
0006
0007
0008
0009 SUBROUTINE LAYERS_CALC_RHS(
0010 I myThid )
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022 IMPLICIT NONE
0023
0024 #include "SIZE.h"
0025 #include "EEPARAMS.h"
0026 #include "PARAMS.h"
0027 #include "GRID.h"
50d8304171 Ryan*0028 #include "FFIELDS.h"
cf336ab6c5 Ryan*0029 #include "LAYERS_SIZE.h"
0030 #include "LAYERS.h"
0031
0032
0033
0034 INTEGER myThid
0035
0036
0037 #ifdef ALLOW_LAYERS
ee16a2cae4 Ryan*0038 #ifdef LAYERS_THERMODYNAMICS
cf336ab6c5 Ryan*0039
0040
0041
8d1543706e Jean*0042
cf336ab6c5 Ryan*0043
0044
8d1543706e Jean*0045
cf336ab6c5 Ryan*0046
0047
0048 INTEGER bi, bj
da9f56e003 Jean*0049 INTEGER i,j,k,kdown,iTracer
0050 _RL fluxfac(2), downfac, tmpfac
0051
50d8304171 Ryan*0052 _RL minusone
0053 PARAMETER (minusOne=-1.)
cf336ab6c5 Ryan*0054
0055
0056
0057 fluxfac(1) = 1.0/(HeatCapacity_Cp*rUnit2mass)
0058 fluxfac(2) = 1.0/rUnit2mass
0059
0060 DO bj = myByLo(myThid), myByHi(myThid)
8d1543706e Jean*0061 DO bi = myBxLo(myThid), myBxHi(myThid)
da9f56e003 Jean*0062
8d1543706e Jean*0063 DO iTracer = 1,2
da9f56e003 Jean*0064 k = 1
8d1543706e Jean*0065
cf336ab6c5 Ryan*0066 DO j=1-OLy,sNy+OLy
0067 DO i=1-OLx,sNx+OLx
50d8304171 Ryan*0068
0069 #ifdef SHORTWAVE_HEATING
0070
0071 IF (iTracer.EQ.1) THEN
0072 layers_surfflux(i,j,k,iTracer,bi,bj) =
0073 & layers_surfflux(i,j,k,iTracer,bi,bj)
0074
0075 & +Qsw(i,j,bi,bj)
0076 ENDIF
da9f56e003 Jean*0077 #endif /* SHORTWAVE_HEATING */
50d8304171 Ryan*0078
cf336ab6c5 Ryan*0079 layers_surfflux(i,j,k,iTracer,bi,bj) =
0080 & layers_surfflux(i,j,k,iTracer,bi,bj)
8d1543706e Jean*0081 & *recip_drF(1)*_recip_hFacC(i,j,1,bi,bj)
cf336ab6c5 Ryan*0082 & *fluxfac(iTracer)
0083 ENDDO
8d1543706e Jean*0084 ENDDO
da9f56e003 Jean*0085
cf336ab6c5 Ryan*0086
0087
0088
0089 DO k=1,Nr
00c7090dc0 Mart*0090 downFac = 1. _d 0
0091
0092 IF ( usingZCoords ) THEN
0093 kdown = MIN(k+1,Nr)
0094 IF ( k.EQ.Nr ) downFac = 0. _d 0
cf336ab6c5 Ryan*0095 ELSE
00c7090dc0 Mart*0096
0097 kdown = MAX(k-1,1)
0098 IF ( k.EQ.1 ) downFac = 0. _d 0
8d1543706e Jean*0099 ENDIF
50d8304171 Ryan*0100 DO j=1-OLy,sNy+OLy-1
0101 DO i=1-OLx,sNx+OLx-1
0102
cf336ab6c5 Ryan*0103 tmpfac = -_recip_hFacC(i,j,k,bi,bj)*recip_drF(k)
0104 & *recip_rA(i,j,bi,bj)*recip_deepFac2C(k)*recip_rhoFacC(k)
0105 layers_dfx(i,j,k,iTracer,bi,bj) = maskInC(i,j,bi,bj) *
8d1543706e Jean*0106 & tmpfac * ( layers_dfx(i+1,j,k,iTracer,bi,bj) -
cf336ab6c5 Ryan*0107 & layers_dfx(i,j,k,iTracer,bi,bj) )
0108 layers_dfy(i,j,k,iTracer,bi,bj) = maskInC(i,j,bi,bj) *
8d1543706e Jean*0109 & tmpfac * ( layers_dfy(i,j+1,k,iTracer,bi,bj) -
0110 & layers_dfy(i,j,k,iTracer,bi,bj) )
cf336ab6c5 Ryan*0111 layers_dfr(i,j,k,iTracer,bi,bj) = tmpfac * rkSign *
0112 & ( layers_dfr(i,j,kdown,iTracer,bi,bj)*downfac -
0113 & layers_dfr(i,j,k,iTracer,bi,bj) )
50d8304171 Ryan*0114
0115 layers_afx(i,j,k,iTracer,bi,bj) = maskInC(i,j,bi,bj) *
0116 & tmpfac * ( layers_afx(i+1,j,k,iTracer,bi,bj) -
0117 & layers_afx(i,j,k,iTracer,bi,bj) )
0118 layers_afy(i,j,k,iTracer,bi,bj) = maskInC(i,j,bi,bj) *
0119 & tmpfac * ( layers_afy(i,j+1,k,iTracer,bi,bj) -
0120 & layers_afy(i,j,k,iTracer,bi,bj) )
0121 layers_afr(i,j,k,iTracer,bi,bj) = tmpfac * rkSign *
0122 & ( layers_afr(i,j,kdown,iTracer,bi,bj)*downfac -
0123 & layers_afr(i,j,k,iTracer,bi,bj) )
0124
0125 #ifdef SHORTWAVE_HEATING
da9f56e003 Jean*0126 IF (iTracer.EQ.1) THEN
50d8304171 Ryan*0127 layers_sw(i,j,k,iTracer,bi,bj) =
0128 & layers_sw(i,j,k,iTracer,bi,bj)
00c7090dc0 Mart*0129 & + Qsw(i,j,bi,bj)*gravitySign
0130 & *( SWFrac3D(i,j,k,bi,bj) - SWFrac3D(i,j,k+1,bi,bj) )
50d8304171 Ryan*0131 & *fluxfac(1)
0132 & *recip_drF(k)*_recip_hFacC(i,j,k,bi,bj)
0133 ENDIF
da9f56e003 Jean*0134 #endif /* SHORTWAVE_HEATING */
0135
cf336ab6c5 Ryan*0136 ENDDO
8d1543706e Jean*0137 ENDDO
cf336ab6c5 Ryan*0138 ENDDO
0139 ENDDO
0140 ENDDO
0141 ENDDO
0142
0143
0144
0145
0146
0147
0148
0149
0150
0151
0152
0153
0154
0155
0156
0157
0158
0159
0160
0161
0162
0163
0164
0165
0166
0167
0168
0169
0170
0171
0172
8d1543706e Jean*0173
cf336ab6c5 Ryan*0174
0175
0176
0177
0178
0179
0180
0181
0182
0183
0184
8d1543706e Jean*0185
cf336ab6c5 Ryan*0186
0187
0188
0189
0190
0191
0192
0193
0194
0195
0196
8d1543706e Jean*0197
cf336ab6c5 Ryan*0198
0199
0200
0201
0202
0203
0204
0205
0206
0207
0208
0209
0210
0211
0212
0213
0214
0215
0216
0217
0218
0219
0220
0221
0222
0223
0224
0225
0226
0227
0228
0229
0230
0231
0232
0233
0234
0235
0236
0237
0238
0239
0240
8d1543706e Jean*0241
cf336ab6c5 Ryan*0242
0243 #endif /* LAYERS_THERMODYNAMICS */
ee16a2cae4 Ryan*0244 #endif /* USE_LAYERS */
cf336ab6c5 Ryan*0245 RETURN
0246 END
0247