Changeset 3489


Ignore:
Timestamp:
2012-10-05T17:10:13+02:00 (8 years ago)
Author:
acc
Message:

Branch: dev_r3385_NOCS04_HAMF; #665. Stage 4 of 2012 development: Tidying of diagnostic output and bug correction in limsbc_2.F90

Location:
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef.xml

    r3294 r3489  
    124124 
    125125   <field id="empmr"        description="Net Upward Water Flux"                                        unit="kg/m2/s"  /> 
    126    <field id="empsmr"       description="concentration/dilution water flux"                            unit="kg/m2/s" /> 
     126   <field id="saltflx"      description="Downward Salt Flux"                                           unit="PSU/m2/s" /> 
    127127   <field id="snowpre"      description="Snow precipitation"                                           unit="kg/m2/s"  /> 
    128128   <field id="runoffs"      description="River Runoffs"                                                unit="Kg/m2/s"  /> 
     
    145145   <field id="qsb_oce"      description="Sensible Downward Heat Flux over open ocean"                  unit="W/m2"     /> 
    146146   <field id="qla_oce"      description="Latent Downward Heat Flux over open ocean"                    unit="W/m2"     /> 
     147   <field id="qhc_oce"      description="Downward Heat Content of E-P over open ocean"                 unit="W/m2"     /> 
    147148   <field id="taum_oce"     description="wind stress module over open ocean"                           unit="N/m2"     /> 
    148149 
     
    173174   <field id="v_imasstr"    description="Sea-ice mass transport along j-axis"                          unit="kg/s"     /> 
    174175 
     176   <!-- available if not defined key_vvl --> 
     177   <field id="emp_x_sst"      description="Concentration/Dilution term on SST"                         unit="kgC/m2/s" /> 
     178   <field id="emp_x_sss"      description="Concentration/Dilution term on SSS"                       unit="kgPSU/m2/s" /> 
    175179   <!-- available key_coupled --> 
    176180   <field id="snow_ao_cea"  description="Snow over ice-free ocean (cell average)"                      unit="kg/m2/s"  /> 
     
    10161020     <field ref="empmr"        name="sowaflup"  /> 
    10171021     <field ref="qsr"          name="soshfldo"  /> 
    1018      <field ref="empsmr"       name="sowaflcd"  /> 
     1022     <field ref="saltflx"      name="sosfldow"  /> 
    10191023     <field ref="qt"           name="sohefldo"  /> 
    10201024     <field ref="mldr10_1"     name="somxl010"  /> 
  • branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90

    r3488 r3489  
    192192 
    193193            fsbbq(ji,jj) = ( 1.0 - ( ifvt + iflt ) ) * fscmbq(ji,jj)     ! store residual heat flux (to put into the ocean at the next time-step) 
    194             zqhc = ( rdq_snw(ji,jj) + rdq_ice(ji,jj) ) * r1_rdtice       ! heat flux due to snow & ice heat content exchange 
    195             ! 
     194            zqhc = ( rdq_snw(ji,jj)                                     & 
     195                 & + rdq_ice(ji,jj) * ( 1.- zswitch) ) * r1_rdtice       ! heat flux due to snow ( & ice heat content,  
     196            !                                                            !           if ice/ocean mass exchange active)  
    196197            qsr  (ji,jj) = zqsr                                          ! solar heat flux  
    197             qns  (ji,jj) = zqns - fdtcn(ji,jj) + zqhc                    ! non solar heat flux 
     198            qns  (ji,jj) = zqns - fdtcn(ji,jj) + zqhc                    ! non solar heat flux  
     199            ! 
    198200            !                          !------------------------------------------! 
    199201            !                          !      mass flux at the ocean surface      ! 
     
    223225            ! salt and mass fluxes : non concentration dilution term in salt flux  and add F-M term in volume flux 
    224226            sfx (ji,jj) = zfsalt +                  zswitch  * zcd   ! salt flux (+ C/D if no ice/ocean mass exchange) 
    225             emp (ji,jj) = zemp   + zemp_snw + ( 1.- zswitch) * zfmm  ! mass flux (- F/M mass flux if no ice/ocean mass exchange) 
     227            emp (ji,jj) = zemp   + zemp_snw + ( 1.- zswitch) * zfmm  ! mass flux (+ F/M mass flux if ice/ocean mass exchange) 
    226228            ! 
    227229         END DO 
  • branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r3488 r3489  
    372372         CALL histdef( nid_T, "sossheig", "Sea Surface Height"                 , "m"      ,   &  ! ssh 
    373373            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    374 >>>>>>>>>>> 
    375       CALL histwrite( nid_T, "sowaflup", it, ( emp-rnf )   , ndim_hT, ndex_hT )   ! upward water flux 
    376       CALL histwrite( nid_T, "sosfldow", it, sfx           , ndim_hT, ndex_hT )   ! downward salt flux  
    377                                                                                   ! (includes virtual salt flux beneath ice  
    378                                                                                   ! in linear free surface case) 
     374         CALL histdef( nid_T, "sowaflup", "Net Upward Water Flux"              , "Kg/m2/s",   &  ! (emp-rnf) 
     375            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
     376         CALL histdef( nid_T, "sosfldow", "downward salt flux"                 , "PSU/m2/s",  &  ! sfx 
     377            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    379378#if ! defined key_vvl 
    380       zw2d(:,:) = emp (:,:) * tsn(:,:,1,jp_tem) * tmask(:,:,1) 
    381       CALL histwrite( nid_T, "sosst_cd", it, zw2d, ndim_hT, ndex_hT )             ! c/d term on sst 
    382       zw2d(:,:) = emp (:,:) * tsn(:,:,1,jp_sal) * tmask(:,:,1) 
    383       CALL histwrite( nid_T, "sosss_cd", it, zw2d, ndim_hT, ndex_hT )             ! c/d term on sss 
    384 #endif 
    385 <<<<<<<<<< 
    386          CALL histdef( nid_T, "sowaflup", "Net Upward Water Flux"              , "Kg/m2/s",   &  ! (emp-rnf) 
    387             &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    388          CALL histdef( nid_T, "sofldow", "downward salt flux"                  , "PSU/m2/s",  &  ! sfx 
    389             &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    390          CALL histdef( nid_T, "sosalflx", "Surface Salt Flux"                  , "Kg/m2/s",   &  ! (emps-rnf) * sn 
    391             &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
     379         CALL histdef( nid_T, "sosst_cd", "Concentration/Dilution term on temperature"        &  ! emp * tsn(:,:,1,jp_tem) 
     380            &                                                                  , "KgC/m2/s",  &  ! sosst_cd 
     381            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
     382         CALL histdef( nid_T, "sosss_cd", "Concentration/Dilution term on salinity"           &  ! emp * tsn(:,:,1,jp_sal) 
     383            &                                                                  , "KgPSU/m2/s",&  ! sosss_cd 
     384            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
     385#endif 
    392386         CALL histdef( nid_T, "sohefldo", "Net Downward Heat Flux"             , "W/m2"   ,   &  ! qns + qsr 
    393387            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
  • branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/DIA/diawri_dimg.h90

    r3402 r3489  
    5454    !!  level 14: qct(:,:)                 equivalent flux due to treshold SST 
    5555    !!  level 15: fbt(:,:)                 feedback term . 
    56     !!  level 16: ( sfx (:,:) - rnf(:,:) ) concentration/dilution water flux 
     56    !!  level 16: ( emp * sss )            concentration/dilution term on salinity 
     57    !!  level 17: ( emp * sst )            concentration/dilution term on temperature 
    5758    !!  level 17: fsalt(:,:)               Ice=>ocean net freshwater 
    5859    !!  level 18: gps(:,:)                 the surface pressure (m). 
     
    6263    !! History:  OPA  ! 1997-02 ( Clipper Group ) dimg files 
    6364    !!            -   ! 2003-12 ( J.M. Molines) f90, mpp output for OPA9.0 
    64     !!   NEMO    1.0  ! 2005-05  (S. Theetten) add sfx fsalt move gps spgu spgv 2 lines below 
     65    !!   NEMO    1.0  ! 2005-05  (S. Theetten) add emps fsalt move gps spgu spgv 2 lines below 
    6566    !!            -   ! 2005-11  (V. Garnier) Surface pressure gradient organization 
    6667    !!---------------------------------------------------------------------- 
     
    107108 
    108109 
    109     inbsel = 17 
     110    inbsel = 18 
    110111 
    111112    IF( inbsel >  jpk ) THEN 
     
    172173       !        fsel(:,:,14) = fsel(:,:,14) + qct(:,:) 
    173174       !        fsel(:,:,15) = fsel(:,:,15) + fbt(:,:) 
    174        fsel(:,:,16) = fsel(:,:,16) + ( sfx (:,:)-rnf(:,:) )  
     175       fsel(:,:,16) = fsel(:,:,16) + ( emp(:,:)*tsn(:,:,1,jp_sal) )  
     176       fsel(:,:,17) = fsel(:,:,17) + ( emp(:,:)*tsn(:,:,1,jp_tem) )  
    175177       ! 
    176178       ! Output of dynamics and tracer fields and selected fields 
     
    240242          !         fsel(:,:,14) =  qct(:,:) 
    241243          !         fsel(:,:,15) =  fbt(:,:) 
    242           fsel(:,:,16) = ( sfx (:,:)-rnf(:,:) ) * tmask(:,:,1)  
     244          fsel(:,:,16) = ( emp(:,:)-tsn(:,:,1,jp_sal) ) * tmask(:,:,1)  
     245          fsel(:,:,17) = ( emp(:,:)-tsn(:,:,1,jp_tem) ) * tmask(:,:,1)  
    243246          ! 
    244247          !         qct(:,:) = 0._wp 
  • branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r3402 r3489  
    392392      CALL iom_put( "qsb_oce", - zqsb )                 ! output downward sensible heat over the ocean 
    393393      CALL iom_put( "qla_oce", - zqla )                 ! output downward latent   heat over the ocean 
     394      CALL iom_put( "qhc_oce",   qns-zqlw+zqsb+zqla )   ! output downward heat content of E-P over the ocean 
    394395      CALL iom_put( "qns_oce",   qns  )                 ! output downward non solar heat over the ocean 
    395396      ! 
Note: See TracChangeset for help on using the changeset viewer.