Ignore:
Timestamp:
2012-03-12T15:44:43+01:00 (9 years ago)
Author:
gm
Message:

Ediag branche: #927 add Kinetic Energy trend diagnostics (trdken.F90)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdvor.F90

    r3318 r3325  
    116116         CALL trd_vor_zint( putrd, pvtrd, jpvor_zdf )                             ! zdf trend including surf./bot. stresses  
    117117         CALL trd_vor_zint( ztswu, ztswv, jpvor_swf )                             ! surface wind stress  
    118       CASE ( jpdyn_bfr ) 
     118      CASE( jpdyn_bfr ) 
    119119         CALL trd_vor_zint( putrd, pvtrd, jpvor_bfr )                             ! Bottom stress 
    120120         ! 
     
    176176      !  ===================================== 
    177177 
    178       SELECT CASE (ktrd)  
    179       ! 
    180       CASE (jpvor_bfr)        ! bottom friction 
     178      SELECT CASE( ktrd )  
     179      ! 
     180      CASE( jpvor_bfr )        ! bottom friction 
    181181         DO jj = 2, jpjm1 
    182182            DO ji = fs_2, fs_jpim1  
     
    188188         END DO 
    189189         ! 
    190       CASE (jpvor_swf)        ! wind stress 
     190      CASE( jpvor_swf )        ! wind stress 
    191191         zudpvor(:,:) = putrdvor(:,:) * fse3u(:,:,1) * e1u(:,:) * umask(:,:,1) 
    192192         zvdpvor(:,:) = pvtrdvor(:,:) * fse3v(:,:,1) * e2v(:,:) * vmask(:,:,1) 
     
    199199    
    200200      ! Curl 
    201       DO ji=1,jpim1 
    202          DO jj=1,jpjm1 
     201      DO ji = 1, jpim1 
     202         DO jj = 1, jpjm1 
    203203            vortrd(ji,jj,ktrd) = (    zvdpvor(ji+1,jj) - zvdpvor(ji,jj)       & 
    204204                 &                - ( zudpvor(ji,jj+1) - zudpvor(ji,jj) )   ) / ( e1f(ji,jj) * e2f(ji,jj) ) 
Note: See TracChangeset for help on using the changeset viewer.