From 526870458f5deff4e92b3aa8b6e5aad28f965806 Mon Sep 17 00:00:00 2001 From: Lorenzo Figini Date: Fri, 10 Feb 2017 10:41:24 +0000 Subject: [PATCH] fixed rhot print in central ray output (fort.4) --- src/gray_params.f90 | 4 ++-- src/graycore.f90 | 2 +- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/src/gray_params.f90 b/src/gray_params.f90 index caec466..0407d4f 100644 --- a/src/gray_params.f90 +++ b/src/gray_params.f90 @@ -232,8 +232,8 @@ contains read(u,*) prfparam%iprof, prfparam%irho ! irho=0,1,2 -> num profiles vs rhot,rhop,psin read(u,*) prfparam%filenm ! psbnd value of psi ( > 1 ) of density boundary - read(u,*) prfparam%psnbnd !, prfparam%sspld - prfparam%sspld=0.001_wp_ + read(u,*) prfparam%psnbnd, prfparam%sspld +! prfparam%sspld=0.001_wp_ ! iscal ne Te scaling 0: nustar=const, 1: n_greenw=const; 2 no rescaling ! factT factn factor for Te&ne scaling read(u,*) prfparam%factte, prfparam%factne, prfparam%iscal diff --git a/src/graycore.f90 b/src/graycore.f90 index 7fcaec1..022fc84 100644 --- a/src/graycore.f90 +++ b/src/graycore.f90 @@ -1959,7 +1959,7 @@ bb: do if(jk.eq.1) then phideg=atan2(yym,xxm)/degree if(psinv>=zero .and. psinv<=one) then - rhot=frhotor(psinv) + rhot=frhotor(sqrt(psinv)) else rhot=1.0_wp_ end if