New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 5794 for branches/2015/dev_r5187_UKMO13_simplification/NEMOGCM/TOOLS/SCOORD_GEN/src/scoord_gen.F90 – NEMO

Ignore:
Timestamp:
2015-10-14T12:22:34+02:00 (8 years ago)
Author:
timgraham
Message:

Minor bug fixes and added namelist file

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5187_UKMO13_simplification/NEMOGCM/TOOLS/SCOORD_GEN/src/scoord_gen.F90

    r5769 r5794  
    4141      !!          s_tanh  (Madec et al 1996) 
    4242      !!                a cosh/tanh function that gives stretched coordinates         
     43      !! 
     44      !! ** History: 2015: Tim Graham - Code created based on online zdf_sco routine 
     45      !! 
    4346      !! 
    4447      !!---------------------------------------------------------------------- 
     
    201204      ! Envelope bathymetry saved in hbatt 
    202205      hbatt(:,:) = zenv(:,:)  
    203       ! TODO - get this section to work 
    204206      IF( ln_eq_taper) THEN 
    205207        CALL READ_GPHIT() 
    206208        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'  
    208210          DO jj = 1, jpj 
    209211             DO ji = 1, jpi 
     
    212214             END DO 
    213215          END DO 
     216        ENDIF 
    214217      ENDIF 
    215218      ! 
     
    298301      ENDIF 
    299302 
    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() 
    301305      CALL check_nf90( nf90_close(ncout) ) 
    302306 
     
    403407 
    404408         z_gsigt3m1 = z_gsigt3 
    405          z_gsiw3m1 = z_gsiw3 
     409         z_gsi3w3m1 = z_gsi3w3 
    406410  
    407411         where (e3t_0   (:,:).eq.0.0)  e3t_0(:,:) = 1.0 
     
    413417         where (e3vw_0  (:,:).eq.0.0)  e3vw_0(:,:) = 1.0 
    414418         
    415          CALL write_netcdf_vars(jk) 
     419         CALL write_netcdf_3d_vars(jk) 
    416420         DO jj = 1, jpj 
    417421            DO ji = 1, jpi 
     
    507511             IF( jk .EQ. 1) THEN 
    508512                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)  
    510514             ELSE 
    511515                z_esigw3(ji,jj) = z_gsigt3(ji,jj) - z_gsigt3m1(ji,jj) 
     
    552556         ! Keep some arrays for next level 
    553557         z_gsigt3m1 = z_gsigt3 
    554          z_gsiw3m1 = z_gsiw3 
     558         z_gsi3w3m1 = z_gsi3w3 
    555559  
    556560         where (e3t_0   (:,:).eq.0.0)  e3t_0(:,:) = 1.0 
     
    562566         where (e3vw_0  (:,:).eq.0.0)  e3vw_0(:,:) = 1.0 
    563567          
    564          CALL write_netcdf_vars(jk) 
     568         CALL write_netcdf_3d_vars(jk) 
    565569 
    566570         DO jj = 1, jpj 
     
    622626      DO jk = 2, jpk 
    623627         z_gsi3w(jk) = z_gsi3w(jk-1) + z_esigw(jk) 
    624       END DO 
    625 !!gm: depuw, depvw can be suppressed (modif in ldfslp) and depw=dep3w can be set (save 3 3D arrays) 
    626       DO jk = 1, jpk 
    627628      END DO 
    628629!!gm: e3uw, e3vw can be suppressed  (modif in dynzdf, dynzdf_iso, zdfbfr) (save 2 3D arrays) 
     
    655656         where (e3vw_0  (:,:).eq.0.0)  e3vw_0(:,:) = 1.0 
    656657  
    657          CALL write_netcdf_vars(jk) 
     658         CALL write_netcdf_3d_vars(jk) 
    658659      ENDDO ! End of loop over jk 
    659660 
Note: See TracChangeset for help on using the changeset viewer.