Running Roms during tropical Cyclone

[prakash22]

Dear all

I am running roms model for tropical cyclone case

I am using 3 hourly forcing from wrf-arw output(wind-stress, pressure, netheatflux, fresh-water flux,rain cloud).
I define my hedder file (.h) as follow:
/*==========Grid Configuration========*/
# define CURVGRID
# define MASKING


/*===========Model Physics=========*/
# define SOLVE3D
# define UV_COR
# define UV_VIS2
# define TS_DIF2
# define UV_ADV
# define LMD_MIXING
# define AVERAGES


/*# define BULK_FLUXES /* Bulk fluxes computation turn OFF */
/*# define EMINUSP /* E-P computation */
/*# undef EMINUSP_SSH
# define ALBEDO /* Albedo equation used for shortwave radiation turn ON */
/*# define DIURNAL_SRFLUX /* Shortwave radiation imposed local diurnal cycle turn OFF */
# define CLOUDS
# define AIR_TEMP
# define ATM_PRESS
/*# define CCS_FLUXES
#define ANA_SMFLUX */
/* Equation of State */
# define SALINITY
# define NONLIN_EOS

/*Numerical Scheme options */

# define TS_U3HADVECTION
# define TS_SVADVECTION
# define DJ_GRADPS
# define MIX_S_UV
# define MIX_S_TS
# define SPLINES

# define MIX_S_TS
/*Vertical Mixing*/
# ifdef LMD_MIXING
# define DIURNAL_SRFLUX
# define SOLAR_SOURCE
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_NONLOCAL
# define LMD_SKPP
# define LMD_BKPP
# endif


/* Surface Boundary Condition or Forcing */
/*# define QCORRECTION
# define SCORRECTION */
/*#define ANA_SMFLUX /* use if analytical surface momentum stress */

/*Lateral Forcing*/
# define CLIMATOLOGY
# ifdef CLIMATOLOGY
# define ZCLIMATOLOGY
# define M3CLIMATOLOGY
# define TCLIMATOLOGY
# endif
/* Bottom Boundary Condition */
# define ANA_BSFLUX
# define ANA_BTFLUX
# define UV_LDRAG
/* Open Boundary Conditions */

# define M3CLM_NUDGING
# define TCLM_NUDGING

/* Input/Output & Diagnostics */
# define UV_LDRAG
# define ANA_FSOBC
# define ANA_M2OBC

/* Lateral Boundary Conditions */
# define WESTERN_WALL
# define NORTHERN_WALL
# undef SOUTHERN_WALL
# undef EASTERN_WALL

# undef NO_LBC_ATT

my question is that, when I am running roms for 4 days cyclone upwelling is happening alone the track before land fall but is should not be happen.
I think that it is happening because of strong current.

Could any one help me

Thanks & regard

Prakash

[kate]

First of all, make sure you are defining what you hope to be defining:

Code: Select all

/*# undef EMINUSP_SSH
# define ALBEDO /* Albedo equation used for shortwave radiation turn ON */

Here, ALBEDO is not defined because C-style comments go from "/*" to "*/". The same is true for SCORRECTION here:

Code: Select all

/*# define QCORRECTION
# define SCORRECTION */

Maybe that's what you want, but I just want to check. Most of us would #undef for clarity, or just delete lines we know we never want.

Second, to check the default (non-storm) response of the model, I would run for some other time.

So, how is it getting the heat flux? You don't have ANA_STFLUX listed, so I assume the model wants to read in a surface heat flux. If so, what is the purpose of providing the air temperature and pressure?

Code: Select all

# define WESTERN_WALL

Is this a really old version or just an old option that ROMS will ignore?

[prakash22]

Hello All Members,

I am running ROMS model for the cyclone case study by using the ECMWF surface forcing data and WRF forcing data sets without activating BULK_FORMULA_PHYSICS, but results appear to me that whatever wind stress I am providing to the model, currents are also showing the same intensity in the model domain that is not the expected result. I tried to change the mixing coefficient values in the include (.in) file as very high, very low, even when I commented also. But the model simulations remain same.


! Harmonic/biharmonic horizontal diffusion of tracer for nonlinear model
! and adjoint-based algorithms: [1:NAT+NPT,Ngrids].

TNU2 == 150.0d0 150.0d0 ! m2/s
TNU4 == 2*500.0d0 ! 2.0d0 ! m4/s

ad_TNU2 == 0.0d0 0.0d0 ! m2/s
ad_TNU4 == 0.0d0 0.0d0 ! m4/s

! Harmonic/biharmonic, horizontal viscosity coefficient for nonlinear model
! and adjoint-based algorithms: [Ngrids].

VISC2 == 100.0d0 ! m2/s
VISC4 == 50.0d0 ! m4/s

ad_VISC2 == 0.0d0 ! m2/s
ad_VISC4 == 0.0d0 ! m4/s

! Vertical mixing coefficients for tracers in nonlinear model and
! basic state scale factor in adjoint-based algorithms: [1:NAT+NPT,Ngrids]

AKT_BAK == 5.0d-6 5.0d-6 ! m2/s

ad_AKT_fac == 1.0d0 1.0d0 ! nondimensional

! Vertical mixing coefficient for momentum for nonlinear model and
! basic state scale factor in adjoint-based algorithms: [Ngrids].

AKV_BAK == 5.0d-5 ! m2/s

ad_AKV_fac == 1.0d0 ! nondimensional




I am very surprised whether model is not able to take these values of mixing and viscosity coefficients either I have wrongly defined some contradictory physics as CPP option in header file or need to define some appropriate OPTION to activate this.


My question is if my currents values are high and causing mare cooling in comparison to observations, how can I control this behaviour?



Next is I am adding my header file options for activated physics for this experiment, can anyone suggest me where I am choosing the wrong physics or I should add more physics (CPP options) to deal with this issue?






/*==========Grid Configuration========*/
# define CURVGRID
# define MASKING


/*===========Model Physics=========*/
# define SOLVE3D
# define UV_COR
# define UV_VIS2
# define TS_DIF2
# define UV_C4ADVECTION
# define UV_VIS4
# define UV_ADV
# define LMD_MIXING
# define AVERAGES



# define CLOUDS
# define ATM_PRESS
/* Equation of State */
# define SALINITY
# define NONLIN_EOS

/*Numerical Scheme options */

# define TS_U3HADVECTION
# define TS_SVADVECTION
# define DJ_GRADPS
# define MIX_S_UV
# define MIX_S_TS
# define SPLINES

# define MIX_S_TS
/*Vertical Mixing*/
# ifdef LMD_MIXING
# define DIURNAL_SRFLUX
# define SOLAR_SOURCE
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_NONLOCAL
# define LMD_SKPP
# define LMD_BKPP
# endif


/*Lateral Forcing*/
# define CLIMATOLOGY
# ifdef CLIMATOLOGY
# define ZCLIMATOLOGY
# define M3CLIMATOLOGY
# define TCLIMATOLOGY
# endif
/* Bottom Boundary Condition */
# define ANA_BSFLUX
# define ANA_BTFLUX
# define UV_LDRAG
/* Open Boundary Conditions */

# define M3CLM_NUDGING
# define TCLM_NUDGING

/* Input/Output & Diagnostics */
# define UV_LDRAG
# define ANA_FSOBC
# define ANA_M2OBC

/* Lateral Boundary Conditions */
# define WESTERN_WALL
# define NORTHERN_WALL
# undef SOUTHERN_WALL
# undef EASTERN_WALL

# undef NO_LBC_ATT

With Regard
Prakash