Changeset 183 for trunk/NEMO/OPA_SRC/DYN
- Timestamp:
- 2004-11-05T15:30:14+01:00 (19 years ago)
- Location:
- trunk/NEMO/OPA_SRC/DYN
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/dynspg_fsc.F90
r119 r183 115 115 z2dt, z2dtg, zraur, znugdt, & ! temporary scalars 116 116 znurau, zssha, zspgu, zspgv, & ! " " 117 z egu, zegv, zgcb, zbtd,& ! " "117 zgcb, zbtd, & ! " " 118 118 ztdgu, ztdgv ! " " 119 119 !!---------------------------------------------------------------------- … … 153 153 zspgu = - grav * ( sshn(ji+1,jj) - sshn(ji,jj) ) / e1u(ji,jj) 154 154 zspgv = - grav * ( sshn(ji,jj+1) - sshn(ji,jj) ) / e2v(ji,jj) 155 zegu = + znurau * ( emp (ji+1,jj) - emp (ji,jj) ) / e1u(ji,jj) 156 zegv = + znurau * ( emp (ji,jj+1) - emp (ji,jj) ) / e2v(ji,jj) 157 spgu(ji,jj) = zspgu + zegu 158 spgv(ji,jj) = zspgv + zegv 155 spgu(ji,jj) = zspgu 156 spgv(ji,jj) = zspgv 159 157 #if defined key_trddyn || defined key_trd_vor 160 158 ! save the surface pressure gradient trends 161 159 utrd2(ji,jj,1) = zspgu 162 160 vtrd2(ji,jj,1) = zspgv 163 ! save the first part of the additional force trend164 utrd2(ji,jj,2) = zegu165 vtrd2(ji,jj,2) = zegv166 161 #endif 167 162 END DO -
trunk/NEMO/OPA_SRC/DYN/dynspg_fsc_atsk.F90
r109 r183 116 116 REAL(wp) :: & ! temporary scalars 117 117 z2dt, z2dtg, zraur, znugdt, znurau, & 118 zssha, zspgu, zspgv, zegu, zegv,&118 zssha, zspgu, zspgv, & 119 119 zgcb, zbtd, ztdgu, ztdgv, zgwgt 120 120 !!---------------------------------------------------------------------- … … 156 156 zspgu = - grav * ( sshn(ji+1,jj) - sshn(ji,jj) ) / e1u(ji,jj) 157 157 zspgv = - grav * ( sshn(ji,jj+1) - sshn(ji,jj) ) / e2v(ji,jj) 158 zegu = + znurau * ( emp (ji+1,jj) - emp (ji,jj) ) / e1u(ji,jj) 159 zegv = + znurau * ( emp (ji,jj+1) - emp (ji,jj) ) / e2v(ji,jj) 160 spgu(ji,jj) = zspgu + zegu 161 spgv(ji,jj) = zspgv + zegv 158 spgu(ji,jj) = zspgu 159 spgv(ji,jj) = zspgv 162 160 #if defined key_trddyn || defined key_trd_vor 163 161 ! save the surface pressure gradient trends 164 162 utrd2(ji,jj,1) = zspgu 165 163 vtrd2(ji,jj,1) = zspgv 166 ! save the first part of the additional force trend167 utrd2(ji,jj,2) = zegu168 vtrd2(ji,jj,2) = zegv169 164 #endif 170 165 END DO
Note: See TracChangeset
for help on using the changeset viewer.