Changeset 731


Ignore:
Timestamp:
2007-11-19T16:10:23+01:00 (14 years ago)
Author:
cetlod
Message:

computation of horizontal divergence needed for vertical active trends, see ticket:12

File:
1 edited

Legend:

Unmodified
Added
Removed
  • tags/nemo_v2_3/NEMO/OFF_SRC/dtadyn.F90

    r612 r731  
    6767      vdta   ,   & ! meridional velocity at two consecutive times 
    6868      wdta   ,   & ! vertical velocity at two consecutive times 
     69#if defined key_trc_diatrd 
     70      hdivdta,   & ! horizontal divergence 
     71#endif 
    6972      avtdta       ! vertical diffusivity coefficient 
     73 
    7074 
    7175#if defined key_ldfslp 
     
    120124   !!---------------------------------------------------------------------- 
    121125   !!   OPA 9.0 , LOCEAN-IPSL  (2005) 
    122    !!   $Header$ 
     126   !!   $Header: /home/opalod/NEMOCVSROOT/NEMO/OFF_SRC/dtadyn.F90,v 1.8 2007/10/12 09:05:29 opalod Exp $ 
    123127   !!   This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
    124128   !!---------------------------------------------------------------------- 
     
    319323                vdta(:,:,:,1)=vdta(:,:,:,2) 
    320324                wdta(:,:,:,1)=wdta(:,:,:,2) 
     325#if defined key_trc_diatrd 
     326                hdivdta(:,:,:,1)=hdivdta(:,:,:,2) 
     327#endif 
    321328                avtdta(:,:,:,1)=avtdta(:,:,:,2) 
    322329                tdta(:,:,:,1)=tdta(:,:,:,2) 
     
    412419                udta(:,:,:,1) = udta(:,:,:,2) 
    413420                vdta(:,:,:,1) = vdta(:,:,:,2) 
    414                 wdta(:,:,:,1)= wdta(:,:,:,2) 
     421                wdta(:,:,:,1) = wdta(:,:,:,2) 
     422#if defined key_trc_diatrd 
     423                hdivdta(:,:,:,1)=hdivdta(:,:,:,2) 
     424#endif 
    415425                avtdta(:,:,:,1) = avtdta(:,:,:,2) 
    416426                tdta(:,:,:,1) = tdta(:,:,:,2) 
     
    508518                    vn(:,:,:)=vdta(:,:,:,2) 
    509519                    wn(:,:,:)=wdta(:,:,:,2) 
     520#if defined key_trc_diatrd 
     521                    hdivn(:,:,:)=hdivdta(:,:,:,2) 
     522#endif 
    510523#if defined key_trc_zdfddm 
    511524                    avs(:,:,:)=avtdta(:,:,:,2) 
     
    570583                    vn(:,:,:) = zweighm1 * vdta(:,:,:,1) + zweigh * vdta(:,:,:,2) 
    571584                    wn(:,:,:) = zweighm1 * wdta(:,:,:,1) + zweigh * wdta(:,:,:,2) 
     585#if defined key_trc_diatrd 
     586                    hdivn(:,:,:) = zweighm1 * hdivdta(:,:,:,1) + zweigh * hdivdta(:,:,:,2) 
     587#endif 
    572588#if defined key_zdfddm 
    573589                    avs(:,:,:)= zweighm1 * avtdta(:,:,:,1) + zweigh * avtdta(:,:,:,2) 
     
    745761      CALL iom_get ( numfl_t, jpdom_data, 'vosaline', zs   (:,:,:), jkenr ) 
    746762      CALL iom_get ( numfl_t, jpdom_data, 'somixhgt', zmld (:,:  ), jkenr ) 
    747       CALL iom_get ( numfl_t, jpdom_data, 'sowaflup', zemp (:,:  ), jkenr ) 
     763      CALL iom_get ( numfl_t, jpdom_data, 'sowaflcd', zemp (:,:  ), jkenr ) 
    748764      CALL iom_get ( numfl_t, jpdom_data, 'soshfldo', zqsr (:,:  ), jkenr ) 
    749765      CALL iom_get ( numfl_t, jpdom_data, 'soicecov', zice (:,:  ), jkenr ) 
     
    850866         zw(:,:,jk) = zw(:,:,jk+1) - fse3t(:,:,jk) * zhdiv(:,:,jk) 
    851867      END DO 
     868 
    852869      wdta(:,:,:,2) = zw(:,:,:) * tmask(:,:,:) 
    853  
     870#if defined key_trc_diatrd 
     871      hdivdta(:,:,:,2) = zhdiv(:,:,:) * tmask(:,:,:) 
     872#endif 
    854873 
    855874      tdta(:,:,:,2)   = zt(:,:,:)   * tmask(:,:,:) 
Note: See TracChangeset for help on using the changeset viewer.