Changeset 2431


Ignore:
Timestamp:
2010-11-25T13:45:32+01:00 (10 years ago)
Author:
cetlod
Message:

Improve the Offline together with the 1D vertical configuration

Location:
branches/nemo_v3_3_beta/NEMOGCM/NEMO
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OFF_SRC/domain.F90

    r2287 r2431  
    222222      ENDIF 
    223223 
    224       n_cla = nn_cla                ! conversion DOCTOR names into model names (this should disappear soon) 
    225  
    226       IF( lk_mpp_rep .AND. n_cla /= 0 )   CALL ctl_stop( ' Reproductibility tests (nbit_cmp=1) require n_cla = 0' ) 
     224      IF( lk_mpp_rep .AND. nn_cla /= 0 )   CALL ctl_stop( ' Reproductibility tests (nbit_cmp=1) require n_cla = 0' ) 
    227225      ! 
    228226   END SUBROUTINE dom_nam 
     
    268266      IF( ln_sco ) ioptio = ioptio + 1 
    269267      IF ( ioptio /= 1 )   CALL ctl_stop( ' none or several vertical coordinate options used' ) 
    270  
    271       IF( nprint == 1 .AND. lwp )   THEN 
    272          WRITE(numout,*) ' MIN val mbathy ', MINVAL( mbathy(:,:) ), ' MAX ', MAXVAL( mbathy(:,:) ) 
    273          WRITE(numout,*) ' MIN val depth t ', MINVAL( fsdept(:,:,:) ),   & 
    274             &                   ' w ',   MINVAL( fsdepw(:,:,:) ), '3w ', MINVAL( fsde3w(:,:,:) ) 
    275          WRITE(numout,*) ' MIN val e3    t ', MINVAL( fse3t(:,:,:) ), ' f ', MINVAL( fse3f(:,:,:) ),  & 
    276             &                   ' u ',   MINVAL( fse3u(:,:,:) ), ' u ', MINVAL( fse3v(:,:,:) ),  & 
    277             &                   ' w ',   MINVAL( fse3w(:,:,:) ) 
    278  
    279          WRITE(numout,*) ' MAX val depth t ', MAXVAL( fsdept(:,:,:) ),   & 
    280             &                   ' w ',   MAXVAL( fsdepw(:,:,:) ), '3w ', MAXVAL( fsde3w(:,:,:) ) 
    281          WRITE(numout,*) ' MAX val e3    t ', MAXVAL( fse3t(:,:,:) ), ' f ', MAXVAL( fse3f(:,:,:) ),  & 
    282             &                   ' u ',   MAXVAL( fse3u(:,:,:) ), ' u ', MAXVAL( fse3v(:,:,:) ),  & 
    283             &                   ' w ',   MAXVAL( fse3w(:,:,:) ) 
    284       ENDIF 
    285268 
    286269   END SUBROUTINE dom_zgr 
  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OFF_SRC/domrea.F90

    r2287 r2431  
    8686         inum4                      ! 'mesh_zgr.nc'  file 
    8787  
    88       REAL(wp), DIMENSION(jpi,jpj) :: & 
    89          zprt 
    90       INTEGER :: ik 
     88      REAL(wp), DIMENSION(jpi,jpj) :: zprt 
     89      REAL(wp) ::   zrefdep         ! depth of the reference level (~10m) 
     90      INTEGER  :: ik 
    9191      !!---------------------------------------------------------------------- 
    9292 
     
    260260         ENDIF 
    261261 
     262!!gm BUG in s-coordinate this does not work! 
     263      ! deepest/shallowest W level Above/Below ~10m 
     264      zrefdep = 10. - ( 0.1*MINVAL(e3w_0) )                          ! ref. depth with tolerance (10% of minimum layer thickness) 
     265      nlb10 = MINLOC( gdepw_0, mask = gdepw_0 > zrefdep, dim = 1 )   ! shallowest W level Below ~10m 
     266      nla10 = nlb10 - 1                                              ! deepest    W level Above ~10m 
     267!!gm end bug 
    262268 
    263269      ! Control printing : Grid informations (if not restart) 
  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OFF_SRC/opa.F90

    r2287 r2431  
    3838   USE dtadyn          ! Lecture and interpolation of the dynamical fields 
    3939   USE stpctl          ! time stepping control            (stp_ctl routine) 
     40 
     41   USE c1d             ! 1D configuration 
    4042 
    4143   USE iom 
     
    335337         END SELECT 
    336338      ENDIF 
     339      ! 
     340      IF( lk_c1d .AND. .NOT. lk_iomput )  & 
     341        CALL ctl_stop( ' The 1D vertical configuration must be used in conjunction',   & 
     342            &          ' with the IOM Input/Output manager. Compile with key_iomput enabled' ) 
     343      ! 
    337344 
    338345   END SUBROUTINE opa_flg 
  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/par_ORCA_R2.h90

    r2379 r2431  
    8282      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom 
    8383      ! jpjzoom =   27   ,    &  !: in data domain indices (160W,60S) 
    84       jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom 
    85       jpjzoom =    7   ,    &  !: in data domain indices (160W,75S) 
     84      ! jpizoom =   61   ,    &  !: left bottom (i,j) indices of the zoom 
     85      ! jpjzoom =    7   ,    &  !: in data domain indices (160W,75S) 
     86      jpizoom =   110   ,    &  !: left bottom (i,j) indices of the zoom 
     87      jpjzoom =   97   ,    &  !: in data domain indices (64W,31.5N) BATS site 
    8688      ! Domain characteristics 
    8789      jperio  =   0            !: lateral cond. type (between 0 and 6) 
  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/PISCES/trcrst_pisces.F90

    r2403 r2431  
    119119      IF(lwp)  WRITE(numout,*) 
    120120 
    121       IF( cp_cfg == "orca" .AND. .NOT. lk_trc_c1d ) THEN      ! ORCA condiguration (not 1D) ! 
     121      IF( cp_cfg == "orca" .AND. .NOT. lk_c1d ) THEN      ! ORCA condiguration (not 1D) ! 
    122122         !                                                    ! --------------------------- ! 
    123123         ! set total alkalinity, phosphate, nitrate & silicate 
  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/TRP/trctrp.F90

    r2287 r2431  
    6060      INTEGER, INTENT( in ) ::  kstp  ! ocean time-step index 
    6161      !! --------------------------------------------------------------------- 
    62       IF( .NOT. lk_trc_c1d ) THEN 
     62      IF( .NOT. lk_c1d ) THEN 
    6363         ! 
    6464                                CALL trc_sbc( kstp )            ! surface boundary condition 
  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/oce_trc.F90

    r2402 r2431  
    3232 
    3333   !* IO manager * 
    34    USE in_out_manager      ! use all the variables 
     34   USE in_out_manager     
     35                           
    3536   !* physical constants * 
    36    USE phycst              ! use all the variables 
     37   USE phycst             
     38 
     39   !* 1D configuration 
     40   USE c1d                                          
    3741 
    3842   !* model domain * 
  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/trcini.F90

    r2402 r2431  
    138138      tra(:,:,:,:) = 0. 
    139139       
    140       IF( ln_zps .AND. .NOT. lk_trc_c1d )   &              ! Partial steps: before horizontal gradient of passive 
     140      IF( ln_zps .AND. .NOT. lk_c1d )   &              ! Partial steps: before horizontal gradient of passive 
    141141      &                     CALL zps_hde( nit000, jptra, trn, gtru, gtrv )       ! tracers at the bottom ocean level 
    142142 
Note: See TracChangeset for help on using the changeset viewer.