Changeset 5461


Ignore:
Timestamp:
2015-06-22T14:40:22+02:00 (5 years ago)
Author:
jchanut
Message:

Move output of wn and sshn in diawri, see ticket #1543

Location:
trunk/NEMOGCM/NEMO/OPA_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r5460 r5461  
    149149         CALL iom_put( "e3w" , fse3w_n(:,:,:) ) 
    150150      ENDIF 
     151 
     152      CALL iom_put( "ssh" , sshn )                 ! sea surface height 
     153      if( iom_use('ssh2') )   CALL iom_put( "ssh2", sshn(:,:) * sshn(:,:) )   ! square of sea surface height 
    151154       
    152155      CALL iom_put( "toce", tsn(:,:,:,jp_tem) )    ! 3D temperature 
     
    196199         END DO 
    197200         CALL iom_put( "sbv", z2d )                ! bottom j-current 
     201      ENDIF 
     202 
     203      CALL iom_put( "woce", wn )                   ! vertical velocity 
     204      IF( iom_use('w_masstr') .OR. iom_use('w_masstr2') ) THEN   ! vertical mass transport & its square value 
     205         ! Caution: in the VVL case, it only correponds to the baroclinic mass transport. 
     206         z2d(:,:) = rau0 * e12t(:,:) 
     207         DO jk = 1, jpk 
     208            z3d(:,:,jk) = wn(:,:,jk) * z2d(:,:) 
     209         END DO 
     210         CALL iom_put( "w_masstr" , z3d )   
     211         IF( iom_use('w_masstr2') )   CALL iom_put( "w_masstr2", z3d(:,:,:) * z3d(:,:,:) ) 
    198212      ENDIF 
    199213 
  • trunk/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

    r4990 r5461  
    2121   USE domvvl          ! Variable volume 
    2222   USE divcur          ! hor. divergence and curl      (div & cur routines) 
    23    USE iom             ! I/O library 
    2423   USE restart         ! only for lrst_oce 
    2524   USE in_out_manager  ! I/O manager 
     
    3130   USE bdy_par          
    3231   USE bdydyn2d        ! bdy_ssh routine 
    33    USE iom 
    3432#if defined key_agrif 
    3533   USE agrif_opa_update 
     
    137135      !                                           !           outputs            ! 
    138136      !                                           !------------------------------! 
    139       CALL iom_put( "ssh" , sshn )   ! sea surface height 
    140       if( iom_use('ssh2') )   CALL iom_put( "ssh2", sshn(:,:) * sshn(:,:) )   ! square of sea surface height 
    141137      ! 
    142138      IF(ln_ctl)   CALL prt_ctl( tab2d_1=ssha, clinfo1=' ssha  - : ', mask1=tmask, ovlap=1 ) 
     
    228224#endif 
    229225      ! 
    230       !                                           !------------------------------! 
    231       !                                           !           outputs            ! 
    232       !                                           !------------------------------! 
    233       CALL iom_put( "woce", wn )   ! vertical velocity 
    234       IF( iom_use('w_masstr') .OR. iom_use('w_masstr2') ) THEN   ! vertical mass transport & its square value 
    235          CALL wrk_alloc( jpi, jpj, z2d )  
    236          CALL wrk_alloc( jpi, jpj, jpk, z3d )  
    237          ! Caution: in the VVL case, it only correponds to the baroclinic mass transport. 
    238          z2d(:,:) = rau0 * e12t(:,:) 
    239          DO jk = 1, jpk 
    240             z3d(:,:,jk) = wn(:,:,jk) * z2d(:,:) 
    241          END DO 
    242          CALL iom_put( "w_masstr" , z3d )   
    243          IF( iom_use('w_masstr2') )   CALL iom_put( "w_masstr2", z3d(:,:,:) * z3d(:,:,:) ) 
    244          CALL wrk_dealloc( jpi, jpj, z2d  )  
    245          CALL wrk_dealloc( jpi, jpj, jpk, z3d )  
    246       ENDIF 
    247       ! 
    248226      IF( nn_timing == 1 )  CALL timing_stop('wzv') 
    249227 
Note: See TracChangeset for help on using the changeset viewer.