Changeset 39 for branches/2017


Ignore:
Timestamp:
01/18/17 11:05:37 (7 years ago)
Author:
vancop
Message:

add tank mass and salt balance

Location:
branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20
Files:
2 added
1 deleted
11 edited

Legend:

Unmodified
Added
Removed

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice.com

    r30 r39  
    182182 
    183183      COMMON/barrowconf/sal_read(11), hi_read(11), hgins, hnins, 
    184      &              tsuins, s_w, oce_flx, num_sal, nday1, 
     184     &              tsuins, oce_flx, num_sal, nday1, 
    185185     &              i_sal 
    186186      

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_brine.f

    r27 r39  
    2020      INCLUDE 'thermo.com' 
    2121      INCLUDE 'bio.com' 
     22      INCLUDE 'tank.com' 
    2223 
    2324!==============================================================================! 

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_output.f

    r29 r39  
    1313      INCLUDE 'bio.com' 
    1414      INCLUDE 'forcing.com' 
     15      INCLUDE 'tank.com' 
    1516 
    1617      CHARACTER(len=10) ::  
    … …  
    237238         WRITE(numout,*) ' fsbp dimension created ' 
    238239 
    239          CALL CF_CREATE_VAR("s_w","Ocean Salinity", 
    240      .                   "kgNaCl/dm3","time","-","-","-" ) 
     240         CALL CF_CREATE_VAR("s_w","Water Salinity", 
     241     .                   "g/kg","time","-","-","-" ) 
    241242         WRITE(numout,*) ' s_w dimension created ' 
     243 
     244         CALL CF_CREATE_VAR("v_w","Water Volume", 
     245     .                   "m3","time","-","-","-" ) 
     246         WRITE(numout,*) ' v_w dimension created ' 
    242247 
    243248         ! Forcing 
    … …  
    621626      CALL CF_WRITE (filenc, 's_w', numit-nstart+1,  
    622627     &               1, 1, 1, REAL(s_w) ) 
     628      CALL CF_WRITE (filenc, 'v_w', numit-nstart+1,  
     629     &               1, 1, 1, REAL(v_w) ) 
    623630 
    624631      ! Forcing 

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_phy_ini.f

    r27 r39  
    1414      INCLUDE 'thermo.com' 
    1515      INCLUDE 'bio.com' 
     16      INCLUDE 'tank.com' 
    1617 
    1718      CHARACTER(len=10) ::  
    … …  
    4041      zeps = 1.0e-10 
    4142      ji = 1 
     43! 
     44!------------------------------------------------------------------------------| 
     45!  1) Initialize tank       
     46!------------------------------------------------------------------------------| 
     47! 
     48 
     49      ! ASIBIA tank parameters 
     50      a_tank = 2.35 * 1.35               ! area of the surface of the tank in m2 
     51 
     52      v_w    = a_tank * 1.               ! initial volume of water in the tank in m3 (without side tank - needs adjustment) 
     53      s_w    = 28.33866                  ! initial water salinity (g/kg) 
     54      ms_w   = rho0 * v_w * s_w / 1000.  ! initial mass of salt in water (kg NaCl) 
    4255 
    4356! 
    4457!------------------------------------------------------------------------------| 
    45 !  1) Case of a NETCDF file 
     58!  2) Initialize sea ice from the netcdf file 
    4659!------------------------------------------------------------------------------| 
    4760! 
    … …  
    148161              t_i_b(ji,layer) = tpw + (- bbb - discrim) / ( 2.0*aaa ) 
    149162           END DO 
    150  
    151163! 
    152164!------------------------------------------------------------------------------| 
    153165! x) Write in output file 
    154166!------------------------------------------------------------------------------| 
     167! 
    155168 
    156169            WRITE(numout,*) ' --- Initial values --- ' 

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_phy_param.f

    r30 r39  
    2929      INCLUDE 'para.com' 
    3030      INCLUDE 'ice.com' 
     31      INCLUDE 'tank.com' 
    3132 
    3233      ! name of the experiment 
    … …  
    108109      rho0    = 1025.0            ! ocean mean density 
    109110      cpw     = 3.99d+03          ! seawater specific heat 
    110       s_w     = 34.0              ! 34. is the control 
    111111 
    112112      visc_br = 1.9e-3            ! dynamic viscosity of water at 0C (kg/(m.s)) - Sharqawy et al, 2010, eq. 22 

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_phy_remap.f

    r27 r39  
    3232      INCLUDE 'thermo.com' 
    3333      INCLUDE 'bio.com' 
     34      INCLUDE 'tank.com' 
    3435 
    3536      ! Local Variables 
    … …  
    578579     &                z_ms_i_fin 
    579580 
     581      dms_i = z_ms_i_fin / 1000. - dms_i         ! 2nd step of computation of the change in mass of salt in sea ice 
     582 
    580583      ! Flux due to bottom and snow ice formation 
    581584      zfgrml = zsh_i_new / ddtb 

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_sal_adv.f

    r30 r39  
    3838      INCLUDE 'ice.com' 
    3939      INCLUDE 'thermo.com' 
     40      INCLUDE 'tank.com' 
    4041 
    4142      REAL(8), DIMENSION(nlay_i) ::   

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_sal_diff.f

    r27 r39  
    3434      INCLUDE 'ice.com' 
    3535      INCLUDE 'thermo.com' 
     36      INCLUDE 'tank.com' 
    3637 
    3738      REAL(8), DIMENSION(nlay_i) ::   
    … …  
    139140     &                    deltaz_i_phy, nlay_i, .FALSE. ) 
    140141      ENDIF ! ln_con 
     142 
     143      dms_i = z_ms_i_ini / 1000.         ! 1st step of computation of the change in mass of salt in sea ice 
    141144 
    142145      IF ( ln_write ) THEN 

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_sal_diff_CW.f

    r27 r39  
    3434      INCLUDE 'ice.com' 
    3535      INCLUDE 'thermo.com' 
     36      INCLUDE 'tank.com' 
    3637 
    3738      REAL(8), DIMENSION(nlay_i) ::   

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_th.f

    r36 r39  
    3131      INCLUDE 'thermo.com' 
    3232      INCLUDE 'bio.com' 
     33      INCLUDE 'tank.com' 
    3334 
    3435      ! Energy conservation 
    … …  
    8384     &                                 q_i_layer_in , 1 , 1 , jl,  
    8485     &                                 nlay_s , nlay_i ) 
     86 
    8587! 
    8688!------------------------------------------------------------------------------- 
    … …  
    150152     &                   layer = 1, nlay_bio ) 
    151153 
     154      !----------------- 
     155      ! CALL water tank 
     156      !----------------- 
     157      CALL wat_tank 
     158      WRITE(numout,*) ' --- ice_th --- ' 
     159      WRITE(numout,*) ' s_w : ', s_w 
     160      WRITE(numout,*) ' -------------- ' 
     161 
    152162!------------------------------------------------------------------------------- 
    153163!  4) Outputs  

  • branches/2017/dev_v3.20_2017_transport_max/SOURCES/source_3.20/ice_th_dh.f

    r34 r39  
    4646      INCLUDE 'ice.com' 
    4747      INCLUDE 'thermo.com' 
     48      INCLUDE 'tank.com' 
    4849 
    4950      ! Local Variables 
Note: See TracChangeset for help on using the changeset viewer.