Strona główna Odkrywaj Pomoc
Zaloguj się
pelletie
/
PARAMOUR_TOOLS
2
0
Forkuj 0
Pliki Problemy 0 Oczekujące zmiany 0 Wiki
Drzewo: 1fcf9f86a1
Gałęzie Tagi
PARAMOUR_TOOLS
/
anom_bias
/
calculate_biases.sh

calculate_biases.sh 3.1 KB
Historia Czysty

12345678910111213141516171819202122232425262728293031323334353637383940 
  1. #!/bin/bash
    2. 

  2. set -u
    3. 

  3. 

  4. # This script takes ORAS5 and EC-Earth historical ocean temperature and salinity conditions and calculates biases from them.
    5. 

  5. # Author: D. Verfaillie
    6. 

  6. 

  7. # Inputs
    8. 

  8. in_folder="/storepelican/dverfail/bsc/anom_bias"
    9. 

  9. oras5_file="ORAS5_preinterp_opa1_grid_T_1985-2014"
    10. 

  10. histo_1m_thetao="thetao_Omon_EC-Earth3_historical_ensmean_gn_198501-201412"
    11. 

  11. histo_1m_so="so_Omon_EC-Earth3_historical_ensmean_gn_198501-201412"
    12. 

  12. 

  13. # Prepare ORAS5 file to have the same grid and format as the EC-Earth files
    14. 

  14. module --force purge
    15. 

  15. module load releases/2020b
    16. 

  16. module load noarch
    17. 

  17. module load NCO
    18. 

  18. module load CDO
    19. 

  19. 

  20. cd ${in_folder}/"oras5_clim"
    21. 

  21. 

  22. cdo remapbil,${in_folder}/oras5_clim/ECEarth_grid_T_new.grid ${oras5_file}.nc ${oras5_file}_remapped.nc
    23. 

  23. ncrename -O -d time_counter,time -d x,i -d y,j -d nvertex,vertices -d deptht,lev -v time_counter,time -v nav_lon,longitude -v nav_lat,latitude -v nav_lon_bnds,longitude_bnds -v nav_lat_bnds,latitude_bnds -v deptht,lev -v conservative_temperature,thetao -v absolute_salinity,so ${oras5_file}_remapped.nc
    24. 

  24. ncap2 -O -s 'defdim("bnds",2);make_bounds(lev,$bnds,"lev_bnds");' ${oras5_file}_remapped.nc ${oras5_file}_remapped_bounds.nc
    25. 

  25. ncks -A -v i,j,lev_bnds ../histo_clim/${histo_1m_thetao}.nc ${oras5_file}_remapped_bounds.nc
    26. 

  26. ncap2 -s 'lev=double(lev); lev_bnds=double(lev_bnds); longitude=float(longitude); longitude_bnds=float(longitude_bnds); latitude=float(latitude); latitude_bnds=float(latitude_bnds);' ${oras5_file}_remapped_bounds.nc ${oras5_file}_remapped_bounds_reformatted.nc
    27. 

  27. ncks -O -x -v so,sossheig ${oras5_file}_remapped_bounds_reformatted.nc thetao_ORAS5_preinterp_opa1_grid_T_1985-2014.nc
    28. 

  28. ncks -O -x -v thetao,sossheig ${oras5_file}_remapped_bounds_reformatted.nc so_ORAS5_preinterp_opa1_grid_T_1985-2014.nc
    29. 

  29. 

  30. # Calculate monthly and climatological biases: histo - oras5 (these biases will then have to be SUBSTRACTED from the hindcasts)
    31. 

  31. ncdiff ../histo_clim/${histo_1m_so}.nc thetao_ORAS5_preinterp_opa1_grid_T_1985-2014.nc ../bias_histo-oras5/thetao_bias_1985-2014_1m.nc
    32. 

  32. ncdiff ../histo_clim/${histo_1m_thetao}.nc so_ORAS5_preinterp_opa1_grid_T_1985-2014.nc ../bias_histo-oras5/so_bias_1985-2014_1m.nc
    33. 

  33. 

  34. cd ${in_folder}/bias_histo-oras5
    35. 

  35. cdo ymonmean thetao_bias_1985-2014_1m.nc thetao_bias_1985-2014_Mm.nc
    36. 

  36. cdo ymonmean so_bias_1985-2014_1m.nc so_bias_1985-2014_Mm.nc
    37. 

  37. 

  38. # Reformat the bias files to have same format as EC-Earth hindcast ocean temperature and salinity files
    39. 

  39. ncrename -O -d time,time_counter -d i,x -d j,y -d vertices,nvertex -d lev,olevel -d bnds,axis_nbounds -v time,time_counter -v time_bnds,time_counter_bounds -v latitude,nav_lat -v latitude_bnds,bounds_nav_lat -v longitude,nav_lon -v longitude_bnds,bounds_nav_lon -v lev,olevel -v lev_bnds,olevel_bounds thetao_bias_1985-2014_Mm.nc thetao_bias_1985-2014_Mm_reformatted.nc
    40. 

  40. ncrename -O -d time,time_counter -d i,x -d j,y -d vertices,nvertex -d lev,olevel -d bnds,axis_nbounds -v time,time_counter -v time_bnds,time_counter_bounds -v latitude,nav_lat -v latitude_bnds,bounds_nav_lat -v longitude,nav_lon -v longitude_bnds,bounds_nav_lon -v lev,olevel -v lev_bnds,olevel_bounds so_bias_1985-2014_Mm.nc so_bias_1985-2014_Mm_reformatted.nc
    41.