Changeset 129 for trunk/NEMO/OPA_SRC/TRD
- Timestamp:
- 2004-07-07T14:40:30+02:00 (20 years ago)
- Location:
- trunk/NEMO/OPA_SRC/TRD
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/TRD/trdvor.F90
r109 r129 12 12 !!---------------------------------------------------------------------- 13 13 !! * Modules used 14 USE oce ! ocean dynamics and tracers variables 14 15 USE dom_oce ! ocean space and time domain variables 15 16 USE oce ! ocean dynamics and tracers variables 17 USE trddyn_oce ! ocean active tracer trend variables 16 USE trddyn_oce ! ocean active tracer trend variables 18 17 USE zdf_oce ! ocean vertical physics 19 18 USE in_out_manager ! I/O manager 20 19 21 20 USE phycst ! Define parameters for the routines 22 USE ldfdyn_oce ! ocean active tracers: lateral physics21 USE ldfdyn_oce ! ocean active tracers: lateral physics 23 22 USE daymod ! calandar 24 23 USE dianam ! build the name of file (routine) 25 USE ldfslp 24 USE ldfslp ! iso-neutral slopes 26 25 USE zdfmxl 27 26 USE ioipsl … … 30 29 IMPLICIT NONE 31 30 PRIVATE 32 !! * Accessibility 31 32 33 !! * Accessibility 33 34 PUBLIC trd_vor ! routine called by step.F90 34 35 … … 37 38 38 39 !! * Module variables 39 40 40 INTEGER,PARAMETER :: jplvor=11 41 42 43 INTEGER :: & 44 nh_t,nmoydpvor, & ! ??? 45 nidvor,nhoridvor, & 46 ndexvor1(jpi*jpj),ndimvor1 41 INTEGER :: & 42 nh_t, nmoydpvor , & 43 nidvor, nhoridvor, & 44 ndexvor1(jpi*jpj), & 45 ndimvor1 47 46 48 47 REAL(wp), DIMENSION(jpi,jpj) :: & 49 vor_avr, & ! average 50 vor_avrb, & ! before vorticity (kt-1) 51 vor_avrbb, & ! vorticity at begining of the nwrite-1 timestep averaging period 52 vor_avrbn, & ! after vorticity at time step after the 53 rotot ! begining of the NWRITE-1 timesteps 54 55 REAL(wp), DIMENSION(jpi,jpj) :: & 56 udpvor , & ! total cumulative trends 57 vdpvor 58 59 REAL(wp), DIMENSION(jpi,jpj,jplvor)::vortrd ! curl of trends 48 vor_avr , & ! average 49 vor_avrb , & ! before vorticity (kt-1) 50 vor_avrbb , & ! vorticity at begining of the nwrite-1 timestep averaging period 51 vor_avrbn , & ! after vorticity at time step after the 52 rotot , & ! begining of the NWRITE-1 timesteps 53 udpvor , & ! total cumulative trends 54 vdpvor ! " " " 55 56 REAL(wp), DIMENSION(jpi,jpj,jplvor):: & !: curl of trends 57 vortrd 60 58 61 59 !! * Substitutions … … 65 63 66 64 !!---------------------------------------------------------------------- 67 !! 200465 !! OPA 9.0 , LODYC-IPS (2004) 68 66 !!---------------------------------------------------------------------- 69 67 … … 71 69 72 70 SUBROUTINE trd_vor( kt ) 73 !!----------------------------------------- vortrd-----------------------------71 !!---------------------------------------------------------------------------- 74 72 !! *** ROUTINE trd_vor *** 75 73 !! … … 152 150 ENDIF 153 151 154 155 152 ! cumulated trends array init 156 153 nmoydpvor = 0 … … 167 164 168 165 IF( idebug /= 0 ) THEN 169 WRITE(numout,*) ' debuging trd_vor: 0. done '170 CALL FLUSH(numout)166 WRITE(numout,*) ' debuging trd_vor: 0. done ' 167 CALL FLUSH(numout) 171 168 ENDIF 172 169 … … 183 180 ! vertically averaged velocity 184 181 DO jk = 1, jpk - 1 185 zun(:,:)=zun(:,:) + e1u(:,:)*un(:,:,jk)*fse3u(:,:,jk)186 zvn(:,:)=zvn(:,:) + e2v(:,:)*vn(:,:,jk)*fse3v(:,:,jk)182 zun(:,:)=zun(:,:) + e1u(:,:)*un(:,:,jk)*fse3u(:,:,jk) 183 zvn(:,:)=zvn(:,:) + e2v(:,:)*vn(:,:,jk)*fse3v(:,:,jk) 187 184 END DO 188 185 … … 193 190 !Curl 194 191 DO ji=1,jpim1 195 DO jj=1,jpjm1 196 197 vor_avr(ji,jj)=((zvn(ji+1,jj)-zvn(ji,jj))- & 198 (zun(ji,jj+1)-zun(ji,jj))) & 199 / ( e1f(ji,jj) * e2f(ji,jj) ) 200 vor_avr(ji,jj)=vor_avr(ji,jj)*fmask(ji,jj,1) 201 END DO 192 DO jj=1,jpjm1 193 vor_avr(ji,jj) = ((zvn(ji+1,jj)-zvn(ji,jj))- & 194 (zun(ji,jj+1)-zun(ji,jj))) & 195 /( e1f(ji,jj) * e2f(ji,jj) ) 196 vor_avr(ji,jj) = vor_avr(ji,jj)*fmask(ji,jj,1) 197 END DO 202 198 END DO 203 199 … … 212 208 ! ================================= 213 209 214 # include "trdvor_ncinit.h90"210 # include "trdvor_ncinit.h90" 215 211 216 212 IF( idebug /= 0 ) THEN … … 219 215 ENDIF 220 216 221 ! =====================================222 ! III vertical integration of 3D trends223 ! =====================================224 ! Beta.V : intergration, no average225 utrd(:,:,:,9)=utrd(:,:,:,4)226 vtrd(:,:,:,9)=vtrd(:,:,:,4)227 228 DO jl=1,jplvor229 230 udpvor(:,:)=0231 vdpvor(:,:)=0232 233 !bottom friction234 IF(jl==jplvor) THEN217 ! ===================================== 218 ! III vertical integration of 3D trends 219 ! ===================================== 220 ! Beta.V : intergration, no average 221 utrd(:,:,:,9)=utrd(:,:,:,4) 222 vtrd(:,:,:,9)=vtrd(:,:,:,4) 223 224 DO jl=1,jplvor 225 226 udpvor(:,:)=0 227 vdpvor(:,:)=0 228 229 !bottom friction 230 IF( jl == jplvor ) THEN 235 231 236 call lbc_lnk( tautrd(:,:,3), 'U' , -1. ) 237 call lbc_lnk( tautrd(:,:,4), 'V' , -1. ) 238 239 DO jj = 2, jpjm1 240 DO ji = fs_2, fs_jpim1 241 242 ikbu = min( mbathy(ji+1,jj), mbathy(ji,jj) ) 243 ikbum1 = max( ikbu-1, 1 ) 244 ikbv = min( mbathy(ji,jj+1), mbathy(ji,jj) ) 245 ikbvm1 = max( ikbv-1, 1 ) 246 247 udpvor(ji,jj)=tautrd(ji,jj,3)*fse3u(ji,jj,ikbum1)*e1u(ji,jj)*umask(ji,jj,ikbum1) 248 vdpvor(ji,jj)=tautrd(ji,jj,4)*fse3v(ji,jj,ikbvm1)*e2v(ji,jj)*vmask(ji,jj,ikbvm1) 232 CALL lbc_lnk( tautrd(:,:,3), 'U' , -1. ) 233 CALL lbc_lnk( tautrd(:,:,4), 'V' , -1. ) 234 235 DO jj = 2, jpjm1 236 DO ji = fs_2, fs_jpim1 237 ikbu = min( mbathy(ji+1,jj), mbathy(ji,jj) ) 238 ikbum1 = max( ikbu-1, 1 ) 239 ikbv = min( mbathy(ji,jj+1), mbathy(ji,jj) ) 240 ikbvm1 = max( ikbv-1, 1 ) 241 242 udpvor(ji,jj)=tautrd(ji,jj,3)*fse3u(ji,jj,ikbum1)*e1u(ji,jj)*umask(ji,jj,ikbum1) 243 vdpvor(ji,jj)=tautrd(ji,jj,4)*fse3v(ji,jj,ikbvm1)*e2v(ji,jj)*vmask(ji,jj,ikbvm1) 244 END DO 245 END DO 246 247 !wind stress 248 ELSE IF( jl == (jplvor-1) ) THEN 249 250 CALL lbc_lnk( tautrd(:,:,1), 'U' , -1. ) 251 CALL lbc_lnk( tautrd(:,:,2), 'V' , -1. ) 252 253 udpvor(:,:)=tautrd(:,:,1)*fse3u(:,:,1)*e1u(:,:)*umask(:,:,1) 254 vdpvor(:,:)=tautrd(:,:,2)*fse3v(:,:,1)*e2v(:,:)*vmask(:,:,1) 255 256 ELSE 257 258 CALL lbc_lnk( utrd(:,:,:,jl), 'U' , -1. ) 259 CALL lbc_lnk( vtrd(:,:,:,jl), 'V' , -1. ) 260 261 !utrd and vtrd terms 262 DO jk = 1,jpk 263 udpvor(:,:)=udpvor(:,:)+utrd(:,:,jk,jl)*fse3u(:,:,jk)*e1u(:,:)*umask(:,:,jk) 264 vdpvor(:,:)=vdpvor(:,:)+vtrd(:,:,jk,jl)*fse3v(:,:,jk)*e2v(:,:)*vmask(:,:,jk) 265 END DO 266 267 ENDIF 268 269 !average except for Beta.V 270 IF (jl/=9) THEN 271 udpvor(:,:) = udpvor(:,:) * hur(:,:) 272 vdpvor(:,:) = vdpvor(:,:) * hvr(:,:) 273 ENDIF 274 275 !Curl 276 DO ji=1,jpim1 277 DO jj=1,jpjm1 278 vortrd(ji,jj,jl)=( vdpvor(ji+1,jj)-vdpvor(ji,jj) & 279 - ( udpvor(ji,jj+1)-udpvor(ji,jj) ) ) & 280 / ( e1f(ji,jj) * e2f(ji,jj) ) 281 END DO 249 282 END DO 250 END DO 251 252 !wind stress 253 ELSE IF(jl==(jplvor-1)) THEN 254 255 call lbc_lnk( tautrd(:,:,1), 'U' , -1. ) 256 call lbc_lnk( tautrd(:,:,2), 'V' , -1. ) 257 258 udpvor(:,:)=tautrd(:,:,1)*fse3u(:,:,1)*e1u(:,:)*umask(:,:,1) 259 vdpvor(:,:)=tautrd(:,:,2)*fse3v(:,:,1)*e2v(:,:)*vmask(:,:,1) 260 261 ELSE 262 263 call lbc_lnk( utrd(:,:,:,jl), 'U' , -1. ) 264 call lbc_lnk( vtrd(:,:,:,jl), 'V' , -1. ) 265 266 !utrd and vtrd terms 267 DO jk = 1,jpk 268 udpvor(:,:)=udpvor(:,:)+utrd(:,:,jk,jl)*fse3u(:,:,jk)*e1u(:,:)*umask(:,:,jk) 269 vdpvor(:,:)=vdpvor(:,:)+vtrd(:,:,jk,jl)*fse3v(:,:,jk)*e2v(:,:)*vmask(:,:,jk) 283 284 vortrd(:,:,9)=vortrd(:,:,9)*hur(:,:) 285 286 !surface mask 287 DO ji=1,jpi 288 DO jj=1,jpj 289 vortrd(ji,jj,jl)=vortrd(ji,jj,jl)*fmask(ji,jj,1) !surface mask 290 END DO 291 END DO 292 270 293 END DO 271 272 ENDIF273 274 !average except for Beta.V275 IF (jl/=9) THEN276 udpvor(:,:) = udpvor(:,:) * hur(:,:)277 vdpvor(:,:) = vdpvor(:,:) * hvr(:,:)278 ENDIF279 280 !Curl281 DO ji=1,jpim1282 DO jj=1,jpjm1283 284 vortrd(ji,jj,jl)=(vdpvor(ji+1,jj)-vdpvor(ji,jj) &285 - (udpvor(ji,jj+1)-udpvor(ji,jj))) &286 / ( e1f(ji,jj) * e2f(ji,jj) )287 288 END DO289 END DO290 291 vortrd(:,:,9)=vortrd(:,:,9)*hur(:,:)292 293 !surface mask294 DO ji=1,jpi295 DO jj=1,jpj296 vortrd(ji,jj,jl)=vortrd(ji,jj,jl)*fmask(ji,jj,1) !surface mask297 END DO298 END DO299 300 END DO301 294 302 295 IF( idebug /= 0 ) THEN … … 321 314 ENDIF 322 315 323 324 316 ! IV.2 cumulated trends over analysis period (kt=2 to nwrite) 325 317 ! ---------------------- … … 329 321 nmoydpvor = nmoydpvor + 1 330 322 DO jl = 1, jplvor 331 IF (jl/=9) THEN332 333 323 IF( jl /= 9 ) THEN 324 rotot(:,:) = rotot(:,:) + vortrd(:,:,jl) 325 ENDIF 334 326 END DO 335 327 ENDIF … … 339 331 CALL FLUSH(numout) 340 332 ENDIF 341 342 343 333 344 334 ! ============================================= … … 352 342 353 343 vor_avrtot(:,:) = ( vor_avr(:,:) - vor_avrbn(:,:) + vor_avrb(:,:) - & 354 vor_avrbb(:,:) ) / (zmean * 2. * rdt) 355 356 IF(idebug /= 0) THEN 357 WRITE(numout,*) ' zmean = ',zmean 358 WRITE(numout,*) ' debuging trd_vor: V.1 done' 359 CALL FLUSH(numout) 360 ENDIF 361 344 vor_avrbb(:,:) ) / (zmean * 2. * rdt) 345 346 IF( idebug /= 0 ) THEN 347 WRITE(numout,*) ' zmean = ',zmean 348 WRITE(numout,*) ' debuging trd_vor: V.1 done' 349 CALL FLUSH(numout) 350 ENDIF 362 351 363 352 ! V.2 compute residual 364 353 ! --------------------- 365 366 354 vor_avrres(:,:) = vor_avrtot(:,:) - rotot(:,:) / zmean 367 ! Boundary conditions 368 355 356 ! Boundary conditions 369 357 CALL lbc_lnk( vor_avrtot, 'F', 1. ) 370 358 CALL lbc_lnk( vor_avrres, 'F', 1. ) … … 374 362 CALL FLUSH(numout) 375 363 ENDIF 376 377 364 378 365 ! V.3 time evolution array swap … … 386 373 ENDIF 387 374 388 375 nmoydpvor=0 389 376 390 377 ENDIF … … 394 381 IF( kt >= nit000+1 ) THEN 395 382 396 397 383 #include "trdvor_ncwrite.h90" 398 384 … … 404 390 ENDIF 405 391 406 IF( MOD( kt - nit000+1, ntrd ) == 0 ) THEN407 408 ENDIF392 IF( MOD( kt - nit000+1, ntrd ) == 0 ) THEN 393 rotot(:,:)=0 394 ENDIF 409 395 410 396 IF( kt == nitend ) THEN 411 CALL histclo( nidvor ) 412 !CLOSE(99) 413 ENDIF 397 CALL histclo( nidvor ) 398 ENDIF 414 399 415 400 END SUBROUTINE trd_vor … … 422 407 CONTAINS 423 408 SUBROUTINE trd_vor( kt ) ! Empty routine 424 WRITE(*,*) kt409 WRITE(*,*) 'trd_vor: You should not have seen this print! error?', kt 425 410 END SUBROUTINE trd_vor 426 411 #endif … … 428 413 !!====================================================================== 429 414 END MODULE trdvor 430 431 432 433 434 435 436 437 438 439 -
trunk/NEMO/OPA_SRC/TRD/trdvor_ncinit.h90
r109 r129 2 2 ! II.1 Define frequency of output and means 3 3 ! ----------------------------------------- 4 5 4 #if defined key_diainstant 6 5 zsto = nwrite*rdt … … 13 12 14 13 IF( kt == nit000 ) THEN 15 16 14 IF(lwp) WRITE (numout,*) ' Diavor: netCDF initialization' 17 15 18 16 ! II.2 Compute julian date from starting date of the run 19 17 ! ------------------------ 20 21 18 CALL ymds2ju( nyear, nmonth, nday, 0.e0, zjulian ) 22 19 IF (lwp) WRITE(numout,*)' ' 23 IF (lwp) WRITE(numout,*)' Date 0 used :',nit000 &20 IF (lwp) WRITE(numout,*)' Date 0 used :',nit000 & 24 21 ,' YEAR ', nyear,' MONTH ', nmonth,' DAY ', nday & 25 22 ,'Julian day : ', zjulian 26 23 27 28 24 ! II.3 Define the T grid trend file (nidvor) 29 25 ! --------------------------------- 30 31 26 CALL dia_nam( clhstnam, ntrd, 'vort' ) ! filename 32 27 IF(lwp) WRITE(numout,*) ' Name of NETCDF file ', clhstnam 33 28 CALL histbeg( clhstnam, jpi, glamf, jpj, gphif,1, jpi, & ! Horizontal grid : glamt and gphit 34 1, jpj, 0, zjulian, rdt, nh_t, nidvor)29 & 1, jpj, 0, zjulian, rdt, nh_t, nidvor) 35 30 CALL wheneq( jpi*jpj, fmask, 1, 1., ndexvor1, ndimvor1 ) ! surface 36 37 31 38 32 ! Declare output fields as netCDF variables 39 40 CALL histdef( nidvor, "sovortPh", cvort//"grad Ph" , "s-2", & ! grad Ph 41 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 42 CALL histdef( nidvor, "sovortEk", cvort//"Energy", "s-2", & ! Energy 43 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 33 CALL histdef( nidvor, "sovortPh", cvort//"grad Ph" , "s-2", & ! grad Ph 34 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 35 CALL histdef( nidvor, "sovortEk", cvort//"Energy", "s-2", & ! Energy 36 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 44 37 CALL histdef( nidvor, "sovozeta", cvort//"rel vorticity", "s-2", & ! rel vorticity 45 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)46 CALL histdef( nidvor, "sovortif", cvort//"coriolis", "s-2", & ! coriolis47 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)48 CALL histdef( nidvor, "sovodifl", cvort//"lat diff ", "s-2", & ! lat diff49 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)50 CALL histdef( nidvor, "sovoadvv", cvort//"vert adv", "s-2", & ! vert adv51 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)52 CALL histdef( nidvor, "sovodifv", cvort//"vert diff" , "s-2", & ! vert diff53 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)54 CALL histdef( nidvor, "sovortPs", cvort//"grad Ps", "s-2", & ! grad Ps55 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)56 CALL histdef( nidvor, "sovortbv", cvort//"Beta V", "s-2", & ! beta.V57 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)58 CALL histdef( nidvor, "sovowind", cvort//"wind stress", "s-2", & ! wind stress59 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)60 CALL histdef( nidvor, "sovobfri", cvort//"bottom friction", "s-2", 61 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)62 CALL histdef( nidvor, "1st_mbre", cvort//"1st mbre", "s-2", & !First membre63 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)64 CALL histdef( nidvor, "sovorgap", cvort//"gap", "s-2", & !gap between 1st and 2 nd mbre65 &jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout)38 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 39 CALL histdef( nidvor, "sovortif", cvort//"coriolis", "s-2", & ! coriolis 40 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 41 CALL histdef( nidvor, "sovodifl", cvort//"lat diff ", "s-2", & ! lat diff 42 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 43 CALL histdef( nidvor, "sovoadvv", cvort//"vert adv", "s-2", & ! vert adv 44 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 45 CALL histdef( nidvor, "sovodifv", cvort//"vert diff" , "s-2", & ! vert diff 46 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 47 CALL histdef( nidvor, "sovortPs", cvort//"grad Ps", "s-2", & ! grad Ps 48 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 49 CALL histdef( nidvor, "sovortbv", cvort//"Beta V", "s-2", & ! beta.V 50 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 51 CALL histdef( nidvor, "sovowind", cvort//"wind stress", "s-2", & ! wind stress 52 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 53 CALL histdef( nidvor, "sovobfri", cvort//"bottom friction", "s-2", & ! bottom friction 54 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 55 CALL histdef( nidvor, "1st_mbre", cvort//"1st mbre", "s-2", & ! First membre 56 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 57 CALL histdef( nidvor, "sovorgap", cvort//"gap", "s-2", & ! gap between 1st and 2 nd mbre 58 & jpi,jpj,nh_t,1,1,1,-99,32,clop,zsto,zout) 66 59 CALL histend( nidvor ) 67 60 -
trunk/NEMO/OPA_SRC/TRD/trdvor_ncwrite.h90
r109 r129 2 2 it= kt-nit000+1 3 3 IF( lwp .AND. MOD( kt, ntrd ) == 0 ) THEN 4 5 4 WRITE(numout,*) ' diavor : write NetCDF fields' 6 5 ENDIF 7 8 6 9 CALL histwrite( nidvor,"sovortPh",it,vortrd(:,:,1),ndimvor1,ndexvor1) ! grad Ph10 CALL histwrite( nidvor,"sovortEk",it,vortrd(:,:,2),ndimvor1,ndexvor1) ! Energy11 CALL histwrite( nidvor,"sovozeta",it,vortrd(:,:,3),ndimvor1,ndexvor1) ! rel vorticity12 CALL histwrite( nidvor,"sovortif",it,vortrd(:,:,4),ndimvor1,ndexvor1) ! coriolis13 CALL histwrite( nidvor,"sovodifl",it,vortrd(:,:,5),ndimvor1,ndexvor1) ! lat diff14 CALL histwrite( nidvor,"sovoadvv",it,vortrd(:,:,6),ndimvor1,ndexvor1) ! vert adv15 CALL histwrite( nidvor,"sovodifv",it,vortrd(:,:,7),ndimvor1,ndexvor1) ! vert diff16 CALL histwrite( nidvor,"sovortPs",it,vortrd(:,:,8),ndimvor1,ndexvor1) ! grad Ps17 CALL histwrite( nidvor,"sovortbv",it,vortrd(:,:,9),ndimvor1,ndexvor1) ! beta.V7 CALL histwrite( nidvor,"sovortPh",it,vortrd(:,:,1),ndimvor1,ndexvor1) ! grad Ph 8 CALL histwrite( nidvor,"sovortEk",it,vortrd(:,:,2),ndimvor1,ndexvor1) ! Energy 9 CALL histwrite( nidvor,"sovozeta",it,vortrd(:,:,3),ndimvor1,ndexvor1) ! rel vorticity 10 CALL histwrite( nidvor,"sovortif",it,vortrd(:,:,4),ndimvor1,ndexvor1) ! coriolis 11 CALL histwrite( nidvor,"sovodifl",it,vortrd(:,:,5),ndimvor1,ndexvor1) ! lat diff 12 CALL histwrite( nidvor,"sovoadvv",it,vortrd(:,:,6),ndimvor1,ndexvor1) ! vert adv 13 CALL histwrite( nidvor,"sovodifv",it,vortrd(:,:,7),ndimvor1,ndexvor1) ! vert diff 14 CALL histwrite( nidvor,"sovortPs",it,vortrd(:,:,8),ndimvor1,ndexvor1) ! grad Ps 15 CALL histwrite( nidvor,"sovortbv",it,vortrd(:,:,9),ndimvor1,ndexvor1) ! beta.V 18 16 CALL histwrite( nidvor,"sovowind",it,vortrd(:,:,10),ndimvor1,ndexvor1) ! wind stress 19 17 CALL histwrite( nidvor,"sovobfri",it,vortrd(:,:,11),ndimvor1,ndexvor1) ! bottom friction
Note: See TracChangeset
for help on using the changeset viewer.