- Timestamp:
- 2020-10-28T10:56:07+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/analyze_output.py
r13682 r13686 135 135 cv_rf_t = cv+'_tol' 136 136 cv_nemo = l_var_ot[jv] 137 print('\n\n ==== Checking flux '+cv_rf_m+' ('+cv_nemo+')!')137 print('\n\n ==== Checking variable '+cv_nemo+' against '+cv_rf_m+' !') 138 138 139 139 F_rf_m = nmp.zeros( Nt ) -
NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/sbcblk_sanity_check.sh
r13682 r13686 70 70 i_si3=0 71 71 FCPP="${SASF_WRK_DIR}/cpp_${CONFIG_BLD}.fcm" 72 if [ ! -f ${FCPP} ]; then echo " Mhhh, we did not find 'cpp_STATION_ASF.fcm' into `dirname ${FCPP}` !"; exit; fi72 if [ ! -f ${FCPP} ]; then echo " Mhhh, we could not find 'cpp_STATION_ASF.fcm' into `dirname ${FCPP}` !"; exit; fi 73 73 ca=`cat ${FCPP} | grep 'key_si3'` 74 74 … … 95 95 96 96 97 if [ ! -d ${DATA_IN_DIR} ]; then echo "PROBLEM!!! => did not find directory 'input_data' with input forcing..."; exit; fi 97 CFGS_SHARED="${NEMO_REF_DIR}/cfgs/SHARED" 98 if [ ! -d ${CFGS_SHARED} ]; then echo "PROBLEM!!! => could not find directory ${CFGS_SHARED} !"; exit; fi 99 100 if [ ! -d ${DATA_IN_DIR} ]; then echo "PROBLEM!!! => could not find directory 'input_data' with input forcing..."; exit; fi 98 101 99 102 cdt_cmpl="`\ls -l ${NEMO_EXE} | cut -d' ' -f 6,7,8`" … … 124 127 ################ 125 128 126 list_xml="iodef.xml file_def_nemo-oce.xml field_def_nemo-oce.xml" 129 list_xml_ref="field_def_nemo-oce.xml domain_def_nemo.xml grid_def_nemo.xml" 130 list_xml_cfg="iodef.xml file_def_nemo-oce.xml" 127 131 fcntxt="context_nemo_OCE.xml" 128 132 if [ ${i_sea_ice} -eq 1 ]; then 129 list_xml+=" file_def_nemo-ice.xml field_def_nemo-ice.xml" 133 list_xml_ref+=" field_def_nemo-ice.xml" 134 list_xml_cfg+=" file_def_nemo-ice.xml" 130 135 fcntxt="context_nemo_OCE+ICE.xml" 131 136 fi … … 136 141 137 142 # All remaining "*.xml" files: 138 for ff in ${list_xml } ; do143 for ff in ${list_xml_cfg} ; do 139 144 if [ ! -f ${SASF_EXPREF}/${ff} ]; then echo " Mhhh, ${ff} not found into ${SASF_EXPREF} !"; exit; fi 140 145 rsync -avPL ${SASF_EXPREF}/${ff} ${PROD_DIR}/ 141 146 done 142 143 144 # Getting reference/defaults files from reference NEMO distro: 145 rdir="${NEMO_REF_DIR}/cfgs/SHARED" 146 for ff in "domain_def_nemo.xml" "grid_def_nemo.xml"; do 147 if [ ! -f ${rdir}/${ff} ]; then echo " Mhhh, ${ff} not found into ${rdir} !"; exit; fi 148 ln -sf ${rdir}/${ff} ${PROD_DIR}/. 147 for ff in ${list_xml_ref} ; do 148 if [ ! -f ${CFGS_SHARED}/${ff} ]; then echo " Mhhh, ${ff} not found into ${CFGS_SHARED} !"; exit; fi 149 rsync -avPL ${CFGS_SHARED}/${ff} ${PROD_DIR}/ 149 150 done 150 151 … … 175 176 # The namelists: 176 177 rsync -avPL ${SASF_EXPREF}/${DIR_NL}/namelist_${scase}_cfg ${PROD_DIR}/namelist_cfg 177 rsync -avPL ${ SASF_EXPREF}/namelist_ref${PROD_DIR}/namelist_ref178 rsync -avPL ${CFGS_SHARED}/namelist_ref ${PROD_DIR}/namelist_ref 178 179 179 180 cd ${PROD_DIR}/ -
NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/MY_SRC/icestp.F90
r13655 r13686 74 74 !!---------------------------------------------------------------------- 75 75 !! NEMO/ICE 4.0 , NEMO Consortium (2018) 76 !! $Id: icestp.F90 136 41 2020-10-19 18:16:58Z clem$76 !! $Id: icestp.F90 13655 2020-10-21 14:15:13Z laurent $ 77 77 !! Software governed by the CeCILL license (see ./LICENSE) 78 78 !!---------------------------------------------------------------------- … … 127 127 ! utau_ice, vtau_ice = surface ice stress [N/m2] 128 128 !------------------------------------------------! 129 CALL ice_sbc_tau( kt, ksbc, utau_ice, vtau_ice ) 130 ! 129 131 !------------------------------------------------------! 130 132 ! --- Thermodynamical coupling with the atmosphere --- ! … … 150 152 !!============================================================= 151 153 152 !IF(lwp) WRITE(numout,*) ''153 !IF(lwp) WRITE(numout,*) 'LOLO:ice_stp()@icestp.F90 => calling "ice_update_flx()" at kt =', kt154 ! CALL ice_update_flx( kt ) ! -- Update ocean surface mass, heat and salt fluxes155 !IF(lwp) WRITE(numout,*) 'LOLO:ice_stp()@icestp.F90 => Done with "ice_update_flx()" at kt =', kt156 !IF(lwp) WRITE(numout,*) ''157 158 154 ! --- case we bypass ice thermodynamics --- ! 159 155 !IF( .NOT. ln_icethd ) THEN ! we suppose ice is impermeable => ocean is isolated from atmosphere … … 196 192 ! --- Ocean time step --- ! 197 193 !-------------------------! 198 !IF( ln_icedyn ) CALL ice_update_tau( kt, uu(:,:,1,Kbb), vv(:,:,1,Kbb) ) ! -- update surface ocean stresses199 !!gm remark, the ocean-ice stress is not saved in ice diag call above ..... find a solution!!!200 194 ! 201 195 IF( ln_timing ) CALL timing_stop('ice_stp') … … 238 232 ! 239 233 CALL ice_sbc_init ! set ice-ocean and ice-atm. coupling parameters 234 ! 235 fr_i (:,:) = at_i(:,:) ! initialisation of sea-ice fraction 240 236 ! 241 237 IF( ln_rstart ) CALL iom_close( numrir ) ! close input ice restart file
Note: See TracChangeset
for help on using the changeset viewer.