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 12555 for branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF – NEMO

Ignore:
Timestamp:
2020-03-13T17:33:02+01:00 (4 years ago)
Author:
charris
Message:

Changes from GO6 package branch (GMED ticket 450):

svn merge -r 11035:11101 ​svn+ssh://charris@forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/UKMO/dev_r5518_GO6_package

Location:
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

    r6486 r12555  
    8181      READ  ( numnam_cfg, namdyn_ldf, IOSTAT = ios, ERR = 902 ) 
    8282902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_ldf in configuration namelist', lwp ) 
    83       IF(lwm) WRITE ( numond, namdyn_ldf ) 
     83      IF(lwm .AND. nprint > 2) WRITE ( numond, namdyn_ldf ) 
    8484 
    8585      IF(lwp) THEN                      ! Parameter print 
    … …  
    159159     nkahm_smag = 1 
    160160#endif 
    161  
     161      ! 
     162      IF(lwp .AND. lflush) CALL flush(numout) 
    162163      ! 
    163164   END SUBROUTINE ldf_dyn_init 
    … …  
    210211      ENDIF 
    211212      ! 
     213      IF(lwp .AND. lflush) CALL flush(numout) 
     214      ! 
    212215   END SUBROUTINE ldf_zpf_1d 
    213216 
    … …  
    252255      ENDIF 
    253256      ! 
     257      IF(lwp .AND. lflush) CALL flush(numout) 
     258      ! 
    254259   END SUBROUTINE ldf_zpf_1d_3d 
    255260 
    … …  
    293298      ENDIF 
    294299      ! 
     300      IF(lwp .AND. lflush) CALL flush(numout) 
     301      ! 
    295302   END SUBROUTINE ldf_zpf_3d 
    296303 

  • branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c2d.h90

    r6486 r12555  
    4141      IF(lwp) WRITE(numout,*) 'ldf_dyn_c2d : 2d lateral eddy viscosity coefficient' 
    4242      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     43      IF(lwp .AND. lflush) CALL flush(numout) 
    4344 
    4445      ! harmonic operator (ahm1, ahm2) : ( T- and F- points) (used for laplacian operators 
    … …  
    8586            CALL prihre(ahm2,jpi,jpj,1,jpi,1,1,jpj,1,1.e-3,numout) 
    8687         ENDIF 
     88         IF(lwp .AND. lflush) CALL flush(numout) 
    8789      ENDIF 
    8890 
    … …  
    119121            CALL prihre(ahm4,jpi,jpj,1,jpi,1,1,jpj,1,1.e-3,numout) 
    120122         ENDIF 
     123         IF(lwp .AND. lflush) CALL flush(numout) 
    121124      ENDIF 
    122125      ! 
    … …  
    160163      IF(lwp) WRITE(numout,*) '~~~~~~  --' 
    161164      IF(lwp) WRITE(numout,*) '        orca ocean configuration' 
     165      IF(lwp .AND. lflush) CALL flush(numout) 
    162166 
    163167      IF( cp_cfg == "orca" .AND. cp_cfz == "antarctic" ) THEN 
    … …  
    316320      ENDIF 
    317321      ! 
     322      IF(lwp .AND. lflush) CALL flush(numout) 
     323      ! 
    318324      CALL wrk_dealloc( jpi   , jpj   , icof  ) 
    319325      ! 
    … …  
    358364      IF(lwp) WRITE(numout,*) '~~~~~~  --' 
    359365      IF(lwp) WRITE(numout,*) '        orca_r1 configuration' 
     366      IF(lwp .AND. lflush) CALL flush(numout) 
    360367 
    361368      IF( cp_cfg == "orca" .AND. cp_cfz == "antarctic" ) THEN 
    … …  
    519526      ENDIF 
    520527      ! 
     528      IF(lwp .AND. lflush) CALL flush(numout) 
     529      ! 
    521530      CALL wrk_dealloc( jpi   , jpj   , icof  ) 
    522531      ! 

  • branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90

    r9975 r12555  
    179179         ENDIF 
    180180      ENDIF 
     181      ! 
     182      IF(lwp .AND. lflush) CALL flush(numout) 
    181183      ! 
    182184      CALL wrk_dealloc( jpk, zcoef ) 
    … …  
    432434      ENDIF 
    433435      ! 
     436      IF(lwp .AND. lflush) CALL flush(numout) 
     437      ! 
    434438      CALL wrk_dealloc( jpi   , jpj   , icof  ) 
    435439      CALL wrk_dealloc( jpk   ,         zcoef ) 

  • branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_smag.F90

    r9583 r12555  
    106106 
    107107 
    108       IF(lwp) WRITE(numout,*) 
    109       IF(lwp) WRITE(numout,*) 'ldf_dyn_smag : 3D lateral eddy viscosity coefficient' 
    110       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
    111       
     108         IF(lwp) THEN 
     109            WRITE(numout,*) 
     110            WRITE(numout,*) 'ldf_dyn_smag : 3D lateral eddy viscosity coefficient' 
     111            WRITE(numout,*) '~~~~~~~~~~~' 
     112            IF(lflush) CALL flush(numout) 
     113         ENDIF 
     114 
    112115      ENDIF 
    113116       
    … …  
    172175 
    173176         ENDDO ! jpk 
    174             ahm1(:,:,jpk) = ahm1(:,:,jpkm1) 
    175             ahm2(:,:,jpk) = ahm2(:,:,jpkm1) 
    176             IF(lwp.and.kt==nit000) WRITE(numout,'(36x," ahm ", 7x)') 
     177 
     178         ahm1(:,:,jpk) = ahm1(:,:,jpkm1) 
     179         ahm2(:,:,jpk) = ahm2(:,:,jpkm1) 
     180 
     181         IF(lwp.and.kt==nit000) THEN 
     182            WRITE(numout,'(36x," ahm ", 7x)')  
     183            IF(lflush) CALL flush(numout) 
     184         ENDIF 
     185 
     186         IF(lwp.and.kt==nit000) THEN 
    177187            DO jk = 1, jpk 
    178  
    179                IF(lwp.and.kt==nit000) WRITE(numout,'(30x,E10.2,8x,i3)') ahm1(jpi/2,jpj/2,jk), jk 
     188               WRITE(numout,'(30x,E10.2,8x,i3)') ahm1(jpi/2,jpj/2,jk), jk 
    180189            END DO 
     190            IF(lflush) CALL flush(numout) 
     191         ENDIF 
     192 
    181193      CALL lbc_lnk( ahm1, 'T', 1. )   ! Lateral boundary conditions on ( ahtt ) 
    182194      CALL lbc_lnk( ahm2, 'F', 1. )   ! Lateral boundary conditions on ( ahtt ) 
    … …  
    245257            ahm4(:,:,jpk) = ahm4(:,:,jpkm1) 
    246258 
    247       DO jk = 1, jpk 
    248       IF(  kt == nit000 ) THEN 
    249  
    250                IF(lwp) WRITE(numout,'(30x,E10.2,8x,i3)') ahm3(jpi/2,jpj/2,jk), jk 
     259      IF(  kt == nit000 .AND. lwp) THEN 
     260         DO jk = 1, jpk 
     261               WRITE(numout,'(30x,E10.2,8x,i3)') ahm3(jpi/2,jpj/2,jk), jk 
     262               IF(lflush) CALL flush(numout) 
     263         END DO 
    251264      ENDIF    
    252       END DO 
     265 
    253266      CALL lbc_lnk( ahm3, 'U', 1. )   ! Lateral boundary conditions 
    254267      CALL lbc_lnk( ahm4, 'V', 1. ) 
    … …  
    269282      IF( lk_mpp )   CALL mpp_max( zdeltat )                 ! max over the global domain 
    270283      ! 
    271       IF( MOD( kt, nwrite ) == 1 .AND. lwp )   WRITE(numout,*) ' ==>> time-step= ',kt,'dynlap:  abs(ahm) max: ', zdeltat 
     284      IF( MOD( kt, nwrite ) == 1 .AND. lwp )  THEN 
     285         WRITE(numout,*) ' ==>> time-step= ',kt,'dynlap:  abs(ahm) max: ', zdeltat 
     286         IF(lflush) CALL flush(numout) 
     287      ENDIF 
    272288    ENDIF 
    273289    If(ln_dynldf_bilap)THEN 
    … …  
    282298      IF( lk_mpp )   CALL mpp_max( zdeltat )                 ! max over the global domain 
    283299      ! 
    284       IF( MOD( kt, nwrite ) == 1 .AND. lwp )   WRITE(numout,*) ' ==>> time-step= ',kt,'dyn_bilap abs(ahm) max: ', zdeltat 
     300      IF( MOD( kt, nwrite ) == 1 .AND. lwp )  THEN 
     301         WRITE(numout,*) ' ==>> time-step= ',kt,'dyn_bilap abs(ahm) max: ', zdeltat 
     302         IF(lflush) CALL flush(numout) 
     303      ENDIF 
    285304      ! 
    286305   ENDIF 

  • branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF/ldfslp.F90

    r6498 r12555  
    414414         IF(lwp) THEN  
    415415            WRITE(numout,*) '          Horizontal mixing in s-coordinate: slope = slope of s-surfaces'  
     416            IF(lflush) CALL flush(numout) 
    416417         ENDIF  
    417418 
    … …  
    451452            IF(lwp) WRITE(numout,*) ' max slop: u',SQRT( MAXVAL(uslp*uslp)), ' v ', SQRT(MAXVAL(vslp)),  &  
    452453               &                             ' wi', sqrt(MAXVAL(wslpi)), ' wj', sqrt(MAXVAL(wslpj))  
     454            IF(lwp .AND. lflush) CALL flush(numout) 
    453455         endif  
    454456   
    … …  
    842844         WRITE(numout,*) 'ldf_slp_init : direction of lateral mixing' 
    843845         WRITE(numout,*) '~~~~~~~~~~~~' 
     846         IF(lflush) CALL flush(numout) 
    844847      ENDIF 
    845848 
    … …  
    865868         IF(ln_sco .AND.  (ln_traldf_hor .OR. ln_dynldf_hor )) THEN 
    866869            IF(lwp)   WRITE(numout,*) '          Horizontal mixing in s-coordinate: slope = slope of s-surfaces' 
     870            IF(lwp .AND. lflush) CALL flush(numout) 
    867871 
    868872            ! geopotential diffusion in s-coordinates on tracers and/or momentum 

  • branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90

    r6486 r12555  
    8585      READ  ( numnam_cfg, namtra_ldf, IOSTAT = ios, ERR = 902 ) 
    8686902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_ldf in configuration namelist', lwp ) 
    87       IF(lwm) WRITE ( numond, namtra_ldf ) 
     87      IF(lwm .AND. nprint > 2) WRITE ( numond, namtra_ldf ) 
    8888 
    8989      IF(lwp) THEN                      ! control print 
    … …  
    175175        IF(lwp.and.ln_traldf_bilap) CALL ctl_stop(' SMAGORINSKY + BILAPLACIAN - UNSTABLE OR NON_CONSERVATIVE' ) 
    176176#endif 
    177  
     177      ! 
     178      IF(lwp .AND. lflush) CALL flush(numout) 
    178179      ! 
    179180   END SUBROUTINE ldf_tra_init 

  • branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_c2d.h90

    r6486 r12555  
    3636         IF(lwp) WRITE(numout,*) ' ldf_tra_c2d : 2D eddy diffusivity and eddy' 
    3737         IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~   --  induced velocity coefficients' 
     38         IF(lwp .AND. lflush) CALL flush(numout) 
    3839      ELSE 
    3940         IF(lwp) WRITE(numout,*) 
    4041         IF(lwp) WRITE(numout,*) ' ldf_tra2d : 2D eddy diffusivity coefficient' 
    4142         IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~   --' 
     43         IF(lwp .AND. lflush) CALL flush(numout) 
    4244      ENDIF 
    4345 
    … …  
    7577            IF(lwp) WRITE(numout,*) '               Constant values used for eddy diffusivity coefficients' 
    7678            IF(lwp) WRITE(numout,*) '               Variation lat/lon only for eddy induced velocity coefficients' 
     79            IF(lwp .AND. lflush) CALL flush(numout) 
    7780         ENDIF 
    7881 

  • branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_c3d.h90

    r6486 r12555  
    4646         IF(lwp) WRITE(numout,*) '               Coefficients set to constant' 
    4747      ENDIF 
    48  
     48    
    4949 
    5050      ! biharmonic operator   (T-point) 
    … …  
    121121 
    122122# endif 
     123 
     124   IF(lwp .AND. lflush) CALL flush(numout) 
     125 
    123126   END SUBROUTINE ldf_tra_c3d 

  • branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_smag.F90

    r9583 r12555  
    9999         IF(lwp) WRITE(numout,*) 
    100100         IF(lwp) WRITE(numout,*) 
     101         IF(lwp .AND. lflush) CALL flush(numout) 
    101102      ENDIF 
    102103 
    … …  
    108109      ahtt(:,:,:) = rn_aht_0 
    109110       IF( ln_traldf_bilap ) THEN 
    110         IF( lwp .AND. kt == nit000) WRITE(numout,* )'ldf_tra_smag :no bilaplacian Smagorinsky diffusivity' 
    111         IF( lwp .AND. kt == nit000) WRITE(numout,* )'ldf_tra_smag :bilaplacian diffusivity set to constant'   
     111        IF( lwp .AND. kt == nit000) THEN 
     112          WRITE(numout,* )'ldf_tra_smag :no bilaplacian Smagorinsky diffusivity' 
     113          WRITE(numout,* )'ldf_tra_smag :bilaplacian diffusivity set to constant'   
     114          IF(lflush) CALL flush(numout) 
     115        ENDIF 
    112116       ENDIF 
    113117 
    … …  
    190194         CALL prihre( ahtw(:,:,1), jpi, jpj, 1, jpi, 1,   & 
    191195            &                                1, jpj, 1, 1.e-1, numout ) 
     196         IF(lflush) CALL flush(numout) 
    192197      ENDIF 
    193198ENDIF 
Note: See TracChangeset for help on using the changeset viewer.