Changeset 5794 for branches/2015/dev_r5187_UKMO13_simplification/NEMOGCM/TOOLS/SCOORD_GEN/src/scoord_gen.F90
- Timestamp:
- 2015-10-14T12:22:34+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5187_UKMO13_simplification/NEMOGCM/TOOLS/SCOORD_GEN/src/scoord_gen.F90
r5769 r5794 41 41 !! s_tanh (Madec et al 1996) 42 42 !! a cosh/tanh function that gives stretched coordinates 43 !! 44 !! ** History: 2015: Tim Graham - Code created based on online zdf_sco routine 45 !! 43 46 !! 44 47 !!---------------------------------------------------------------------- … … 201 204 ! Envelope bathymetry saved in hbatt 202 205 hbatt(:,:) = zenv(:,:) 203 ! TODO - get this section to work204 206 IF( ln_eq_taper) THEN 205 207 CALL READ_GPHIT() 206 208 IF( MINVAL( gphit(:,:) ) * MAXVAL( gphit(:,:) ) <= 0. ) THEN 207 CALL ctl_warn( ' s-coordinates are tapered in vicinity of the Equator' )209 WRITE(*,*) 's-coordinates are tapered in vicinity of the Equator' 208 210 DO jj = 1, jpj 209 211 DO ji = 1, jpi … … 212 214 END DO 213 215 END DO 216 ENDIF 214 217 ENDIF 215 218 ! … … 298 301 ENDIF 299 302 300 CALL check_nf90( nf90_put_var(ncout, var_ids(11), mbathy) ) 303 !Write all 2D variables to output file 304 CALL write_netcdf_2d_vars() 301 305 CALL check_nf90( nf90_close(ncout) ) 302 306 … … 403 407 404 408 z_gsigt3m1 = z_gsigt3 405 z_gsi w3m1 = z_gsiw3409 z_gsi3w3m1 = z_gsi3w3 406 410 407 411 where (e3t_0 (:,:).eq.0.0) e3t_0(:,:) = 1.0 … … 413 417 where (e3vw_0 (:,:).eq.0.0) e3vw_0(:,:) = 1.0 414 418 415 CALL write_netcdf_ vars(jk)419 CALL write_netcdf_3d_vars(jk) 416 420 DO jj = 1, jpj 417 421 DO ji = 1, jpi … … 507 511 IF( jk .EQ. 1) THEN 508 512 z_esigw3(ji,jj ) = 2.0 * (z_gsigt3(ji,jj ) - z_gsigw3(ji,jj )) 509 z_gsi3w3(ji,jj) = 0.5 * z_esigw3(ji,jj) 513 z_gsi3w3(ji,jj) = 0.5 * z_esigw3(ji,jj) 510 514 ELSE 511 515 z_esigw3(ji,jj) = z_gsigt3(ji,jj) - z_gsigt3m1(ji,jj) … … 552 556 ! Keep some arrays for next level 553 557 z_gsigt3m1 = z_gsigt3 554 z_gsi w3m1 = z_gsiw3558 z_gsi3w3m1 = z_gsi3w3 555 559 556 560 where (e3t_0 (:,:).eq.0.0) e3t_0(:,:) = 1.0 … … 562 566 where (e3vw_0 (:,:).eq.0.0) e3vw_0(:,:) = 1.0 563 567 564 CALL write_netcdf_ vars(jk)568 CALL write_netcdf_3d_vars(jk) 565 569 566 570 DO jj = 1, jpj … … 622 626 DO jk = 2, jpk 623 627 z_gsi3w(jk) = z_gsi3w(jk-1) + z_esigw(jk) 624 END DO625 !!gm: depuw, depvw can be suppressed (modif in ldfslp) and depw=dep3w can be set (save 3 3D arrays)626 DO jk = 1, jpk627 628 END DO 628 629 !!gm: e3uw, e3vw can be suppressed (modif in dynzdf, dynzdf_iso, zdfbfr) (save 2 3D arrays) … … 655 656 where (e3vw_0 (:,:).eq.0.0) e3vw_0(:,:) = 1.0 656 657 657 CALL write_netcdf_ vars(jk)658 CALL write_netcdf_3d_vars(jk) 658 659 ENDDO ! End of loop over jk 659 660
Note: See TracChangeset
for help on using the changeset viewer.