search

ShunLiu-NOAA / regional_workflow

Other
0 stars 0 forks source link

FV3LAMDAX without LBC expansion #54

Open ShunLiu-NOAA opened 2 years ago

ShunLiu-NOAA commented 2 years ago

document the change of removing LBC expansion fv3lamdax. It will be inverse step to add this to RRFS.

workflow changes: I. in exregional_gsianl.sh, do not copy "_new" file. !This file contains horizontal grid information !Now reads file with enlarged grid specs cp $fv3_case/grid_spec_new.nc fv3_grid_spec

! This file contains 3d fields u,v,w,dz,T,delp, and 2d sfc geopotential phis ! Now reads enlarged file cp $fv3_case/fv_core.res.tile1_new.nc fv3_dynvars

! This file contains 3d tracer fields sphum, liq_wat, o3mr ! Now reads enlarged file cp $fv3_case/fv_tracer.res.tile1_new.nc fv3_tracer

! This file contains surface fields (vert dims of 3, 4, and 63) ! Now reads enlarged file cp $fv3_case/sfc_data_new.nc fv3_sfcdata

  1. in input.nml, regional_bcs_from_gsi = .false. write_restart_with_bcs = .false.

  2. turn off "regional_gsibcprep" in drive_fv3sar_da.xml

  3. ! Put remaining analysis files in ANLdir (defined in J-job) mv fv3_akbk $ANLdir/fv_core.res.nc mv coupler.res $ANLdir/coupler.res

    mv fv3_sfcdata $ANLdir/sfc_data.nc

    cp $COMOUT/anl.tm06/gfs_ctrl.nc $ANLdir/.

  4. use ncat in exregional_gsianl.sh !run ncatted to fix checksum aborts in the new code ncatted -a checksum,,d,, fv_core.res.tile1.nc fv_core.res.tile1.nc_save ncatted -a checksum,,d,, fv_tracer.res.tile1.nc fv_tracer.res.tile1.nc_save rm fv_core.res.tile1.nc fv_tracer.res.tile1.nc mv fv_core.res.tile1.nc_save $ANLdir/fv_core.res.tile1.nc mv fv_tracer.res.tile1.nc_save $ANLdir/fv_tracer.res.tile1.nc

  5. change the number of boundary files in exregional_forecast.sh numbndy_check=$NHRS+2 to numbndy_check=$NHRS+1 change if [ $numbndy -ne 3 ] ; to if [ $numbndy -ne 2 ] ;

ShunLiu-NOAA commented 2 years ago

Evaluation result: https://docs.google.com/presentation/d/1zRaloxVQb5W7XNM2gCTMXVCs4Nt4SEW5x4Mau4LoAnQ/edit?usp=sharing

ShunLiu-NOAA commented 2 years ago

EMC parallel is using the function of LBC expansion. The function is not in RRFS prototype yet. Let's discuss LBC related issues in DA and discuss if we should use this in RRFS. Eric added LBC expansion to EMC's parallel, I have a comparison test. It is attached here.

https://docs.google.com/presentation/d/1zRaloxVQb5W7XNM2gCTMXVCs4Nt4SEW5x4Mau4LoAnQ/edit?usp=sharing

Eric will present: add domain boundary rows to regional FV3 DA https://docs.google.com/presentation/d/1AR4E2B_-EWEcJVRsuDqJ-Z5u3VquVuzlOrbbA0aAcU4/edit?usp=sharing