Prepare script generating sic/sst data

ln -s /bdd/ERA5/NETCDF/GLOBAL_025/hourly/AN_SF/2024/ci.2024* .
ln -s /bdd/ERA5/NETCDF/GLOBAL_025/hourly/AN_SF/2024/sstk.2024* .

for file in ci.2024* ; do
    cdo -daymean $file day_$file
done
ncrcat day_ci.2024* siconc_ERA5_GLOBAL025_2024_daily.nc
# added 
ncrename -v ci,siconc ....nc
ncks --fix_rec_dmn time siconc_ERA5_GLOBAL025_2024_daily.nc siconc_ERA5_GLOBAL025_2024_daily_tfix.nc
ncap2 -s 'siconc=siconc*100.' siconc_ERA5_GLOBAL025_2024_daily_tfix.nc hadisst_siconc.nc
ncatted -a units,siconc,o,c,"%" hadisst_siconc.nc

for file in sstk.2024* ; do
    cdo -daymean $file day_$file
done
ncrcat day_sstk.2024* sstk_ERA5_GLOBAL025_2024_daily.nc
ncks --fix_rec_dmn time sstk_ERA5_GLOBAL025_2024_daily.nc sstk_ERA5_GLOBAL025_2024_daily_tfix.nc
ncap2 -s 'sstk=sstk-273.15' sstk_ERA5_GLOBAL025_2024_daily_tfix.nc hadisst_tos.nc
ncrename -v sstk,tos hadisst_tos.nc
ncatted -a units,tos,o,c,"C" hadisst_tos.nc

Normally, the version v2 is used which is converted to netcdf4: ncks -4 --no_rec_dmn time $dirold/${file} -o $dirnew/${file}

Some required modifications on the lmdz side:

<!-- Case with amip files -->
<!--  
      <file id="amipbc_sic_1x1" name="amipbc_sic_1x1" >
        <field id="sic_reg" name="sicbcs" domain_ref="domain_limit_amip" axis_ref="time_sic" operation="instant" freq_offset="1ts"/>
      </file>
 
      <file id="amipbc_sst_1x1" name="amipbc_sst_1x1" >
         <field id="sst_reg" name="tosbcs" domain_ref="domain_limit_amip" axis_ref="time_sst" operation="instant" freq_offset="1ts">this+273.15</field> 
      </file>
-->

<!-- Case with HadISST files -->
<!--  
      <file id="siconc" name="hadisst_siconc" >
        <field id="sic_reg" name="sicbcs" domain_ref="domain_limit_amip" axis_ref="time_sic" operation="instant" freq_offset="1ts"/>
      </file>
 
      <file id="tos" name="hadisst_tos" >
         <field id="sst_reg" name="tosbcs" domain_ref="domain_limit_amip" axis_ref="time_sst" operation="instant" freq_offset="1ts">this+273.15</field> 
      </file>
-->

<!-- Case with ERA5 files -->
      <file id="ci" name="ci" >
        <field id="sic_reg" name="siconc" domain_ref="domain_limit_amip" axis_ref="time_sic" operation="instant" freq_offset="1ts"/>
      </file>

      <file id="sstk" name="sstk" >
         <field id="sst_reg" name="sstk" domain_ref="domain_limit_amip" axis_ref="time_sst" operation="instant" freq_offset="1ts"/>
      </file>