17 subroutine cu_ntiedtke_pre_run (flag_init, flag_restart, kdt, fhour, dtp, t, q, prevst, prevsq, &
18 forcet, forceq, errmsg, errflg)
24 logical,
intent(in) :: flag_init
25 logical,
intent(in) :: flag_restart
26 integer,
intent(in) :: kdt
27 real(kind_phys),
intent(in) :: fhour
28 real(kind_phys),
intent(in) :: dtp
29 real(kind_phys),
intent(in) :: t(:,:)
30 real(kind_phys),
intent(in) :: q(:,:)
31 real(kind_phys),
intent(in),
optional :: prevst(:,:)
32 real(kind_phys),
intent(in),
optional :: prevsq(:,:)
33 real(kind_phys),
intent(out),
optional :: forcet(:,:)
34 real(kind_phys),
intent(out),
optional :: forceq(:,:)
35 character(len=*),
intent(out) :: errmsg
36 integer,
intent(out) :: errflg
39 real(kind=kind_phys) :: dtdyn
48 if(flag_init .and. .not.flag_restart)
then
52 dtdyn=3600.0*(fhour)/kdt
54 forcet(:,:)=(t(:,:) - prevst(:,:))/dtp
55 forceq(:,:)=(q(:,:) - prevsq(:,:))/dtp
57 forcet(:,:)=(t(:,:) - prevst(:,:))/dtdyn
58 forceq(:,:)=(q(:,:) - prevsq(:,:))/dtdyn