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 710 for trunk/NEMO/OPA_SRC/DIA/diawri_dimg.h90 – NEMO

Ignore:
Timestamp:
2007-10-10T17:36:03+02:00 (17 years ago)
Author:
smasson
Message:

tau(x|y) -> (u|v)tau, see ticket:7

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/DIA/diawri_dimg.h90

    r699 r710  
    4141    !!    To be tested with a lot of procs !!!! 
    4242    !! 
    43     !!  level 1:  taux(:,:) * umask(:,:,1) zonal stress in N.m-2 
    44     !!  level 2:  tauy(:,:) * vmask(:,:,1) meridional stress in N. m-2 
     43    !!  level 1:  utau(:,:) * umask(:,:,1) zonal stress in N.m-2 
     44    !!  level 2:  vtau(:,:) * vmask(:,:,1) meridional stress in N. m-2 
    4545    !!  level 3:   qt  (:,:)               total heat flux (W/m2) 
    4646    !!  level 4:   emp (:,:)               E-P flux (mm/day) 
     
    167167       sm(:,:,:)=sm(:,:,:) + sn (:,:,:) 
    168168       ! 
    169        fsel(:,:,1 ) = fsel(:,:,1 ) + taux(:,:) * umask(:,:,1) 
    170        fsel(:,:,2 ) = fsel(:,:,2 ) + tauy(:,:) * vmask(:,:,1) 
     169       fsel(:,:,1 ) = fsel(:,:,1 ) + utau(:,:) * umask(:,:,1) 
     170       fsel(:,:,2 ) = fsel(:,:,2 ) + vtau(:,:) * vmask(:,:,1) 
    171171       fsel(:,:,3 ) = fsel(:,:,3 ) + qt  (:,:)  
    172172       fsel(:,:,4 ) = fsel(:,:,4 ) + emp (:,:) 
     
    231231          fsel(:,:,20)= spgv(:,:) 
    232232#endif 
    233           ! mask mean field with tmask except taux tauy (1,2) 
     233          ! mask mean field with tmask except utau vtau (1,2) 
    234234          DO jk=3,inbsel 
    235235            fsel(:,:,jk)=fsel(:,:,jk)*tmask(:,:,1) 
     
    256256          fsel(:,:,:) = 0._wp 
    257257          ! 
    258           fsel(:,:,1 ) = taux(:,:) * umask(:,:,1) 
    259           fsel(:,:,2 ) = tauy(:,:) * vmask(:,:,1) 
     258          fsel(:,:,1 ) = utau(:,:) * umask(:,:,1) 
     259          fsel(:,:,2 ) = vtau(:,:) * vmask(:,:,1) 
    260260          fsel(:,:,3 ) = qt  (:,:) * tmask(:,:,1) 
    261261          fsel(:,:,4 ) = emp (:,:) * tmask(:,:,1) 
Note: See TracChangeset for help on using the changeset viewer.