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 13759 – NEMO

Changeset 13759


Ignore:
Timestamp:
2020-11-09T18:04:19+01:00 (4 years ago)
Author:
emanuelaclementi
Message:

minor changes to include writing statements and update n. of max coupling fields with oasis - tickets #2155 #2339

Location:
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/SBC/cpl_oasis3.F90

    r13710 r13759  
    6666   INTEGER                    ::   nsnd         ! total number of fields sent  
    6767   INTEGER                    ::   ncplmodel    ! Maximum number of models to/from which NEMO is potentialy sending/receiving data 
    68    INTEGER, PUBLIC, PARAMETER ::   nmaxfld=60   ! Maximum number of coupling fields 
     68   INTEGER, PUBLIC, PARAMETER ::   nmaxfld=62   ! Maximum number of coupling fields 
    6969   INTEGER, PUBLIC, PARAMETER ::   nmaxcat=5    ! Maximum number of coupling fields 
    7070   INTEGER, PUBLIC, PARAMETER ::   nmaxcpl=5    ! Maximum number of coupling fields 
  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/SBC/sbccpl.F90

    r13710 r13759  
    358358         WRITE(numout,*)'                      - mesh          = ', sn_snd_crtw%clvgrd  
    359359      ENDIF 
    360  
     360      IF( lwp .AND. ln_wave) THEN                        ! control print 
     361      WRITE(numout,*)'      surface waves:' 
     362         WRITE(numout,*)'      Significant wave heigth         = ', TRIM(sn_rcv_hsig%cldes  ), ' (', TRIM(sn_rcv_hsig%clcat  ), ')' 
     363         WRITE(numout,*)'      Wave to oce energy flux         = ', TRIM(sn_rcv_phioc%cldes ), ' (', TRIM(sn_rcv_phioc%clcat ), ')' 
     364         WRITE(numout,*)'      Surface Stokes drift grid u     = ', TRIM(sn_rcv_sdrfx%cldes ), ' (', TRIM(sn_rcv_sdrfx%clcat ), ')' 
     365         WRITE(numout,*)'      Surface Stokes drift grid v     = ', TRIM(sn_rcv_sdrfy%cldes ), ' (', TRIM(sn_rcv_sdrfy%clcat ), ')' 
     366         WRITE(numout,*)'      Mean wave period                = ', TRIM(sn_rcv_wper%cldes  ), ' (', TRIM(sn_rcv_wper%clcat  ), ')' 
     367         WRITE(numout,*)'      Mean wave number                = ', TRIM(sn_rcv_wnum%cldes  ), ' (', TRIM(sn_rcv_wnum%clcat  ), ')' 
     368         WRITE(numout,*)'      Stress frac adsorbed by waves   = ', TRIM(sn_rcv_wstrf%cldes ), ' (', TRIM(sn_rcv_wstrf%clcat ), ')' 
     369         WRITE(numout,*)'      Neutral surf drag coefficient   = ', TRIM(sn_rcv_wdrag%cldes ), ' (', TRIM(sn_rcv_wdrag%clcat ), ')' 
     370         WRITE(numout,*)'      Charnock coefficient            = ', TRIM(sn_rcv_charn%cldes ), ' (', TRIM(sn_rcv_charn%clcat ), ')' 
     371         WRITE(numout,*)' Transport associated to Stokes drift grid u = ', TRIM(sn_rcv_tusd%cldes ), ' (', TRIM(sn_rcv_tusd%clcat ), ')' 
     372         WRITE(numout,*)' Transport associated to Stokes drift grid v = ', TRIM(sn_rcv_tvsd%cldes ), ' (', TRIM(sn_rcv_tvsd%clcat ), ')' 
     373         WRITE(numout,*)'      Bernouilli pressure head        = ', TRIM(sn_rcv_bhd%cldes   ), ' (', TRIM(sn_rcv_bhd%clcat  ), ')' 
     374         WRITE(numout,*)'Wave to ocean momentum flux and Net wave-supported stress = ', TRIM(sn_rcv_taw%cldes ), ' (', TRIM(sn_rcv_taw%clcat ), ')' 
     375         WRITE(numout,*)'      Surface current to waves        = ', TRIM(sn_snd_crtw%cldes  ), ' (', TRIM(sn_snd_crtw%clcat  ), ')' 
     376         WRITE(numout,*)'                      - referential   = ', sn_snd_crtw%clvref 
     377         WRITE(numout,*)'                      - orientation   = ', sn_snd_crtw%clvor 
     378         WRITE(numout,*)'                      - mesh          = ', sn_snd_crtw%clvgrd 
     379      ENDIF 
    361380      !                                   ! allocate sbccpl arrays 
    362381      IF( sbc_cpl_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'sbc_cpl_alloc : unable to allocate arrays' ) 
     
    10821101      !   initialisation of the coupler  ! 
    10831102      ! ================================ ! 
    1084  
    10851103      CALL cpl_define(jprcv, jpsnd, nn_cplmodel) 
    10861104       
  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/SBC/sbcwave.F90

    r13710 r13759  
    353353 
    354354      READ  ( numnam_cfg, namsbc_wave, IOSTAT = ios, ERR = 902 ) 
    355 902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namsbc_wave in configurationnamelist' ) 
     355902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namsbc_wave in configuration namelist' ) 
    356356      IF(lwm) WRITE ( numond, namsbc_wave ) 
    357357      ! 
    358358      IF(lwp) THEN 
    359359         WRITE(numout,*) '   Namelist namsbc_wave' 
    360          WRITE(numout,*) '      Stokes drift                               ln_sdw       = ', ln_sdw 
    361          WRITE(numout,*) '      Stokes Coriolis & tracer advection terms   ln_stcor     = ', ln_stcor 
    362          WRITE(numout,*) '      wave modified ocean stress                 ln_tauoc     = ', ln_tauoc 
    363          WRITE(numout,*) '      neutral drag coefficient (CORE bulk only)  ln_cdgw      = ', ln_cdgw 
    364          WRITE(numout,*) '      charnock coefficient (IFS bulk only)       ln_charn     = ', ln_charn 
    365          WRITE(numout,*) '      Test with constant wave fields             ln_wave_test = ', ln_wave_test 
     360         WRITE(numout,*) '      Stokes drift                                  ln_sdw = ', ln_sdw 
     361         WRITE(numout,*) '      Breivik 2016                       ln_breivikFV_2016 = ', ln_breivikFV_2016 
     362         WRITE(numout,*) '      Stokes Coriolis & tracer advection terms    ln_stcor = ', ln_stcor 
     363         WRITE(numout,*) '      Vortex Force                         ln_vortex_force = ', ln_vortex_force 
     364         WRITE(numout,*) '      Bernouilli Head Pressure                ln_bern_srfc = ', ln_bern_srfc 
     365         WRITE(numout,*) '      wave modified ocean stress                  ln_tauoc = ', ln_tauoc 
     366         WRITE(numout,*) '      neutral drag coefficient (CORE bulk only)    ln_cdgw = ', ln_cdgw 
     367         WRITE(numout,*) '      charnock coefficient                        ln_charn = ', ln_charn 
     368         WRITE(numout,*) '      Stress modificated by wave                    ln_taw = ', ln_taw 
     369         WRITE(numout,*) '      TKE flux from wave                          ln_phioc = ', ln_phioc 
     370         WRITE(numout,*) '      Surface shear with Stokes drift           ln_stshear = ', ln_stshear 
     371         WRITE(numout,*) '      Test with constant wave fields          ln_wave_test = ', ln_wave_test 
    366372      ENDIF 
    367373 
  • NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/ZDF/zdftke.F90

    r13710 r13759  
    408408 
    409409      IF ( cpl_phioc .and. ln_phioc )  THEN 
    410          SELECT CASE (nn_bc_surf)  ! Dirichlet Boundary Condition using surface TKE flux from waves 
    411  
    412          CASE ( 0 ) 
     410         SELECT CASE (nn_bc_surf)  ! Boundary Condition using surface TKE flux from waves 
     411 
     412         CASE ( 0 ) ! Dirichlet BC 
    413413 
    414414            DO_2D( 0, 0, 0, 0 )    ! en(1)   = rn_ebb taum / rho0  (min value rn_emin0) 
     
    420420            END_2D 
    421421 
    422          CASE ( 1 ) 
     422         CASE ( 1 ) ! Neumann BC 
    423423            DO_2D( 0, 0, 0, 0 ) 
    424424               IF ( phioc(ji,jj) < 0 )  phioc(ji,jj) = 0._wp 
     
    769769         WRITE(numout,*) '      Langmuir cells parametrization              ln_lc     = ', ln_lc 
    770770         WRITE(numout,*) '         coef to compute vertical velocity of LC     rn_lc  = ', rn_lc 
     771         IF ( cpl_phioc .and. ln_phioc )  THEN 
     772            SELECT CASE( nn_bc_surf)             ! Type of scaling under sea-ice 
     773            CASE( 0 )   ;   WRITE(numout,*) '  nn_bc_surf=0  ==>>>   DIRICHLET Surface Boundary Condition using surface TKE flux from waves' 
     774            CASE( 1 )   ;   WRITE(numout,*) '  nn_bc_surf=1  ==>>>   NEUMANN Surface Boundary Condition using surface TKE flux from waves' 
     775            END SELECT 
     776         ENDIF  
    771777         WRITE(numout,*) '      test param. to add tke induced by wind      nn_etau   = ', nn_etau 
    772778         WRITE(numout,*) '          type of tke penetration profile            nn_htau   = ', nn_htau 
Note: See TracChangeset for help on using the changeset viewer.