Changeset 2528 for trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diafwb.F90
- Timestamp:
- 2010-12-27T18:33:53+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diafwb.F90
- Property svn:eol-style deleted
r1581 r2528 8 8 !! 9.0 ! 05-11 (V. Garnier) Surface pressure gradient organization 9 9 !!---------------------------------------------------------------------- 10 #if ( defined key_orca_r2 || defined key_orca_r4 )&& ! defined key_coupled11 !!---------------------------------------------------------------------- 12 !! "key_orca_r 2 or 4"10 #if ( defined key_orca_r1 || defined key_orca_r2 || defined key_orca_r4 ) && ! defined key_coupled 11 !!---------------------------------------------------------------------- 12 !! "key_orca_r1 or 2 or 4" 13 13 !!---------------------------------------------------------------------- 14 14 !!---------------------------------------------------------------------- … … 30 30 LOGICAL, PUBLIC, PARAMETER :: lk_diafwb = .TRUE. !: fresh water budget flag 31 31 32 REAL(wp) :: a_ emp, &32 REAL(wp) :: a_fwf , & 33 33 & a_sshb, a_sshn, a_salb, a_saln 34 34 REAL(wp), DIMENSION(4) :: a_flxi, a_flxo, a_temi, a_temo, a_sali, a_salo … … 38 38 # include "vectopt_loop_substitute.h90" 39 39 !!---------------------------------------------------------------------- 40 !! OPA 9.0 , LOCEAN-IPSL (2006)40 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 41 41 !! $Id$ 42 !! Software governed by the CeCILL licence ( modipsl/doc/NEMO_CeCILL.txt)42 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 43 43 !!---------------------------------------------------------------------- 44 44 … … 59 59 REAL(wp) :: ztemi(4), ztemo(4), zsali(4), zsalo(4), zflxi(4), zflxo(4) 60 60 REAL(wp) :: zt, zs, zu 61 REAL(wp) :: zsm0, z empnew61 REAL(wp) :: zsm0, zfwfnew 62 62 !!---------------------------------------------------------------------- 63 63 … … 65 65 zsm0 = 34.72654 66 66 67 ! To compute emp mean value mean emp67 ! To compute fwf mean value mean fwf 68 68 69 69 IF( kt == nit000 ) THEN 70 70 71 a_ emp= 0.e071 a_fwf = 0.e0 72 72 a_sshb = 0.e0 ! valeur de ssh au debut de la simulation 73 73 a_salb = 0.e0 ! valeur de sal au debut de la simulation … … 87 87 ENDIF 88 88 89 a_ emp = SUM( e1t(:,:) * e2t(:,:) * emp (:,:) * tmask_i(:,:) )90 IF( lk_mpp ) CALL mpp_sum( a_ emp) ! sum over the global domain89 a_fwf = SUM( e1t(:,:) * e2t(:,:) * ( emp(:,:)-rnf(:,:) ) * tmask_i(:,:) ) 90 IF( lk_mpp ) CALL mpp_sum( a_fwf ) ! sum over the global domain 91 91 92 92 IF( kt == nitend ) THEN … … 95 95 zarea = 0.e0 96 96 zvol = 0.e0 97 z empnew = 0.e097 zfwfnew = 0.e0 98 98 ! Mean sea level at nitend 99 99 a_sshn = SUM( e1t(:,:) * e2t(:,:) * sshn(:,:) * tmask_i(:,:) ) … … 115 115 116 116 ! Conversion in m3 117 a_ emp = a_emp* rdttra(1) * 1.e-3117 a_fwf = a_fwf * rdttra(1) * 1.e-3 118 118 119 ! empcorrection to bring back the mean ssh to zero120 z empnew = a_sshn / ( ( nitend - nit000 + 1 ) * rdt ) * 1.e3 / zarea119 ! fwf correction to bring back the mean ssh to zero 120 zfwfnew = a_sshn / ( ( nitend - nit000 + 1 ) * rdt ) * 1.e3 / zarea 121 121 122 122 ENDIF … … 162 162 ii0 = 140 ; ii1 = 140 163 163 ij0 = 102 ; ij1 = 102 164 ! ! ======================= 165 CASE ( 1 ) ! ORCA_R1 configurations 166 ! ! ======================= 167 ii0 = 283 ; ii1 = 283 168 ij0 = 200 ; ij1 = 200 164 169 ! ! ======================= 165 170 CASE DEFAULT ! ORCA R05 or R025 … … 204 209 ii0 = 137 ; ii1 = 137 205 210 ij0 = 101 ; ij1 = 102 211 ! ! ======================= 212 CASE ( 1 ) ! ORCA_R1 configurations 213 ! ! ======================= 214 ii0 = 282 ; ii1 = 282 215 ij0 = 200 ; ij1 = 200 206 216 ! ! ======================= 207 217 CASE DEFAULT ! ORCA R05 or R025 … … 247 257 ij0 = 88 ; ij1 = 88 248 258 ! ! ======================= 259 CASE ( 1 ) ! ORCA_R1 configurations 260 ! ! ======================= 261 ii0 = 331 ; ii1 = 331 262 ij0 = 176 ; ij1 = 176 263 ! ! ======================= 249 264 CASE DEFAULT ! ORCA R05 or R025 250 265 ! ! ======================= … … 289 304 ij0 = 116 ; ij1 = 116 290 305 ! ! ======================= 306 CASE ( 1 ) ! ORCA_R1 configurations 307 ! ! ======================= 308 ii0 = 297 ; ii1 = 297 309 ij0 = 230 ; ij1 = 230 310 ! ! ======================= 291 311 CASE DEFAULT ! ORCA R05 or R025 292 312 ! ! ======================= … … 362 382 WRITE(inum,*) 363 383 WRITE(inum,*) 'Net freshwater budget ' 364 WRITE(inum,9010) ' emp = ',a_emp, ' m3 =', a_emp/(FLOAT(nitend-nit000+1)*rdttra(1)) * 1.e-6,' Sv'384 WRITE(inum,9010) ' fwf = ',a_fwf, ' m3 =', a_fwf /(FLOAT(nitend-nit000+1)*rdttra(1)) * 1.e-6,' Sv' 365 385 WRITE(inum,*) 366 386 WRITE(inum,9010) ' zarea =',zarea
Note: See TracChangeset
for help on using the changeset viewer.