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 714 for trunk/NEMO – NEMO

Changeset 714 for trunk/NEMO


Ignore:
Timestamp:
2007-10-11T11:58:02+02:00 (16 years ago)
Author:
smasson
Message:

replace qt by qns, see ticket:11

Location:
trunk/NEMO
Files:
7 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/C1D_SRC/diawri1d.F90

    r712 r714  
    210210         CALL histdef( nid_T, "sosalflx", "Surface Salt Flux"                  , "Kg/m2/s",   &  ! emps * sn 
    211211            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    212          CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux"             , "W/m2"   ,   &  ! qt 
     212         CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux"             , "W/m2"   ,   &  ! qsr + qns 
    213213            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    214214         CALL histdef( nid_T, "soshfldo", "Shortwave Radiation"                , "W/m2"   ,   &  ! qsr 
     
    373373      zw2d(:,:) = emps(:,:) * sn(:,:,1) * tmask(:,:,1) 
    374374      CALL histwrite( nid_T, "sosalflx", it, zw2d          , ndim_hT, ndex_hT )   ! c/d salt flux 
    375       CALL histwrite( nid_T, "sohefldo", it, qt            , ndim_hT, ndex_hT )   ! total heat flux 
     375      CALL histwrite( nid_T, "sohefldo", it, qsr + qns     , ndim_hT, ndex_hT )   ! total heat flux 
    376376      CALL histwrite( nid_T, "soshfldo", it, qsr           , ndim_hT, ndex_hT )   ! solar heat flux 
    377377      CALL histwrite( nid_T, "somxl010", it, hmlp          , ndim_hT, ndex_hT )   ! mixed layer depth 
  • trunk/NEMO/C1D_SRC/icestp1d.F90

    r710 r714  
    214214      IF( kt == nit000 ) THEN      
    215215         qsr    (:,:) = 0.e0 
    216          qt     (:,:) = 0.e0 
     216         qns    (:,:) = 0.e0 
    217217         qrp    (:,:) = 0.e0 
    218218         emp    (:,:) = 0.e0 
     
    232232      ! ----------------- 
    233233       
    234       qt  (:,:) = fnsolar(:,:) + fsolar(:,:)     ! non solar heat flux + solar flux 
     234      qns (:,:) = fnsolar(:,:)                    ! non solar heat flux 
    235235      qsr (:,:) = fsolar(:,:)                     ! solar flux 
    236236       
  • trunk/NEMO/C1D_SRC/step1d.F90

    r710 r714  
    152152         CALL prt_ctl(tab2d_1=emp    , clinfo1=' emp  -   : ', mask1=tmask, ovlap=1) 
    153153         CALL prt_ctl(tab2d_1=emps   , clinfo1=' emps -   : ', mask1=tmask, ovlap=1) 
    154          CALL prt_ctl(tab2d_1=qt     , clinfo1=' qt   -   : ', mask1=tmask, ovlap=1) 
     154         CALL prt_ctl(tab2d_1=qns    , clinfo1=' qns  -   : ', mask1=tmask, ovlap=1) 
    155155         CALL prt_ctl(tab2d_1=qsr    , clinfo1=' qsr  -   : ', mask1=tmask, ovlap=1) 
    156156         CALL prt_ctl(tab2d_1=runoff , clinfo1=' runoff   : ', mask1=tmask, ovlap=1) 
  • trunk/NEMO/OPA_SRC/DIA/diawri.F90

    r712 r714  
    260260         CALL histdef( nid_T, "sosalflx", "Surface Salt Flux"                  , "Kg/m2/s",   &  ! emps * sn 
    261261            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    262          CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux"             , "W/m2"   ,   &  ! qt 
     262         CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux"             , "W/m2"   ,   &  ! qns + qsr 
    263263            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    264264         CALL histdef( nid_T, "soshfldo", "Shortwave Radiation"                , "W/m2"   ,   &  ! qsr 
     
    430430      zw2d(:,:) = emps(:,:) * sn(:,:,1) * tmask(:,:,1) 
    431431      CALL histwrite( nid_T, "sosalflx", it, zw2d          , ndim_hT, ndex_hT )   ! c/d salt flux 
    432       CALL histwrite( nid_T, "sohefldo", it, qt            , ndim_hT, ndex_hT )   ! total heat flux 
     432      CALL histwrite( nid_T, "sohefldo", it, qns + qsr     , ndim_hT, ndex_hT )   ! total heat flux 
    433433      CALL histwrite( nid_T, "soshfldo", it, qsr           , ndim_hT, ndex_hT )   ! solar heat flux 
    434434      CALL histwrite( nid_T, "somxl010", it, hmlp          , ndim_hT, ndex_hT )   ! mixed layer depth 
     
    638638 
    639639      ! Write all fields on T grid 
    640       CALL histwrite( id_i, "votemper", 1, tn    , jpi*jpj*jpk, idex )    ! now temperature 
    641       CALL histwrite( id_i, "vosaline", 1, sn    , jpi*jpj*jpk, idex )    ! now salinity 
    642 #if defined key_dynspg_rl 
    643       CALL histwrite( id_i, "sobarstf", 1, bsfn  , jpi*jpj    , idex )    ! barotropic streamfunction 
     640      CALL histwrite( id_i, "votemper", 1, tn      , jpi*jpj*jpk, idex )    ! now temperature 
     641      CALL histwrite( id_i, "vosaline", 1, sn        , jpi*jpj*jpk, idex )    ! now salinity 
     642#if defined key_dynspg_rl 
     643      CALL histwrite( id_i, "sobarstf", 1, bsfn     , jpi*jpj    , idex )    ! barotropic streamfunction 
    644644#else 
    645       CALL histwrite( id_i, "sossheig", 1, sshn  , jpi*jpj    , idex )    ! sea surface height 
    646 #endif 
    647       CALL histwrite( id_i, "vozocrtx", 1, un    , jpi*jpj*jpk, idex )    ! now i-velocity 
    648       CALL histwrite( id_i, "vomecrty", 1, vn    , jpi*jpj*jpk, idex )    ! now j-velocity 
    649       CALL histwrite( id_i, "vovecrtz", 1, wn    , jpi*jpj*jpk, idex )    ! now k-velocity 
    650       CALL histwrite( id_i, "sowaflup", 1, emp   , jpi*jpj    , idex )    ! freshwater budget 
    651       CALL histwrite( id_i, "sohefldo", 1, qt    , jpi*jpj    , idex )    ! total heat flux 
    652       CALL histwrite( id_i, "soshfldo", 1, qsr   , jpi*jpj    , idex )    ! total heat flux 
    653       CALL histwrite( id_i, "soicecov", 1, freeze, jpi*jpj    , idex )    ! ice cover 
    654       CALL histwrite( id_i, "sozotaux", 1, utau  , jpi*jpj    , idex )    ! i-wind stress 
    655       CALL histwrite( id_i, "sometauy", 1, vtau  , jpi*jpj    , idex )    ! j-wind stress 
     645      CALL histwrite( id_i, "sossheig", 1, sshn     , jpi*jpj    , idex )    ! sea surface height 
     646#endif 
     647      CALL histwrite( id_i, "vozocrtx", 1, un       , jpi*jpj*jpk, idex )    ! now i-velocity 
     648      CALL histwrite( id_i, "vomecrty", 1, vn       , jpi*jpj*jpk, idex )    ! now j-velocity 
     649      CALL histwrite( id_i, "vovecrtz", 1, wn       , jpi*jpj*jpk, idex )    ! now k-velocity 
     650      CALL histwrite( id_i, "sowaflup", 1, emp      , jpi*jpj    , idex )    ! freshwater budget 
     651      CALL histwrite( id_i, "sohefldo", 1, qsr + qns, jpi*jpj    , idex )    ! total heat flux 
     652      CALL histwrite( id_i, "soshfldo", 1, qsr      , jpi*jpj    , idex )    ! solar heat flux 
     653      CALL histwrite( id_i, "soicecov", 1, freeze   , jpi*jpj    , idex )    ! ice cover 
     654      CALL histwrite( id_i, "sozotaux", 1, utau     , jpi*jpj    , idex )    ! i-wind stress 
     655      CALL histwrite( id_i, "sometauy", 1, vtau     , jpi*jpj    , idex )    ! j-wind stress 
    656656 
    657657      ! 3. Close the file 
  • trunk/NEMO/OPA_SRC/DIA/diawri_dimg.h90

    r710 r714  
    4343    !!  level 1:  utau(:,:) * umask(:,:,1) zonal stress in N.m-2 
    4444    !!  level 2:  vtau(:,:) * vmask(:,:,1) meridional stress in N. m-2 
    45     !!  level 3:   qt  (:,:)               total heat flux (W/m2) 
    46     !!  level 4:   emp (:,:)               E-P flux (mm/day) 
     45    !!  level 3:  qsr + qns                total heat flux (W/m2) 
     46    !!  level 4:  emp (:,:)               E-P flux (mm/day) 
    4747    !!  level 5:  tb  (:,:,1)-sst          model SST -forcing sst (degree C) 
    4848    !!  level 6:  bsfb(:,:)         streamfunction (m**3/s) 
     
    169169       fsel(:,:,1 ) = fsel(:,:,1 ) + utau(:,:) * umask(:,:,1) 
    170170       fsel(:,:,2 ) = fsel(:,:,2 ) + vtau(:,:) * vmask(:,:,1) 
    171        fsel(:,:,3 ) = fsel(:,:,3 ) + qt  (:,:)  
     171       fsel(:,:,3 ) = fsel(:,:,3 ) + qsr (:,:) + qns  (:,:)  
    172172       fsel(:,:,4 ) = fsel(:,:,4 ) + emp (:,:) 
    173173       fsel(:,:,5 ) = fsel(:,:,5 ) + tb  (:,:,1) - sst(:,:) 
     
    258258          fsel(:,:,1 ) = utau(:,:) * umask(:,:,1) 
    259259          fsel(:,:,2 ) = vtau(:,:) * vmask(:,:,1) 
    260           fsel(:,:,3 ) = qt  (:,:) * tmask(:,:,1) 
     260          fsel(:,:,3 ) = (qsr (:,:) + qnr (:,:)) * tmask(:,:,1) 
    261261          fsel(:,:,4 ) = emp (:,:) * tmask(:,:,1) 
    262262          fsel(:,:,5 ) = (tb  (:,:,1) -sst(:,:)) *tmask(:,:,1) 
  • trunk/NEMO/OPA_SRC/TRA/trasbc.F90

    r708 r714  
    132132#endif 
    133133            IF( lk_vvl) THEN 
    134                zta = ro0cpr * ( qt(ji,jj) - qsr(ji,jj) ) * zse3t &   ! temperature : heat flux 
     134               zta = ro0cpr * qns(ji,jj) * zse3t &                   ! temperature : heat flux 
    135135                &    - emp(ji,jj) * zsrau * tn(ji,jj,1)  * zse3t     ! & cooling/heating effet of EMP flux 
    136136               zsa = 0.e0                                            ! No salinity concent./dilut. effect 
    137137            ELSE 
    138                zta = ro0cpr * ( qt(ji,jj) - qsr(ji,jj) ) * zse3t     ! temperature : heat flux 
     138               zta = ro0cpr * qns(ji,jj) * zse3t     ! temperature : heat flux 
    139139               zsa = emps(ji,jj) * zsrau * sn(ji,jj,1)   * zse3t     ! salinity :  concent./dilut. effect 
    140140            ENDIF 
  • trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90

    r710 r714  
    459459            zBosol(ji,jj) = grav * zthermal * qsr(ji,jj) 
    460460            ! Non radiative surface buoyancy force 
    461             zBo   (ji,jj) = grav * zthermal * ( qt(ji,jj) - qsr(ji,jj) ) -  grav * zhalin * emp(ji,jj) 
     461            zBo   (ji,jj) = grav * zthermal * qns(ji,jj) -  grav * zhalin * emp(ji,jj) 
    462462            ! Surface Temperature flux for non-local term 
    463             wt0(ji,jj) = - qt(ji,jj) * ro0cpr * tmask(ji,jj,1) 
     463            wt0(ji,jj) = - ( qsr(ji,jj) + qns(ji,jj) )* ro0cpr * tmask(ji,jj,1) 
    464464            ! Surface salinity flux for non-local term 
    465465            ws0(ji,jj) = - ( emp(ji,jj) * sn(ji,jj,1) * rcs ) * tmask(ji,jj,1) 
Note: See TracChangeset for help on using the changeset viewer.