rosimarwx / MPAS-A_aquaplanet

Modified MPAS-Atmosphere with aquaplanet capability
Other
1 stars 0 forks source link

What is the format of sounding.input #1

Open Rock-Dai opened 2 years ago

Rock-Dai commented 2 years ago

Hello,

I am trying to using your aquaplanet model. I have been run into a problem of the lack of sounding.input. I have see the reference and noticed that you used the sounding at every grid point. I wonder the format of the sounding.input in this model, and then I can generate by myself.

rosimarwx commented 2 years ago

Thanks so much for bringing this to my attention! I will add the sounding.input file to the repository. In the meantime, here is an example sounding.input file. I also include apeo3.dat, which you will need to run the model with symmetric ozone.

On Fri, Mar 18, 2022 at 3:52 AM Rock-Dai @.***> wrote:

Hello,

I am trying to using your aquaplanet model. I have been run into a problem of the lack of sounding.input. I have see the reference and noticed that you used the sounding at every grid point. I wonder the format of the sounding.input in this model, and then I can generate by myself.

— Reply to this email directly, view it on GitHub https://github.com/rosimarwx/MPAS-A_aquaplanet/issues/1, or unsubscribe https://github.com/notifications/unsubscribe-auth/AHTDMZ6EI2M2RXAJCLCEEKTVARG4FANCNFSM5RBK6UAQ . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you are subscribed to this thread.Message ID: @.***>

--


Rosimar Rios-Berrios, PhD (she/her/hers) Scientist I, Mesoscale and Microscale Meteorology National Center for Atmospheric Research https://rosimarwx.github.io/


rosimarwx commented 2 years ago

I have updated the repository with an example sounding.input file. I also added example namelists per an offline request. Let me know if this solves your issues.

Rock-Dai commented 2 years ago

Thanks for providing the input data, and it really worked and helped me finish the run of initial data(init_atmosphere_model). However, I meet some trouble in the model run(atmosphere_model), and I find the reason is the wind speed calculated in boundary layer module. The wind speed equals to 0 would lead ross to 0, too, and the model would throw a divided by 0 and stop when calculated the brcr_sbro. I found that the code is from the WRF code base. Did you have the same problem before? Do you think this problem is caused by my wrong configuration or the WRF code itself is wrong?

/src/core_atmosphere/physics/physics_wrf/module_bl_ysu.F:938-945 do i = its,ite if((.not.stable(i)).and.((xland(i)-1.5).ge.0))then wspd10 = u10(i)u10(i) + v10(i)v10(i) wspd10 = sqrt(wspd10) ross = wspd10 / (coriznt(i)) brcr_sbro(i) = min(0.16(1.e-7*ross)**(-0.18),.3) endif enddo

BTW, may I know the compiler you used for compiling?

Rock-Dai commented 2 years ago

u_error_v2

Fig.1: U in nVertLevels 1

Horizontal normal velocity

Fig.2: Horizontal Normal Velocity

Hello, sorry for bothering again. I have successfully run the model and have a 2-month output. However, I have a severe problem about the data convert. I read the tutorial of MPAS from NCAR online, and use their program to convert the original history file to the lat-lon orthogonal grid output. I have successfully converted several variables. However, when I converted the u, v data, I found that there is a strange sub-grid scale pattern. I display the pattern in fig.1, and I also use NCL to plot in small region to show the detailed pattern. I know that the u, v vector data is set on the edges of hexagonal cell, which is different from other variables that set in the center of the cell, and this may be the key to the problem. After several try, I still can't solve the problem. Have you ever met this problem before? Did I omit some key process in data convert?

Tutorial Link: https://www2.mmm.ucar.edu/projects/mpas/tutorial/Boulder2019/index.html
Convert Program: https://github.com/mgduda/convert_mpas