Changeset 12756


Ignore:
Timestamp:
2020-04-16T15:34:34+02:00 (8 months ago)
Author:
clem
Message:

fix xml again

Location:
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/tests/CANAL
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/tests/CANAL/EXPREF/file_def_nemo-oce.xml

    r9572 r12756  
    1515     <field field_ref="soce" />  
    1616     <field field_ref="ssh"  /> 
    17      <field field_ref="salgrad"  /> 
    18      <field field_ref="ke_zint"  /> 
     17     <field field_ref="socegrad"  /> 
     18     <field field_ref="eken_int"  /> 
    1919     <field field_ref="relvor"  /> 
    2020     <field field_ref="potvor"  /> 
  • NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/tests/CANAL/MY_SRC/diawri.F90

    r12206 r12756  
    230230      IF( iom_use('logavs') )   CALL iom_put( "logavs", LOG( MAX( 1.e-20_wp, avs(:,:,:) ) ) ) 
    231231 
    232       IF ( iom_use("salgrad") .OR. iom_use("salgrad2") ) THEN 
     232      IF ( iom_use("socegrad") .OR. iom_use("socegrad2") ) THEN 
    233233         z3d(:,:,jpk) = 0. 
    234234         DO jk = 1, jpkm1 
     
    244244         END DO 
    245245         CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 
    246          CALL iom_put( "salgrad2",  z3d )          ! square of module of sal gradient 
     246         CALL iom_put( "socegrad2",  z3d )          ! square of module of sal gradient 
    247247         z3d(:,:,:) = SQRT( z3d(:,:,:) ) 
    248          CALL iom_put( "salgrad" ,  z3d )          ! module of sal gradient 
     248         CALL iom_put( "socegrad" ,  z3d )          ! module of sal gradient 
    249249      ENDIF 
    250250          
     
    299299            END DO 
    300300         END DO 
    301          CALL iom_put( "salt2c", rau0 * z2d )          ! vertically integrated salt content (PSU*kg/m2) 
    302       ENDIF 
    303       ! 
    304       IF ( iom_use("eken") ) THEN 
     301         CALL iom_put( "salt2c", rau0 * z2d )          ! vertically integrated squared salt content (PSU*kg/m2) 
     302      ENDIF 
     303      ! 
     304      IF ( iom_use("eken") .OR. iom_use("eken_int") ) THEN 
    305305         z3d(:,:,jpk) = 0._wp  
    306306         DO jk = 1, jpkm1 
    307             DO jj = 2, jpj 
    308                DO ji = 2, jpi 
     307            DO jj = 2, jpjm1 
     308               DO ji = 2, jpim1 
    309309                  zztmpx = 0.5 * ( un(ji-1,jj  ,jk) + un(ji,jj,jk) ) 
    310310                  zztmpy = 0.5 * ( vn(ji  ,jj-1,jk) + vn(ji,jj,jk) ) 
     
    315315         CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 
    316316         CALL iom_put( "eken", z3d )                 ! kinetic energy 
    317       ENDIF 
    318  
    319       IF ( iom_use("ke") .or. iom_use("ke_zint") ) THEN 
    320          ! 
    321          z3d(:,:,jpk) = 0._wp 
    322          z3d(1,:, : ) = 0._wp 
    323          z3d(:,1, : ) = 0._wp 
    324          DO jk = 1, jpkm1 
    325             DO jj = 2, jpj 
    326                DO ji = 2, jpi 
    327                   z3d(ji,jj,jk) = 0.25_wp * ( un(ji  ,jj,jk) * un(ji  ,jj,jk) * e1e2u(ji  ,jj) * e3u_n(ji  ,jj,jk)  & 
    328                      &                      + un(ji-1,jj,jk) * un(ji-1,jj,jk) * e1e2u(ji-1,jj) * e3u_n(ji-1,jj,jk)  & 
    329                      &                      + vn(ji,jj  ,jk) * vn(ji,jj  ,jk) * e1e2v(ji,jj  ) * e3v_n(ji,jj  ,jk)  & 
    330                      &                      + vn(ji,jj-1,jk) * vn(ji,jj-1,jk) * e1e2v(ji,jj-1) * e3v_n(ji,jj-1,jk)  )  & 
    331                      &                    * r1_e1e2t(ji,jj) / e3t_n(ji,jj,jk) * tmask(ji,jj,jk) 
    332                END DO 
    333             END DO 
    334          END DO 
    335           
    336          CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 
    337          CALL iom_put( "ke", z3d ) ! kinetic energy 
    338317 
    339318         z2d(:,:)  = 0._wp  
     
    341320            DO jj = 1, jpj 
    342321               DO ji = 1, jpi 
    343                   z2d(ji,jj) = z2d(ji,jj) + e3t_n(ji,jj,jk) * z3d(ji,jj,jk) * tmask(ji,jj,jk) 
    344                END DO 
    345             END DO 
    346          END DO 
    347          CALL iom_put( "ke_zint", z2d )   ! vertically integrated kinetic energy 
    348  
     322                  z2d(ji,jj) = z2d(ji,jj) + e3t_n(ji,jj,jk) * z3d(ji,jj,jk) * e1e2t(ji,jj) * tmask(ji,jj,jk) 
     323               END DO 
     324            END DO 
     325         END DO 
     326         CALL iom_put( "eken_int", z2d )   ! vertically integrated kinetic energy 
    349327      ENDIF 
    350328      ! 
     
    358336               DO ji = 1, fs_jpim1   ! vector opt. 
    359337                  z3d(ji,jj,jk) = (  e2v(ji+1,jj  ) * vn(ji+1,jj  ,jk) - e2v(ji,jj) * vn(ji,jj,jk)    & 
    360                      &              - e1u(ji  ,jj+1) * un(ji  ,jj+1,jk) + e1u(ji,jj) * un(ji,jj,jk)  ) * r1_e1e2f(ji,jj) 
     338                     &             - e1u(ji  ,jj+1) * un(ji  ,jj+1,jk) + e1u(ji,jj) * un(ji,jj,jk)  ) * r1_e1e2f(ji,jj) 
    361339               END DO 
    362340            END DO 
Note: See TracChangeset for help on using the changeset viewer.