Problems with npzd_Powell

[gabrielanalini]

Hi!

I am studying the response of phytoplankton community to eddy activity in Brazilian Southeast coast.

My hydrodinamic model is working fine, the problem happens with the npzd_Powell model. I made tests with both homogeneus initial conditions and idealized npzd fields. There seems to be no response or interaction between the npzd compartments.The vertical structure for physical variables works, but for the biological parameters not.

I made a run using the upwelling testcase with npzd_Powell active and everything seems to be fine, with instant response from the compartments.

Is there something else to set in the CP or .in for the coupling works?

I'm sending the .h and .in images and some plots.

I’d appreciate any help with this.

Thanks in advance!

[kate]

First, you should find out why the colorbar for the first two vertical slice plots wants to extend to 1.e37. Are those points in the land mask?

[wilkin]

Your log file shows the following:

Metrics information for Grid 01:
===============================

Minimum X-grid spacing, DXmin = 0.00000000E+00 km
Maximum X-grid spacing, DXmax = 2.97235216E+00 km
Minimum Y-grid spacing, DYmin = 0.00000000E+00 km
Maximum Y-grid spacing, DYmax = 2.97199962E+00 km
Minimum Z-grid spacing, DZmin = 0.00000000E+00 m
Maximum Z-grid spacing, DZmax = 3.91620893E+02 m

That's odd that you have zero as the minimum grid spacings in all dimensions.

Have you applied special/missing values in the land mask points of the grid metrics? I don't think ROMS will like that.

[gabrielanalini]

Thank you very much for the replies Kate and Wilkin.

Kate, there is no land mask on the biological initial fields, I'm working on this.


Wilkin, I don't know why this is happening. I used to run my experiments in another machine, using MPI, and the minimum grid-spacing was fine. Now, im using OpenMP in a new machine, but I know people who use OpenMP as well without problems. This could be a problem?

[chaoyu1985]

Thank you very much for the replies Kate and Wilkin.

Kate, there is no land mask on the biological initial fields, I'm working on this.


Wilkin, I don't know why this is happening. I used to run my experiments in another machine, using MPI, and the minimum grid-spacing was fine. Now, im using OpenMP in a new machine, but I know people who use OpenMP as well without problems. This could be a problem?

Hi, i am runing npzd as well. How did u give the boundary condition (i.e.alkalinity,NH4,SDeN). In my study area,NH4 data is absent. I can not find any other reanalysis data from internet. So,I used a constent value for all grids and boundary condition. the model can not run sucessuflly. How did you set the boundary condition in the time series. did they vary with time ? did they match with the changes of the chlorophyll and phytoplankton? thank you so much for your help.