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 15434 for NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/DIA/diaregmean.F90 – NEMO

Ignore:
Timestamp:
2021-10-22T12:32:56+02:00 (3 years ago)
Author:
hadjt
Message:

Region mean and NOOS cross-sections

Use rcp (Specific heat capacity) from nemo physical constant

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/DIA/diaregmean.F90

    r15411 r15434  
    1010   USE in_out_manager  ! I/O units 
    1111   USE iom             ! I/0 library 
     12   USE phycst          ! physical constants 
    1213   USE diapea          ! PEA 
    1314   USE zdfmxl          ! MLD 
     
    589590        DEALLOCATE(  tmp1mat ) 
    590591 
    591         tmp_field_HSVM_mat(:,:,1) = (zwtmbT(:,:,6)*tmask(:,:,1)*3850.) !heat 4200 is value for FW, 3850 is the value for sea water.  
     592        !tmp_field_HSVM_mat(:,:,1) = (zwtmbT(:,:,6)*tmask(:,:,1)*3850.) !heat 3850/4200 is value for FW, 3850 is the value for sea 
     593        tmp_field_HSVM_mat(:,:,1) = (zwtmbT(:,:,6)*tmask(:,:,1)*rcp) !heat 4200 is value for FW, 3850 is the value for sea water.  
    592594        tmp_field_HSVM_mat(:,:,2) = (zwtmbS(:,:,6)*tmask(:,:,1))       !salt 
    593595        tmp_field_HSVM_mat(:,:,3) = (zwtmb1(:,:,5)*tmask(:,:,1))       !vol 
     
    692694            name_SBC_mat(7) = 'rnf' 
    693695 
    694             tmp_field_SBC_mat(:,:,8) = tmp_field_SBC_mat(:,:,8) + (emp*tmask(:,:,1)*tsn(:,:,1,jp_tem)*3850.) 
     696            !tmp_field_SBC_mat(:,:,8) = tmp_field_SBC_mat(:,:,8) + (emp*tmask(:,:,1)*tsn(:,:,1,jp_tem)*3850.) 
     697            tmp_field_SBC_mat(:,:,8) = tmp_field_SBC_mat(:,:,8) + (emp*tmask(:,:,1)*tsn(:,:,1,jp_tem)*rcp) 
    695698            name_SBC_mat(8) = 'empheat' 
    696             tmp_field_SBC_mat(:,:,9) = tmp_field_SBC_mat(:,:,9) + (rnf*tmask(:,:,1)*tsn(:,:,1,jp_tem)*3850.) 
     699            !tmp_field_SBC_mat(:,:,9) = tmp_field_SBC_mat(:,:,9) + (rnf*tmask(:,:,1)*tsn(:,:,1,jp_tem)*3850.) 
     700            tmp_field_SBC_mat(:,:,9) = tmp_field_SBC_mat(:,:,9) + (rnf*tmask(:,:,1)*tsn(:,:,1,jp_tem)*rcp) 
    697701            name_SBC_mat(9) = 'rnfheat' 
    698702 
Note: See TracChangeset for help on using the changeset viewer.