WRF Timeseries Output Inverse Density 'alt'

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

J

Jslaten

New member
Hello,

I am currently attempting at extracting the instantaneous inverse density as a timeseries with the TSLIST option. I have been following the guide here :

https://forum.mmm.ucar.edu/phpBB3/viewtopic.php?t=8898

I believe I may be missing something in my own version with this error during compilation ,
================================================================
Code: 
wrf_timeseries.f90:38:74:
       ts_file_endings = (/ 'UU', 'VV', 'PH', 'TH', 'QV' ,'WW', 'PR', 'QP', 'RHO' /)
                                                                          1
Error: Different CHARACTER lengths (2/3) in array constructor at (1)
wrf_timeseries.f90:236:58:

                   WRITE(ts_profile_filename(i-1:i),'(A2)') ts_file_endings(j)
                                                          1
Error: Function ‘ts_file_endings’ at (1) has no IMPLICIT type
wrf_timeseries.f90:228:42:

                DO j=1,SIZE(ts_file_endings)
                                          1
====================================================

Here QP is a new variable we have previously added with no problem so likely that isn't the error. I'm not really comfortable with FORTRAN but I believe the error comes from this line that is added in the #EM_CORE block,

Code: 
 grid%ts_rho_profile(n,i,k) = grid%alt(ix,k,iy)


I can provide any additional changes I've made if that helps spot the problem.


Thank You,
Jesse
 
Update 18 Jan 2021

A colleague of mine was able to diagnose the issue I was having prior, as it turns out the ts_file_endings is expecting only characters of length 2 where as my variable name was 'RHO' it should have read 'RH'. The correct line is as such

Code: 
      ts_file_endings = (/ 'UU', 'VV', 'PH', 'TH', 'QV' ,'WW', 'PR', 'QP', 'RH' /)

Also for ease of life we can store the value of 'alt' density as follows

Code: 
             grid%ts_rho_profile(n,i,k) = 1/grid%alt(ix,k,iy)

After this fix I am able to compile ... Sadly my new issue is the associated CSV files are not saving during runtime :^)


* Update 19 Jan 2021

I found the solution to my problem but it's rather specific to my own WRF build so I'm not sure how much it would help others.



- Jesse
 
Jesse,
Thanks for the detailed description of the problem and solution. We appreciate your kindness to share this information with the community.
 
You must log in or register to reply here.
Back
Top